Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,154,601 A→T 5.5% S196C (AGC→TGC)  mdtA → multidrug efflux pump membrane fusion protein MdtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,154,6010AT 5.5% 351.1 / 5.9 108S196C (AGC→TGC) mdtAmultidrug efflux pump membrane fusion protein MdtA
Reads supporting (aligned to +/- strand):  ref base A (51/50);  new base T (2/4);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 108 total aligned reads displayed.

ACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGC  >  NC_000913/2154455‑2154744
                                                                                                                                                  |                                                                                                                                               
aCCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGcc                                                                                                                                              <  1:1064846/150‑1 (MQ=255)
 ccTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCg                                                                                                                                             >  2:880531/1‑150 (MQ=255)
   tGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGa                                                                                                                                           <  2:74519/150‑1 (MQ=255)
       gcgTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTAc                                                                                                                                       <  1:40255/150‑1 (MQ=255)
           tATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCa                                                                                                                                   >  1:587240/1‑150 (MQ=255)
            aTCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAc                                                                                                                                  >  2:787166/1‑150 (MQ=255)
             tCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAcc                                                                                                                                 >  2:542964/1‑150 (MQ=255)
             tCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAcc                                                                                                                                 <  1:1522610/150‑1 (MQ=255)
              caacaaCTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCa                                                                                                                                >  2:1141165/1‑150 (MQ=255)
               aacaacTGTCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCAAACAGGCGCTGGTCAGTGAAACCGAAGGCAACATTAAGGCTGATGAAGCAAGAGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAg                                                                                                                               >  2:1123741/1‑150 (MQ=255)
               aacaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAg                                                                                                                               >  1:1163535/1‑150 (MQ=255)
                acaacTCGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGt                                                                                                                              <  1:13408/150‑1 (MQ=255)
                 caacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTc                                                                                                                             <  2:587240/150‑1 (MQ=255)
                  aacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGTGCGGTAATCCGGCTCCAGTCg                                                                                                                            <  2:106745/150‑1 (MQ=255)
                    cTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGTGCGGTAATCCGGCTCCAGTCGag                                                                                                                          <  1:561683/150‑2 (MQ=255)
                     tGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGTGCGGTAATCCGGCTCCAGTCGagt                                                                                                                         <  1:953126/150‑3 (MQ=255)
                        cAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTc                                                                                                                      >  2:388564/1‑150 (MQ=255)
                        cAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTc                                                                                                                      >  2:1066969/1‑150 (MQ=255)
                           aaaCCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTcgcg                                                                                                                   >  2:1343520/1‑150 (MQ=255)
                             aCCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGtt                                                                                                                 >  1:745828/1‑150 (MQ=255)
                               cAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTgg                                                                                                               >  2:1299482/1‑150 (MQ=255)
                                 aTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtc                                                                                                             <  2:707673/150‑1 (MQ=255)
                                    tcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCa                                                                                                          <  2:686652/150‑1 (MQ=255)
                                      gTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGTGcgg                                                                                                                                             <  2:243498/113‑3 (MQ=255)
                                      gTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAg                                                                                                        <  1:354859/150‑1 (MQ=255)
                                          cccGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAgg                                                                                                    >  1:199813/1‑150 (MQ=255)
                                            cGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGtt                                                                                                  <  1:223571/150‑1 (MQ=255)
                                            cGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGtt                                                                                                  >  2:855216/1‑150 (MQ=255)
                                                 ggAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGt                                                                                             >  1:775327/1‑150 (MQ=255)
                                                 ggAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGt                                                                                             <  1:757906/150‑1 (MQ=255)
                                                  gAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGtt                                                                                            <  1:782286/150‑1 (MQ=255)
                                                    gCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTgg                                                                                          >  1:1219288/1‑150 (MQ=255)
                                                      tGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTgg                                                                                                               <  2:1440063/127‑1 (MQ=255)
                                                      tGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTgg                                                                                                               >  1:1440063/1‑127 (MQ=255)
                                                      tGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTa                                                                                        >  2:1521401/1‑150 (MQ=255)
                                                      tGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTa                                                                                        <  2:855465/150‑1 (MQ=255)
                                                      tGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTa                                                                                        <  2:68776/150‑1 (MQ=255)
                                                       ggATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTaa                                                                                       >  2:1413756/1‑150 (MQ=255)
                                                        gATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAAc                                                                                      >  1:591542/1‑150 (MQ=255)
                                                        gATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAAc                                                                                      >  2:864146/1‑150 (MQ=255)
                                                            cccAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCaaa                                                                                  >  1:631783/1‑150 (MQ=255)
                                                              cAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTACCCAAAtc                                                                                >  2:417515/1‑150 (MQ=255)
                                                                 cAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTcc                                                                             <  1:714358/150‑1 (MQ=255)
                                                                   ggCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAg                                                                           <  2:865289/150‑1 (MQ=255)
                                                                     cgcTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtg                                                                         >  2:584246/1‑150 (MQ=255)
                                                                      gcTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtgg                                                                        <  1:393976/150‑1 (MQ=255)
                                                                      gcTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtgg                                                                        <  2:686755/150‑1 (MQ=255)
                                                                        tGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGTTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtggtg                                                                      <  2:1312464/150‑1 (MQ=255)
                                                                         ggTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGa                                                                     >  1:1019155/1‑150 (MQ=255)
                                                                          gTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGAt                                                                    >  1:156244/1‑150 (MQ=255)
                                                                            cAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATac                                                                  >  2:104872/1‑150 (MQ=255)
                                                                              gagAAACAGAAGGCCCCATTAAGGCTGATGAAGCAAGCGTAGCCAGCGCCCAGCTGCAACTCGACAGCAGCCAGATAACCGCACCAGTCGATCGTCGCGTTGGTCTCAAGCACGTTGATGTTCGTAACAAACTCTCCACTGCTGATacaa                                                                >  2:344469/3‑148 (MQ=255)
                                                                              gTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATacca                                                                >  2:1408122/1‑150 (MQ=255)
                                                                                gAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATaccacc                                                              >  2:649242/1‑150 (MQ=255)
                                                                                gAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATaccacc                                                              >  2:636701/1‑150 (MQ=255)
                                                                                  aaCCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCgg                                                            >  1:1014104/1‑150 (MQ=255)
                                                                                  aaCCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCgg                                                            <  2:666015/150‑1 (MQ=255)
                                                                                   atctaaGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCTCGCAGATGCAACTGGACTGGAGCCGTATTACCGCACCATTCGATGGTCGCGTTGGTCTCAAGCATGTTGATTTTGGTAACCAAATCTCCAGTGGTGATACCACCggg                                                           <  2:855340/146‑1 (MQ=255)
                                                                                    ccGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCACACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGa                                                          >  2:805756/1‑150 (MQ=255)
                                                                                          gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                    >  1:400645/1‑150 (MQ=255)
                                                                                          gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                    >  2:739254/1‑150 (MQ=255)
                                                                                          gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                    <  2:156244/150‑1 (MQ=255)
                                                                                           cACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGAt                                                                    <  2:1339252/133‑1 (MQ=255)
                                                                                           cACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGAt                                                                    >  1:1339252/1‑133 (MQ=255)
                                                                                           cACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtggt                                                   >  2:1255205/1‑150 (MQ=255)
                                                                                           cACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtggt                                                   <  2:199813/150‑1 (MQ=255)
                                                                                             ccATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGa                                                 <  2:64218/150‑1 (MQ=255)
                                                                                              cATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGAt                                                >  1:1175539/1‑150 (MQ=255)
                                                                                                 tAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCAc                                             <  1:1343520/150‑1 (MQ=255)
                                                                                                    ggCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtc                                                                                                             >  1:1016749/1‑83 (MQ=255)
                                                                                                    ggCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtc                                                                                                             <  2:1016749/83‑1 (MQ=255)
                                                                                                        gatgaAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACg                                      <  2:1119803/150‑1 (MQ=255)
                                                                                                        gatgaAGCAAGCGTTGCCAGCCCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACg                                      >  1:351038/1‑150 (MQ=255)
                                                                                                          tgaAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCa                                    <  2:591542/150‑1 (MQ=255)
                                                                                                             agcaagcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCCCCCAGACGCATcc                                 >  2:1333226/1‑150 (MQ=255)
                                                                                                             agcaagcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATcc                                 <  1:1066969/150‑1 (MQ=255)
                                                                                                               caagcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTa                               <  1:864146/150‑1 (MQ=255)
                                                                                                                aagcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTAt                              <  1:1413756/150‑1 (MQ=255)
                                                                                                                 agcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATc                             <  1:1255205/150‑1 (MQ=255)
                                                                                                                  gcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCg                            <  1:1170289/150‑1 (MQ=255)
                                                                                                                    gTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGAt                          >  1:576673/1‑150 (MQ=255)
                                                                                                                     ttGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGGAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGAtt                         >  1:1432629/1‑150 (MQ=255)
                                                                                                                       gCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTa                       >  1:585581/1‑150 (MQ=255)
                                                                                                                       gCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTa                       >  1:1308393/1‑150 (MQ=255)
                                                                                                                         cAGCGCGCCGCTGCAACTCGACTGGCGCCGGCTTACCGCACCCGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTCACCAAATCTCCAGTGGTGATACCACCGGGATCGGGGTGCGCCCCCGGCCGCCTCCTATGGATTTAGt                     >  2:474368/1‑150 (MQ=255)
                                                                                                                         cAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGt                     <  2:833725/150‑1 (MQ=255)
                                                                                                                           gcgcgcAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCt                   <  1:1328073/150‑1 (MQ=255)
                                                                                                                            cgcgcAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCtt                  <  1:636701/150‑1 (MQ=255)
                                                                                                                              cgcAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTa                <  1:1159627/150‑1 (MQ=255)
                                                                                                                               gcAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTAc               >  2:1186667/1‑150 (MQ=255)
                                                                                                                                cAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTAcc              <  2:1148168/150‑1 (MQ=255)
                                                                                                                                cAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTAcc              <  1:1141165/150‑1 (MQ=255)
                                                                                                                                cAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGCGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTAcc              <  2:1115570/150‑1 (MQ=255)
                                                                                                                                 aGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTAccc             <  2:351038/150‑1 (MQ=255)
                                                                                                                                  gCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTTGTCTCAAGCAGGTTGATGTTTGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCt            <  2:1163535/150‑1 (MQ=255)
                                                                                                                                       ccatCGACTGGTGCGGTAATCCGGCTCCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCgg       >  1:1410273/4‑150 (MQ=255)
                                                                                                                                       aaCTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCgg       <  2:1219288/150‑1 (MQ=255)
                                                                                                                                         cTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGaa     <  1:444259/150‑1 (MQ=255)
                                                                                                                                           cGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtggtg                                                                      >  2:1313940/1‑83 (MQ=255)
                                                                                                                                            gACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGc  <  1:880531/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGC  >  NC_000913/2154455‑2154744

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: ‑