Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,565 G→T V218V (GTG→GTT radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,5650GT96.8% 174.1 / ‑5.3 63V218V (GTG→GTTradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (26/35);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGC  >  NC_000913/4626324‑4626805
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
cAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCCCTTCGCCAAAACGCGGGGGGTGGCGGTTGTCTTGGTTGGGCACGtt                                                                                                                                                                                                                                             >  2:404747/1‑250 (MQ=255)
  aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTATCATGGTTGGACACGTAAc                                                                                                                                                                                                                                           <  1:321269/251‑1 (MQ=255)
     cgcgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCa                                                                                                                                                                                                                                         >  1:129521/1‑250 (MQ=255)
              cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAAGACGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAAAGTCATCGCCTGGCAGCGTGGCGCAGGTGAGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAAAGCGCGGTTGGGCGATTGTCATGGGTAGACACGTAACCCAAGAttgctc                                                                                                                                                                                                                               >  1:460729/1‑246 (MQ=255)
                          ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGGTGGACACGTAACCAAAGATGGCTCGCTGGCGGGccc                                                                                                                                                                                                                   >  1:359096/1‑251 (MQ=255)
                                        cccgtccgatatggAAACCAGCATCGTGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGAGGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAac                                                                                                                                                                                                    <  2:200376/239‑1 (MQ=255)
                                                 gtCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATc                                                                                                                                                                                           <  1:111345/250‑1 (MQ=255)
                                                    ggAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTTAAGGGCTGGAacagaacgggag                                                                                                                                                                                          >  1:200376/1‑240 (MQ=255)
                                                    ggAAACCAGCATCGAGCACATCTGCCTGTGTGCAGAAGAAGAGCAACCGAAGCTGAGGGTAATTGACGCGACCCAGTGGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTTAAACTGCGGCTTATCTTACGCGCTTCCCCAAAACGCGCGGTGTGGCGAGTGTCATGGTTGGACTCGTAACCAAAGATGGCTCGCTGGCTGGCCCTAAAGTGCTGGAACACTGTATCgcct                                                                                                                                                                                       <  2:425107/251‑4 (MQ=255)
                                                            gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                                                                                                                                               <  1:338391/251‑1 (MQ=255)
                                                              aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCAAATGGCGGATGGACAGTCATAGCCTGGCAGCGTGGCGCAGGGGCGGGAAACGGAGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGAGGGGCGATTGTCATGGGTGGAAACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGAGCTGGAAAACTGTATCGACAGTTCGGGGCt                                                                                                                                                                             >  1:176099/1‑251 (MQ=255)
                                                                      gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                             <  2:235181/227‑1 (MQ=255)
                                                                      gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                             >  1:235181/1‑227 (MQ=255)
                                                                           gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGTTCGGTGCTTTTGGATGGCg                                                                                                                                                                 >  1:76635/1‑249 (MQ=255)
                                                                              agATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGcc                                                                                                                                                            <  1:40660/250‑1 (MQ=255)
                                                                                aTTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTAACGCATGTTCGGTGCTTTTGGATGGCGATGCCgg                                                                                                                                                          >  1:347776/1‑250 (MQ=255)
                                                                                           agagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                                                                                                                                               <  1:148646/251‑1 (MQ=255)
                                                                                            gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                                                              <  1:447822/251‑1 (MQ=255)
                                                                                              ccaccgtaaGCTGATGGTAATTGACGCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTACCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCa                                                                                                                                            <  1:127369/244‑1 (MQ=255)
                                                                                                  ccGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGAATCCCGGTTTCGCActtt                                                                                                                                        >  1:423422/1‑248 (MQ=255)
                                                                                                    gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTTGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTGCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActttac                                                                                                                                      <  1:104364/251‑6 (MQ=255)
                                                                                                      aGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCCGGCGGGGCCCAAAATGCTGGAAAACTGTATCGCCTTTTTGGTGCTTTTTGGTGGCGATGCCGACTCCCCTTTTCGCActtttcgc                                                                                                                                    >  2:144475/1‑244 (MQ=255)
                                                                                                          gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGcc                                                                                                                                <  1:380788/251‑1 (MQ=255)
                                                                                                               tAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCATGGTTGGCCACGTAACCAAAGATGGCTCGCGGGCGGGCCCTAAAGTGCTGGAACACTGTATAGACTGTTTGGGGGGTTTGGGTGGGGATGCCGACTCCCCtttttgccctttaccccccccaca                                                                                                                            >  2:9160/1‑230 (MQ=255)
                                                                                                                    gACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                                                                                                       <  2:305367/250‑1 (MQ=255)
                                                                                                                                gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                                                                          <  1:298075/251‑1 (MQ=255)
                                                                                                                                gTGATGCATATGGCGGATGTACAGACATCGCCTGGCAGCGGGGAGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCCCGCTGGCGGGCCCTTAAGTGCTGGGACACTGGTACGCATGTTCGGTGCCTTTGGAAGGGCGATCAGCCTCCCGTTTTCCCAACTTACGCCGCCATAAAAACCCCttaggggaggg                                                                                                          >  2:396544/1‑242 (MQ=255)
                                                                                                                                   aTGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGaa                                                                                                       <  1:35723/251‑1 (MQ=255)
                                                                                                                                     gCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATg                                                                                                     <  2:372209/251‑1 (MQ=255)
                                                                                                                                       atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGGGCTATTGAATGGCGTTGCCGACTCCCGGTTTCGCACTTTACGCAGCCATAAAAACCGCTTCAGGGCGGTGAATg                                                                                                     >  1:457861/1‑249 (MQ=255)
                                                                                                                                         atGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGGGGGCCCTAAAGAGGGGGAGCACTGTTTCGACTGGTCGGGGCGTTTGGATGGGGAGGCCGGCTCCCGTTTTTGCAATTTACGCCGCCCTAAAAACCCCTTCgggccgggggatgggat                                                                                                 >  2:292824/1‑236 (MQ=255)
                                                                                                                                             cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGGGGGCGGTGCCGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCggctcgggcgggggaatgggctgggg                                                                                             >  2:45882/1‑226 (MQ=255)
                                                                                                                                                    accgTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAACTGACTCGCTTCTCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGTTCGCTGGCGGGCCCTAAAGTTCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATGAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcg                                                                                      <  2:176099/248‑1 (MQ=255)
                                                                                                                                                               cTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGa                                                                           >  1:51192/1‑250 (MQ=255)
                                                                                                                                                                      gggggCGCTGGGGCGTGAAACGGTGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTTGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGc                                                                   <  2:16368/248‑1 (MQ=255)
                                                                                                                                                                        tGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATAGTCATGGATGGACACGTAACCAAAGATGGGTCGCTGGCGGGCCCTATAGTGCTGGACCACTGTATCGACTGGTCGGGGCTTTTGGATGGCGATGCCTACTCCCCGTTTCTCACTTTTCGCAGCCATAAAAAACGCTTCGGCGCGGTGGATGAGACGTGCGACTTCGCCTAGGACCGAACAgggctg                                                                  >  2:293936/1‑248 (MQ=255)
                                                                                                                                                                            gcAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                                                             >  1:137365/1‑251 (MQ=255)
                                                                                                                                                                               ggTGCATGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAg                                                          <  1:9771/251‑1 (MQ=255)
                                                                                                                                                                                   cGTGAAACGGCGGCTTATCGTACGCGCTTCGCCAAATCGCGCGGTGTGTCTATTTTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTGCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCTGCCATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAg                                                      <  1:38551/251‑1 (MQ=255)
                                                                                                                                                                                     tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGACCACTGTATCGCCTGTTCGGTGCTTTTGGGGGGCGATGCCGACTCCCGTTTTCGCCCTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGGCTGAGCTGGGCGTCTTCGCGATGCCCGAACAGGGGCTGCGTggaggcaac                                                     >  2:359403/1‑242 (MQ=255)
                                                                                                                                                                                     tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                       >  1:197324/1‑248 (MQ=255)
                                                                                                                                                                                     tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                       <  2:197324/248‑1 (MQ=255)
                                                                                                                                                                                           ggcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGTAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCtt                                              <  1:175006/251‑1 (MQ=255)
                                                                                                                                                                                            gtggcttacctGACGCGCCTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTTGACACGTAACCAAAGATGGCGCGCTGGCGTGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTTAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                             <  2:293833/241‑1 (MQ=255)
                                                                                                                                                                                            gcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                             <  1:200037/251‑1 (MQ=255)
                                                                                                                                                                                              ggtttagctGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTGGCGGGCGGGCCCTAAAGTGCTGGAACACGGTATCGACTGTGCGGTGCTTTTTGATGGCGATGCCGACTCCCGTTGTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCgg                                           <  2:162350/243‑1 (MQ=255)
                                                                                                                                                                                                       tacgcgattCGCCAAAACTCGCGGTGGGGCGATTGTCATGGGTGGATACGGAACCAAAGATGGTCCGCTGGTGGGCCGTAAAGTGCTGGACCACTGTATCGACTGTACGTTGTTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGTCCGAACAGGGTCTGCGTGAATACAGCAATATTTCGGCAATTTTc                                   <  2:394924/243‑1 (MQ=255)
                                                                                                                                                                                                       gacgcGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtt                                      <  1:241998/243‑1 (MQ=255)
                                                                                                                                                                                                          gcgcgtCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGGTGGGGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTaa                               <  1:61136/246‑1 (MQ=255)
                                                                                                                                                                                                           cGCTTCGCAAAAACGCGGGGTGTGGCGATTGTCATGGTGGGGCGCGTACCCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACGCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAg                              <  2:114656/251‑1 (MQ=255)
                                                                                                                                                                                                             cTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGCTCGGTGCCTTTGGATGGCGCTGCCGACTACCGTTTTCGCCATTTACGCACACATAAAAACCACTTCGGCGCGGGTAATGAGCAGGACTGCTTCGCGATAGCCCACCCGGGGCTGCGTGATGTAAGACACACTTAGGCAATTTTCTTTAGTc                            >  2:405223/1‑251 (MQ=255)
                                                                                                                                                                                                               tCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTc                            <  1:144475/249‑1 (MQ=255)
                                                                                                                                                                                                                  ccAAAACGCGCGGTGTGGCGATTGTAATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgga                       <  2:331747/251‑2 (MQ=255)
                                                                                                                                                                                                                    aaaaCGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGCTGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCCCGATGGCCGAACAGGGGCTCCGTGAAGTCAGCAACCCTTCGGACATTTTTTTAAGTCgcgggg                      >  2:389398/1‑248 (MQ=255)
                                                                                                                                                                                                                         ggggggTGTGGGGATTGTCATGGTTGGACACGGAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                <  2:9683/247‑1 (MQ=255)
                                                                                                                                                                                                                         gggcgGTGTGGCGATTGTCAGGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCTCTAAAGTGCTGGAACACTGTATCGACGGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTGCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                <  2:469532/249‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAgg               >  1:413023/1‑250 (MQ=255)
                                                                                                                                                                                                                          cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGGACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGt               <  1:32672/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGAc            <  1:221702/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAAACGCTTCGGCGCGGTGAATGGGCTGGGCGTCTTCGCGAAGACCGAAAAGGGGCTGCGGGACGTCAGCAACCCTTCGGCAATTTTCTTAAGGCgcggagtgggagggccctccggccg   >  2:414251/1‑230 (MQ=255)
                                                                                                                                                                                                                                        tGTCATGGTTGGACACGTAACCAACGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTGTCTTAATTCGCGGAGATGAAGTGACCTCCGGCAGc  <  1:483537/250‑1 (MQ=255)
                                                                                                                                                                                                                                        tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa             <  2:102891/239‑1 (MQ=255)
                                                                                                                                                                                                                                        tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa             >  1:102891/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
CAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGC  >  NC_000913/4626324‑4626805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: