Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 557,866 G→T 5.0% intergenic (‑125/‑346) folD ← / → sfmA bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase/putative fimbrial protein SfmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913557,8660GT 5.0% 367.9 / 3.2 109intergenic (‑125/‑346)folD/sfmAbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase/putative fimbrial protein SfmA
Reads supporting (aligned to +/- strand):  ref base G (48/55);  new base T (2/4);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAATTTACCGCTTCA  >  NC_000913/557720‑558015
                                                                                                                                                  |                                                                                                                                                     
cGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCAc                                                                                                                                                    <  2:981959/150‑1 (MQ=255)
    aataatCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTc                                                                                                                                                >  2:1203769/1‑150 (MQ=255)
      taatCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGc                                                                                                                                              >  2:555054/1‑150 (MQ=255)
       aatCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGACGGTGCGCCTCACCTTCGca                                                                                                                                             <  2:591993/150‑2 (MQ=255)
        atCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGACGGTGCGCCTCACCTTCGcaa                                                                                                                                            <  2:1627459/150‑3 (MQ=255)
        atCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTg                                                                                                                                            <  1:1163776/150‑1 (MQ=255)
        atCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTg                                                                                                                                            <  1:1169593/150‑1 (MQ=255)
        atCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTg                                                                                                                                            <  1:130832/150‑1 (MQ=255)
          cTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTgag                                                                                                                                          <  1:468248/150‑1 (MQ=255)
            ttGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGACGGTGCGCCTCACCTTCGcaagagg                                                                                                                                        <  2:1218373/150‑7 (MQ=255)
               ctgcCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCTGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGa                                                                                                                                     <  1:1489497/150‑1 (MQ=255)
               ctgcCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGa                                                                                                                                     >  1:140391/1‑150 (MQ=255)
                tgcCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGaa                                                                                                                                    >  1:136679/1‑150 (MQ=255)
                 gcCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGaaa                                                                                                                                   >  1:783629/1‑150 (MQ=255)
                   cATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGc                                                                                                                                 >  1:954320/1‑150 (MQ=255)
                   cATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGc                                                                                                                                 <  2:1378212/150‑1 (MQ=255)
                    aTCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCt                                                                                                                                <  1:1347871/150‑1 (MQ=255)
                     tCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCtt                                                                                                                               >  1:96354/1‑150 (MQ=255)
                        gagagagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAt                                                                                                                            <  2:783629/150‑1 (MQ=255)
                        gagagagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAt                                                                                                                            <  1:1126063/150‑1 (MQ=255)
                         agagagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAtt                                                                                                                           >  2:1001228/1‑150 (MQ=255)
                         agagagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAtt                                                                                                                           >  2:168281/1‑150 (MQ=255)
                         agagagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAtt                                                                                                                           <  1:1049837/150‑1 (MQ=255)
                            gagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAttttt                                                                                                                        >  1:1420236/1‑150 (MQ=255)
                            gagGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCAttttt                                                                                                                        <  2:1075409/150‑1 (MQ=255)
                               gATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGt                                                                                                                     <  1:413567/150‑1 (MQ=255)
                                 ttCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTcc                                                                                                                   >  2:1052915/1‑150 (MQ=255)
                                  tCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCa                                                                                                                  <  1:1001228/150‑1 (MQ=255)
                                  tCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCa                                                                                                                  <  2:307469/150‑1 (MQ=255)
                                        tGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATgg                                                                                                            <  2:140852/150‑1 (MQ=255)
                                         gTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGc                                                                                                           >  2:1243082/1‑150 (MQ=255)
                                           tACGTAGATCGAAGGGTATGCGCCTATTTTGTCAGAAGCGTGTCGCGCTTTCAGGTTTCGTTTCATATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGt                                                                                                         <  2:177456/150‑1 (MQ=255)
                                           tACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGt                                                                                                         <  1:1191946/150‑1 (MQ=255)
                                             cGTCGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTg                                                                                                       <  2:1479457/150‑1 (MQ=255)
                                             cGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTg                                                                                                       >  1:911548/1‑150 (MQ=255)
                                             cGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTg                                                                                                       >  1:301821/1‑150 (MQ=255)
                                             cGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTg                                                                                                       >  1:704480/1‑150 (MQ=255)
                                             cGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTg                                                                                                       <  1:760641/150‑1 (MQ=255)
                                             cGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGGCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTg                                                                                                       <  1:1110342/150‑1 (MQ=255)
                                                 gATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTaaa                                                                                                   >  1:144204/1‑150 (MQ=255)
                                                 gATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTaaa                                                                                                   >  1:987582/1‑150 (MQ=255)
                                                  aacGAAGGGGATGCGCCTATTTAGTCAGAAGCGCGGCGCGCTGTCAGCGATCGTTTCAGATTTAACGCGTGCAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAAt                                                                                                  <  2:1563505/148‑1 (MQ=255)
                                                   tCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATc                                                                                                 >  1:739522/1‑150 (MQ=255)
                                                     gAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTg                                                                                               <  1:900376/150‑1 (MQ=255)
                                                          ggATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGCTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACt                                                                                          >  2:1253684/1‑150 (MQ=255)
                                                           gATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTGTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTGATTTTTCGTCCATGATGGCGTTGTAAATCTGGAAGTg                                                                                         <  2:856230/150‑1 (MQ=255)
                                                              gcgcCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGatt                                                                                      >  2:346619/1‑150 (MQ=255)
                                                                gcCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGattta                                                                                    >  1:578011/1‑150 (MQ=255)
                                                                   tATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGatttattt                                                                                 <  1:1598940/150‑1 (MQ=255)
                                                                    aTTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTc                                                                                >  1:1315809/1‑150 (MQ=255)
                                                                     ttttGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTcc                                                                               <  1:26704/150‑1 (MQ=255)
                                                                        tGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTg                                                                            <  2:386259/150‑1 (MQ=255)
                                                                         gTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGt                                                                           <  1:1189964/150‑1 (MQ=255)
                                                                          tCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTc                                                                          >  2:211852/1‑150 (MQ=255)
                                                                            aGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTa                                                                        <  2:80673/150‑1 (MQ=255)
                                                                              aaGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAg                                                                      >  1:1215240/1‑150 (MQ=255)
                                                                               aGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAgg                                                                     >  1:812245/1‑150 (MQ=255)
                                                                               aGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAgg                                                                     >  2:722115/1‑150 (MQ=255)
                                                                                  ggggCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGAtt                                                                  <  2:911548/150‑1 (MQ=255)
                                                                                   gggCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTa                                                                 >  1:374314/1‑150 (MQ=255)
                                                                                    ggCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTaa                                                                <  2:301821/150‑1 (MQ=255)
                                                                                     gcgcgcTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAg                                                               <  1:77668/150‑1 (MQ=255)
                                                                                           tGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAAtt                                                         >  1:453198/1‑150 (MQ=255)
                                                                                            gTCAGGTTTCTTTTCTGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGATGCTCAGCTTCGCTGAGACTGAAAGCTTGATTTTTCGTCCATGATGGCGTTGTTAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAAttt                                                        <  1:746060/150‑1 (MQ=255)
                                                                                               aGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAg                                                     <  1:72201/150‑1 (MQ=255)
                                                                                                cgtttccattcAGATTTATCGCGTGAAGAGAACTCTTGCGAAGGTGAGGCGAAACGTAGATGAGACTGAAAGCTTCATTTTTAGTACATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGa                                                    >  2:715089/2‑150 (MQ=255)
                                                                                                     cgtttcAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATAc                                               <  1:429172/150‑1 (MQ=255)
                                                                                                      gtttcAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATAcc                                              <  2:738267/150‑1 (MQ=255)
                                                                                                      gtttcAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATAcc                                              >  2:339824/1‑150 (MQ=255)
                                                                                                      gtttcAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATAcc                                              <  1:171028/150‑1 (MQ=255)
                                                                                                          cAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTCAGGATTAAGATAATTTAAGAAATACCTGAc                                          <  2:1020781/150‑1 (MQ=255)
                                                                                                          cAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGAc                                          <  2:525747/150‑1 (MQ=255)
                                                                                                          cAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGAc                                          <  1:26235/150‑1 (MQ=255)
                                                                                                          cAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGAc                                          >  2:1125185/1‑150 (MQ=255)
                                                                                                          cAGATTTATCGCGTGAAGAGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTAATTGTATAAGGATTAAGATAATTTAAGAAATACCTGAc                                          >  2:41913/1‑150 (MQ=255)
                                                                                                            gATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACaa                                        >  1:1582591/1‑150 (MQ=255)
                                                                                                             aTTTATCGCGGGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTTGACCTGATTTAATTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAat                                       <  2:12220/150‑1 (MQ=255)
                                                                                                              tttATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAata                                      <  1:1052915/150‑1 (MQ=255)
                                                                                                               ttATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAatat                                     <  2:957098/150‑1 (MQ=255)
                                                                                                                tATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTATAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAAtata                                    <  1:363632/150‑1 (MQ=255)
                                                                                                                    gcgTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATaaaaa                                <  1:874628/150‑1 (MQ=255)
                                                                                                                    gcgTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATaaaaa                                >  1:894250/1‑150 (MQ=255)
                                                                                                                        gAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAAt                            >  1:647796/1‑150 (MQ=255)
                                                                                                                          aGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAAttt                          >  1:561781/1‑150 (MQ=255)
                                                                                                                            cGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTc                        <  2:50269/150‑1 (MQ=255)
                                                                                                                             gACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCa                       >  1:244654/1‑150 (MQ=255)
                                                                                                                               gtctcagCGACGGTGCGCCTCACCTTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGc                     >  1:1311973/7‑150 (MQ=255)
                                                                                                                                tctcagCGACGGTGCGCCTCACCTTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGcc                    >  1:391067/6‑150 (MQ=255)
                                                                                                                                ctctTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGcc                    <  2:181830/150‑1 (MQ=255)
                                                                                                                                ctctTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGcc                    <  2:1119145/150‑1 (MQ=255)
                                                                                                                                  ctTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGatt                                                                                      >  1:1024193/1‑82 (MQ=255)
                                                                                                                                  ctTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGatt                                                                                      <  2:1024193/82‑1 (MQ=255)
                                                                                                                                    tGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGt                >  1:762842/1‑150 (MQ=255)
                                                                                                                                    tGCGAAGGTGAGGCGCACCGTCGCTGAGACAGAAAGCTTCATTTTTAGTACATGATGGCGTTGTAAATATGGAACTGATTTATTTGCTTGTCTAAGGATTAAGATAATTTAAGAAAAAACTGACAAAATAAAAAGAATTTTAACCCTGGt                >  2:985803/1‑150 (MQ=255)
                                                                                                                                     gCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTa               >  1:917278/1‑150 (MQ=255)
                                                                                                                                       gAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAAt             <  2:1393249/150‑1 (MQ=255)
                                                                                                                                          ggTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAATTTa          >  1:196809/1‑150 (MQ=255)
                                                                                                                                                gcgcACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAATTTACCGCtt    >  2:149629/1‑150 (MQ=255)
                                                                                                                                                 cgcACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAATTTACCGCTTc   >  2:406920/1‑150 (MQ=255)
                                                                                                                                                  gcACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAATTTACCGCTTCa  >  2:745308/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTCGCTGAGACTGAAAGCTTCATTTTTCGTCCATGATGGCGTTGTAAATCTGGAACTGATTTATTTCCTTGTCTAAGGATTAAGATAATTTAAGAAATACCTGACAATATAAAAAGAATTTTCAGCCTGGTAATTTACCGCTTCA  >  NC_000913/557720‑558015

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