Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,794,150 A→T 5.1% D55V (GAC→GTC)  rfaD → ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,1500AT 5.1% 426.3 / 8.1 136D55V (GAC→GTC) rfaDADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (67/61);  new base T (2/5);  total (70/66)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 137 total aligned reads displayed.

CGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCC  >  NC_000913/3794007‑3794299
                                                                                                                                               |                                                                                                                                                     
cgcgGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCt                                                                                                                                                 <  2:857039/150‑1 (MQ=255)
cgcgGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCt                                                                                                                                                 <  1:55115/150‑1 (MQ=255)
 gcgGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTcc                                                                                                                                                  <  1:584371/148‑1 (MQ=255)
 gcgGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTcc                                                                                                                                                  >  2:584371/1‑148 (MQ=255)
 gcgGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTg                                                                                                                                                <  2:481590/150‑1 (MQ=255)
 gcgGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTg                                                                                                                                                <  2:791391/150‑1 (MQ=255)
     gCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATcc                                                                                                                                            <  2:1035750/150‑1 (MQ=255)
     gCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATcc                                                                                                                                            >  1:940231/1‑150 (MQ=255)
     gCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATcc                                                                                                                                            <  2:491205/150‑1 (MQ=255)
      cTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATCAGGAAGTCTTCCTTATCCa                                                                                                                                           <  2:440395/150‑1 (MQ=255)
      cTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCCGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCa                                                                                                                                           <  1:873553/150‑1 (MQ=255)
        ttATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGAGGATCTGAATATCGCAGACTATATGGATCAGGAAGTCTTCCTTATCCata                                                                                                                                         <  2:3655/150‑3 (MQ=255)
         tATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGAt                                                                                                                                        >  2:1011708/1‑150 (MQ=255)
             ggCAGCAACATCGTTAAAGCCCTGAATTATAATTGCATGAACTATATTCTGGTGGTGGACAACCTGATAGACGGCACCAAGTTTGTTAACTTGTTGTATCTGAATATCTCAGACTATATGGATCAGGAAGTCTTCCTTATCCatatagtc                                                                                                                                    <  2:1369314/150‑5 (MQ=255)
              gcagcaACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtgg                                                                                                                                   <  2:1005098/150‑1 (MQ=255)
              gcagcaACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtgg                                                                                                                                   <  1:178107/150‑1 (MQ=255)
               cagcaACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtggc                                                                                                                                  <  1:150886/150‑1 (MQ=255)
                agcaACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtggcg                                                                                                                                 >  1:255234/1‑149 (MQ=255)
                  caACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGTTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtggctgg                                                                                                                               <  2:610448/150‑1 (MQ=255)
                    aCATCGTTAAAGCCCTGAATGATGAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCg                                                                                                                             >  2:1244371/1‑150 (MQ=255)
                    aCATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCg                                                                                                                             <  2:1151249/150‑1 (MQ=255)
                     cATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGa                                                                                                                            >  2:1197790/1‑150 (MQ=255)
                       tCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATGTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag                                                                                                                          <  1:1319325/150‑1 (MQ=255)
                       tCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag                                                                                                                          >  2:602802/1‑150 (MQ=255)
                       tCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag                                                                                                                          >  1:281672/1‑150 (MQ=255)
                           tAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGtt                                                                                                                      >  2:1404000/1‑150 (MQ=255)
                            aaaGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTc                                                                                                                     >  1:1256962/1‑150 (MQ=255)
                                cccTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCg                                                                                                                 <  2:376842/150‑1 (MQ=255)
                                   tGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATg                                                                                                              <  1:2182/150‑1 (MQ=255)
                                   tGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATg                                                                                                              <  1:431414/150‑1 (MQ=255)
                                    gAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGt                                                                                                             >  1:731405/1‑150 (MQ=255)
                                    gAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGt                                                                                                             >  1:1374552/1‑150 (MQ=255)
                                      aTGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATCAGGAAGTCTTCCTTATCCatatagt                                                                                                                                     <  1:106579/124‑4 (MQ=39)
                                      aTGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATCAGGAAGTCTTCCTTATCCatatagt                                                                                                                                     >  2:106579/1‑121 (MQ=39)
                                              ggCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAtt                                                                                                   >  2:1001749/1‑150 (MQ=255)
                                                cATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAtttt                                                                                                 >  1:1025031/1‑150 (MQ=255)
                                                 aTCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTc                                                                                                >  2:56704/1‑150 (MQ=255)
                                                  tCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTcc                                                                                               <  1:1191746/150‑1 (MQ=255)
                                                  tCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTcc                                                                                               <  2:414308/150‑1 (MQ=255)
                                                   cACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa                                                                                              >  2:574413/1‑150 (MQ=255)
                                                   cACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACATGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa                                                                                              >  2:643843/1‑150 (MQ=255)
                                                    aCCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCAc                                                                                             >  1:1126515/1‑150 (MQ=255)
                                                    aCCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCAc                                                                                             >  1:201120/1‑150 (MQ=255)
                                                      cGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGa                                                                                           <  2:1126515/150‑1 (MQ=255)
                                                      cGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGa                                                                                           <  2:201120/150‑1 (MQ=255)
                                                                gtggtggACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGc                                                                                 <  2:731405/150‑1 (MQ=255)
                                                                    tggACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCtt                                                                             >  1:1326526/1‑150 (MQ=255)
                                                                     ggACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTc                                                                            <  1:196977/150‑1 (MQ=255)
                                                                       aCAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTcca                                                                          >  2:467752/1‑150 (MQ=255)
                                                                        cAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTccac                                                                         <  1:1011708/150‑1 (MQ=255)
                                                                         aaCCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTccacc                                                                        <  1:127481/150‑1 (MQ=255)
                                                                          aCCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTccacca                                                                       >  1:786251/1‑150 (MQ=255)
                                                                          aCCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGAGATTTTCAACGAAGGCGCGTGCTCTTccacca                                                                       >  1:51677/1‑150 (MQ=255)
                                                                           ccTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCcaccac                                                                      >  1:12430/1‑150 (MQ=255)
                                                                           ccTGAAAGACGGAAACAATTTTGTGAACCTGGTGGATCTGAATATCGCAGACTAAATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTACGGAGATGTCGAAGCGATATTCCACGAAGGCGCGTGCTCTTCcaccac                                                                      >  2:926846/1‑150 (MQ=255)
                                                                            cTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCaccacc                                                                     <  1:452985/150‑1 (MQ=255)
                                                                            cTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCaccacc                                                                     >  1:871889/1‑150 (MQ=255)
                                                                              gAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGa                                                                   >  2:377201/1‑150 (MQ=255)
                                                                                aaGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGt                                                                 >  2:587683/1‑150 (MQ=255)
                                                                                 aGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTg                                                                <  1:1275772/150‑1 (MQ=255)
                                                                                    cGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGa                                                             >  2:349235/1‑150 (MQ=255)
                                                                                      gCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACg                                                           >  1:482014/1‑150 (MQ=255)
                                                                                        aCCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATTGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCTAAGAGTTCGGCGATGTCGAATCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGc                                                         <  2:118884/150‑1 (MQ=255)
                                                                                        aCCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGc                                                         <  2:281672/150‑1 (MQ=255)
                                                                                          cAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCaa                                                       >  1:1115219/1‑150 (MQ=255)
                                                                                           aaGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAg                                                      <  1:8615/150‑1 (MQ=255)
                                                                                            aGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGt                                                     <  2:789477/150‑1 (MQ=255)
                                                                                            aGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGt                                                     <  1:973775/150‑1 (MQ=255)
                                                                                             ttttGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGta                                                    <  1:587683/149‑1 (MQ=255)
                                                                                             gTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCAACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGtt                                                    >  2:1138722/1‑149 (MQ=255)
                                                                                              tttGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGtat                                                   <  2:186852/150‑1 (MQ=255)
                                                                                               ttGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGtata                                                  >  1:355291/1‑150 (MQ=255)
                                                                                                     aCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATgg                                            >  2:1074500/1‑150 (MQ=255)
                                                                                                      ccTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGa                                           <  2:871889/150‑1 (MQ=255)
                                                                                                       cTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGAt                                          <  1:1244371/150‑1 (MQ=255)
                                                                                                        tggtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATa                                         <  2:1256962/150‑1 (MQ=255)
                                                                                                        tggtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATa                                         >  1:1032134/1‑150 (MQ=255)
                                                                                                         ggtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATaa                                        >  2:1359868/1‑150 (MQ=255)
                                                                                                           tggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATaaca                                      >  2:1024564/1‑150 (MQ=255)
                                                                                                             gATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATaacaac                                    >  2:253470/1‑150 (MQ=255)
                                                                                                                cTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTAt                                 <  2:940231/150‑1 (MQ=255)
                                                                                                                 tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCTGCGATGTCGAAGCTAATTAAAACGAAGGCTCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTAta                                >  1:326384/1‑149 (MQ=255)
                                                                                                                 tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATc                                >  1:1192288/1‑150 (MQ=255)
                                                                                                                  gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCa                               >  2:1298715/1‑150 (MQ=255)
                                                                                                                  gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCa                               <  1:796139/150‑1 (MQ=255)
                                                                                                                     tatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATgg                                            <  2:231598/134‑1 (MQ=255)
                                                                                                                     tatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATgg                                            >  1:231598/1‑134 (MQ=255)
                                                                                                                         gCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTcc                        >  2:711606/1‑150 (MQ=255)
                                                                                                                          cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCa                       <  2:1048485/150‑1 (MQ=255)
                                                                                                                            gaatataTGGATAAGGAAGACTTCCTGATCCAGATTATGACTGGCGAAGAGTTAGGCGATGTCGAAGCGATTTTACACGAAGCCGCGTGCTCTTACACCACCGAGTGTGACGGCAAGTATATAATGGATAACAACTATCAATACTCCaat                     >  2:849841/4‑149 (MQ=255)
                                                                                                                                atatGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAgc                 <  2:173780/150‑1 (MQ=255)
                                                                                                                                   tGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAgctgc              <  2:1032134/150‑1 (MQ=255)
                                                                                                                                   tGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAgctgc              >  1:804695/1‑150 (MQ=255)
                                                                                                                                    ggATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTAGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAgctgct             >  1:1324269/1‑150 (MQ=255)
                                                                                                                                     gATAAGGAAGCCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGctgctg            >  1:1135037/1‑150 (MQ=255)
                                                                                                                                      aTAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCtgctgc           <  1:643843/150‑1 (MQ=255)
                                                                                                                                        aatgAAGACTTCCTGATCGAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCac         <  2:461127/147‑1 (MQ=255)
                                                                                                                                         aGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAgg                    >  1:585409/1‑138 (MQ=255)
                                                                                                                                         aGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCact        >  2:52646/1‑150 (MQ=255)
                                                                                                                                               aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGcc  <  2:316562/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCC  >  NC_000913/3794007‑3794299

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: