Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,324,959 C→T 5.5% intergenic (‑48/+42) folP ← / ← ftsH 7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,324,9590CT 5.5% 1284.3 / 38.7 379intergenic (‑48/+42)folP/ftsH7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base C (169/189);  new base T (11/10);  total (180/199)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 394 total aligned reads displayed.

TGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGA  >  NC_000913/3324811‑3325108
                                                                                                                                                    |                                                                                                                                                      
tGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGccc                                                                                                                                                       >  1:4406295/1‑150 (MQ=255)
 gCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGcccc                                                                                                                                                      >  1:616018/1‑150 (MQ=255)
  cGTGCCACCATCCGAAAGGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑g                                                                                                                                                    >  1:1004810/1‑150 (MQ=255)
    tGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑ggg                                                                                                                                                  <  2:1816642/150‑1 (MQ=255)
     gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑gggg                                                                                                                                                 <  1:3599965/150‑1 (MQ=255)
       caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCAGTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCC‑GGGGtt                                                                                                                                               <  2:5958107/150‑1 (MQ=255)
        accaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGtt                                                                                                                                               <  1:911603/149‑1 (MQ=255)
          caTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCAGTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAACGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTc                                                                                                                                            <  2:6398142/150‑1 (MQ=255)
             ccGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGCGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGCGCACGCCCC‑GGGGGTTTCGGt                                                                                                                                         >  1:84912/1‑150 (MQ=255)
               gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTAc                                                                                                                                       >  1:4807548/1‑150 (MQ=255)
                aaaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACa                                                                                                                                      <  1:1400277/150‑1 (MQ=255)
                 aaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACaa                                                                                                                                     <  1:6093151/150‑1 (MQ=255)
                  aaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACaaa                                                                                                                                    >  2:1132278/1‑150 (MQ=255)
                     gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCTGAGCACGCCCC‑GGGGTTTTCGGTACAAATAc                                                                                                                                 <  1:1761115/150‑1 (MQ=255)
                     gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATAc                                                                                                                                 >  1:1101115/1‑150 (MQ=255)
                      aaTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGTTGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTAATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACa                                                                                                                                <  1:1435260/150‑1 (MQ=255)
                      aaTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACa                                                                                                                                <  2:2209640/150‑1 (MQ=255)
                      aaTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACa                                                                                                                                <  1:3098782/150‑1 (MQ=255)
                        tCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGt                                                                                                                              >  2:1619299/1‑150 (MQ=255)
                         cAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTc                                                                                                                             >  1:2305715/1‑150 (MQ=255)
                                   gTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCg                                                                                                                   <  1:6654830/150‑1 (MQ=255)
                                   gTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCg                                                                                                                   >  2:2417397/1‑150 (MQ=255)
                                    ttGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCgg                                                                                                                  <  2:4424013/150‑1 (MQ=255)
                                      gAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGa                                                                                                                >  2:3019400/1‑150 (MQ=255)
                                       aGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                               <  1:4141026/150‑1 (MQ=255)
                                          aTCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCGGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAcctg                                                                                                            >  1:71247/1‑149 (MQ=255)
                                            ccccATGACGTGAGGAGGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCAGGTATGAATTGATAAGCAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATCCAGTCATCTGATGCGGGAacttac                                                                                                          <  1:5434809/150‑1 (MQ=255)
                                              ccATTACGTGAGGATGGCTAAGGTCCTGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttactt                                                                                                        <  1:1675320/150‑1 (MQ=255)
                                              ccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGACAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttactt                                                                                                        <  2:3482115/150‑1 (MQ=255)
                                              ccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttactt                                                                                                        <  2:5272829/150‑1 (MQ=255)
                                              ccAGTACGTGAGGATGGCTACGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttactt                                                                                                        <  1:2478726/150‑1 (MQ=255)
                                                aTTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGt                                                                                                      >  2:2509521/1‑150 (MQ=255)
                                                 ttACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTc                                                                                                     <  1:3496758/150‑1 (MQ=255)
                                                 ttACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTc                                                                                                     <  2:2159608/150‑1 (MQ=255)
                                                 ttACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTc                                                                                                     <  2:5538066/150‑1 (MQ=255)
                                                    cGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGcc                                                                                                  <  1:5685999/150‑1 (MQ=255)
                                                    cGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGcc                                                                                                  >  1:2659371/1‑150 (MQ=255)
                                                       gAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTaa                                                                                               <  2:2602523/150‑1 (MQ=255)
                                                        aGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAAc                                                                                              >  1:1167387/1‑150 (MQ=255)
                                                          gATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTg                                                                                            <  1:390153/150‑1 (MQ=255)
                                                          gATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTg                                                                                            <  1:3067217/150‑1 (MQ=255)
                                                          gATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTg                                                                                            >  1:796322/1‑150 (MQ=255)
                                                            tGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGct                                                                                          <  2:627932/150‑1 (MQ=255)
                                                             ggCTAAGGTCCAGTGAAGTACCCTGGGCAAAGTGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctc                                                                                         <  2:3788427/150‑1 (MQ=255)
                                                              gCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                                        <  1:3816118/150‑1 (MQ=255)
                                                              gCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                                        >  2:4575334/1‑150 (MQ=255)
                                                              gCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                                        >  2:293644/1‑150 (MQ=255)
                                                               cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAAGTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTGTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTg                                                                                       >  2:4192342/1‑150 (MQ=255)
                                                                 aaGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGAc                                                                                     <  2:2901913/150‑1 (MQ=255)
                                                                   ggTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACAt                                                                                   <  1:3237227/150‑1 (MQ=255)
                                                                     tCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATgg                                                                                 >  2:5844029/1‑150 (MQ=255)
                                                                     tCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATgg                                                                                 <  1:4575334/150‑1 (MQ=255)
                                                                      ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                                <  1:1526871/150‑1 (MQ=255)
                                                                      ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                                <  2:4917311/150‑1 (MQ=255)
                                                                      ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                                <  2:6675995/150‑1 (MQ=255)
                                                                      ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                                <  1:1797461/150‑1 (MQ=255)
                                                                       cAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtg                                                                               >  2:4711337/1‑150 (MQ=255)
                                                                       cAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtg                                                                               <  1:2396102/150‑1 (MQ=255)
                                                                        aGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGGCATGgtgt                                                                              >  2:764081/1‑150 (MQ=255)
                                                                        aGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                              >  2:970250/1‑150 (MQ=255)
                                                                        aGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                              >  2:1292792/1‑150 (MQ=255)
                                                                              gTACCCTGGGCAAAGAGTATCATGATGTTATCCCTGGTATGAATAGAGAAGAAAAAACCCCGGAGCACGCCCC‑GCGGATTTCGGTACAAATACAGTCATCTGATTCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGCTGTTACACg                                                                        >  2:5998352/1‑150 (MQ=255)
                                                                               tACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCgg                                                                       <  1:189540/150‑1 (MQ=255)
                                                                                  ccTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGtt                                                                    >  1:6015456/1‑150 (MQ=255)
                                                                                   cTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTc                                                                   <  2:2305715/150‑1 (MQ=255)
                                                                                      ggCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGc                                                                >  1:4495634/1‑150 (MQ=255)
                                                                                             agTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgcg                                                         >  2:3153928/1‑150 (MQ=255)
                                                                                               tttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGt                                                       <  1:5330233/150‑1 (MQ=255)
                                                                                                  cATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTCCCCGGGTTCGGCGTACGCGGTtca                                                    >  1:6540846/1‑150 (MQ=255)
                                                                                                   atgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcat                                                   >  2:4049492/1‑150 (MQ=255)
                                                                                                    tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                                  <  2:4965079/150‑1 (MQ=255)
                                                                                                    tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                                  <  2:1229052/150‑1 (MQ=255)
                                                                                                    tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                                  <  1:3439762/150‑1 (MQ=255)
                                                                                                    tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                                  <  1:656827/150‑1 (MQ=255)
                                                                                                     gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                                  <  2:3076648/150‑1 (MQ=255)
                                                                                                        gTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTCCGCGGTTCATCAAcc                                              >  1:6394558/1‑150 (MQ=255)
                                                                                                          tATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAg                                                                                                                               <  1:3662676/67‑1 (MQ=255)
                                                                                                          tATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAg                                                                                                                               >  2:3662676/1‑67 (MQ=255)
                                                                                                                tGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagg                                      <  2:3351153/150‑1 (MQ=255)
                                                                                                                tGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagg                                      >  1:3145957/1‑150 (MQ=255)
                                                                                                                  gTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgaggag                                    >  1:1609201/1‑150 (MQ=255)
                                                                                                                  gTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgaggag                                    >  1:5000762/1‑150 (MQ=255)
                                                                                                                   tATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCTGGTTCGGCGTACGCGGTTCATAAACCGGACGAGGAGc                                   <  2:189840/150‑1 (MQ=255)
                                                                                                                        aTTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg                              >  2:3526685/1‑150 (MQ=255)
                                                                                                                         ttGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                             <  1:5495219/150‑1 (MQ=255)
                                                                                                                            aTAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACt                          >  1:2933480/1‑150 (MQ=255)
                                                                                                                              aaGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAc                        >  1:2612203/1‑150 (MQ=255)
                                                                                                                              aaGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAc                        >  2:26824/1‑150 (MQ=255)
                                                                                                                               acaaaaaaCCCCGGAGCACGCCCC‑GGGGATTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGGTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc                       >  1:4750507/3‑150 (MQ=255)
                                                                                                                                gAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCa                      >  2:5869547/1‑150 (MQ=255)
                                                                                                                                 aaaaaaCCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCAt                     >  1:5883242/1‑150 (MQ=255)
                                                                                                                                   aaaaCCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAc                        <  2:1492479/145‑1 (MQ=255)
                                                                                                                                     aaCCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTc                 <  1:5843586/150‑1 (MQ=255)
                                                                                                                                        cccGGGGCGTGCTCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtc                                                     <  2:5700622/111‑1 (MQ=255)
                                                                                                                                        cccGGGGCGTGCTCC‑GGGGTTTTAGGTACAAATAAAGTCATCTGATGCGGGAACTTACTTGTCGCCTAAATGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc              >  1:21926/1‑150 (MQ=255)
                                                                                                                                         ccGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTCCCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCa             >  1:2769735/1‑150 (MQ=255)
                                                                                                                                               cACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAAtt        >  2:1010269/1‑149 (MQ=255)
                                                                                                                                                 cGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGtt     <  2:4221702/150‑1 (MQ=255)
                                                                                                                                                   cccc‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAg   >  2:519465/1‑150 (MQ=255)
                                                                                                                                                    ccc‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGa  >  1:1339729/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCC‑GGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGA  >  NC_000913/3324811‑3325108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑