Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,161,057 C→T A891V (GCG→GTG)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,161,0570CT96.1% 352.9 / 2.8 103A891V (GCG→GTG) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (49/50);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 103 total aligned reads displayed.

CGCAGGTGTTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGCGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGACGCCGC  >  NC_000913/2160910‑2161203
                                                                                                                                                   |                                                                                                                                                  
cgcAGGTGTTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCt                                                                                                                                                 <  2:497709/151‑1 (MQ=255)
  cAGGTGTTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGt                                                                                                                                               <  1:323220/151‑1 (MQ=255)
     gtgtTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTgg                                                                                                                                            >  2:1663024/1‑151 (MQ=255)
       gtTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGgcg                                                                                                                                          <  2:1328883/151‑1 (MQ=255)
       gtTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGgc                                                                                                                                           >  2:1004802/1‑150 (MQ=255)
          ccAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTTCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTg                                                                                                                                       >  1:797684/1‑151 (MQ=255)
          ccAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCt                                                                                                                                        <  2:1562983/150‑1 (MQ=255)
            aGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTgg                                                                                                                                      >  1:1349078/1‑150 (MQ=255)
            aGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGa                                                                                                                                     <  1:1229122/151‑1 (MQ=255)
            aGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAATGCTTTTgggcct                                                                                                                                        >  1:847669/1‑144 (MQ=255)
                gaCGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTg                                                                                                                                 >  2:190966/1‑151 (MQ=255)
                  cGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGtt                                                                                                                               <  2:357080/151‑1 (MQ=255)
                  cGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGtt                                                                                                                               <  2:955306/151‑1 (MQ=255)
                      gAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAAt                                                                                                                           <  2:63051/151‑1 (MQ=255)
                           cGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGcccc                                                                                                                      >  2:167398/1‑151 (MQ=255)
                            gCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCg                                                                                                                     <  2:557430/151‑1 (MQ=255)
                              aGGTGATCCTGATTCGCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGtt                                                                                                                   <  1:1654180/151‑1 (MQ=255)
                              aGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGtt                                                                                                                   <  1:524576/151‑1 (MQ=255)
                              aGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGtt                                                                                                                   <  1:1000702/151‑1 (MQ=255)
                                    tCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCt                                                                                                             >  2:314899/1‑151 (MQ=255)
                                     ccTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCt                                                                                                             <  1:1336471/150‑1 (MQ=255)
                                     ccTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTa                                                                                                            <  2:255965/151‑1 (MQ=255)
                                        gATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTaa                                                                                                           <  1:470329/149‑1 (MQ=255)
                                          ttATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCg                                                                                                        >  2:1001025/1‑150 (MQ=255)
                                           tATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGcc                                                                                                      >  2:1633970/1‑151 (MQ=255)
                                            aTCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGcc                                                                                                      <  2:1392083/150‑1 (MQ=255)
                                             tCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCt                                                                                                    >  1:865835/1‑151 (MQ=255)
                                              cgccgcCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTg                                                                                                   <  2:535594/151‑1 (MQ=255)
                                                ccgccATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGAt                                                                                                 <  2:232969/151‑1 (MQ=255)
                                                ccgccATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGAt                                                                                                 <  2:192430/151‑1 (MQ=255)
                                                   tcATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAgg                                                                                              >  1:316654/2‑151 (MQ=255)
                                                     aTCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGa                                                                                            <  2:949199/151‑1 (MQ=255)
                                                      tCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGAt                                                                                           >  2:683138/1‑151 (MQ=255)
                                                        gCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATc                                                                                          >  2:440409/1‑150 (MQ=255)
                                                         ccACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCa                                                                                         <  2:160019/150‑1 (MQ=255)
                                                         ccACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCAt                                                                                        <  1:56220/151‑1 (MQ=255)
                                                         ccACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCAt                                                                                        >  1:509009/1‑151 (MQ=255)
                                                         ccACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCAt                                                                                        <  1:334541/151‑1 (MQ=255)
                                                         ccACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCAt                                                                                        >  1:1267321/1‑151 (MQ=255)
                                                          cACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATg                                                                                       >  1:1339667/1‑151 (MQ=255)
                                                          cACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATg                                                                                       >  1:1677882/1‑151 (MQ=255)
                                                           aCGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGc                                                                                      >  2:372501/1‑151 (MQ=255)
                                                            cGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCt                                                                                     <  2:367370/151‑1 (MQ=255)
                                                             ggTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCt                                                                                     >  2:1279024/1‑150 (MQ=255)
                                                               tgtATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTAtt                                                                                  >  1:630005/1‑151 (MQ=255)
                                                                  atatCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTaa                                                                                <  1:436334/150‑1 (MQ=255)
                                                                  atatCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTaa                                                                                <  2:813116/150‑1 (MQ=255)
                                                                  atatCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAAt                                                                               >  2:72009/1‑151 (MQ=255)
                                                                   tatCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATc                                                                              >  1:1036/1‑151 (MQ=255)
                                                                      cGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCgg                                                                            <  2:1608243/150‑1 (MQ=255)
                                                                      cGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGc                                                                           <  1:330405/151‑1 (MQ=255)
                                                                      cGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGCGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGc                                                                           >  2:1065374/1‑151 (MQ=255)
                                                                         gCTGGGTATCCTTTACGAGAGTTACGTACATCCGCGGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCAt                                                                         <  1:1624042/150‑1 (MQ=255)
                                                                            gggTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGt                                                                      <  2:957508/150‑1 (MQ=255)
                                                                              gTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTga                                                                    >  1:450592/1‑150 (MQ=255)
                                                                                aTCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTgaaga                                                                 <  1:72006/151‑1 (MQ=255)
                                                                                  ccTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTgaagaa                                                                <  1:1172440/150‑1 (MQ=255)
                                                                                   cTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTgaagaa                                                                <  2:243325/149‑1 (MQ=255)
                                                                                    tttACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGaaaa                                                              >  2:525599/1‑150 (MQ=255)
                                                                                     ttACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGaaaaa                                                             <  2:1698111/150‑1 (MQ=255)
                                                                                       aCGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACg                                                           <  2:1267417/150‑1 (MQ=255)
                                                                                        cGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGCGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGcc                                                         >  1:116302/1‑151 (MQ=255)
                                                                                         gagaGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGcc                                                         >  2:703355/1‑150 (MQ=255)
                                                                                         gagaGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGcc                                                         >  1:537862/1‑150 (MQ=255)
                                                                                          agagTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCAt                                                       >  1:75921/1‑151 (MQ=255)
                                                                                              ttACGTACATCCGCTGACGATTCTCTCCGCCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTat                                                    <  2:628268/150‑1 (MQ=255)
                                                                                                aCGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTatga                                                  <  2:420046/150‑1 (MQ=255)
                                                                                                aCGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTatga                                                  <  2:1067107/150‑1 (MQ=255)
                                                                                                  gTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATgg                                               <  1:745852/151‑1 (MQ=255)
                                                                                                     cATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGCGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCg                                            <  2:1501137/151‑1 (MQ=255)
                                                                                                      aTCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCg                                            >  1:1226401/1‑150 (MQ=255)
                                                                                                       tCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGAt                                          >  2:236517/1‑151 (MQ=255)
                                                                                                       tCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGAt                                          >  1:1401708/1‑151 (MQ=255)
                                                                                                        ccGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGAtt                                         <  2:1619654/151‑1 (MQ=255)
                                                                                                         cGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGGGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTGTTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGGCGAttt                                        >  1:243994/1‑151 (MQ=255)
                                                                                                         cGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGAttt                                        >  1:1593035/1‑151 (MQ=255)
                                                                                                           cTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTg                                      >  2:944703/1‑151 (MQ=255)
                                                                                                            tGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGAtttt                                       >  1:48346/1‑149 (MQ=255)
                                                                                                            tGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGAtttt                                       <  2:48352/149‑1 (MQ=255)
                                                                                                             gACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTg                                      <  1:1329232/149‑1 (MQ=255)
                                                                                                                 aTTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGGGCTGTTCAATGCCCCGTTCAGCCTAAACGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGGCGGTTTTTCGCTTg                                >  1:1716804/1‑151 (MQ=255)
                                                                                                                 aTTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCACCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTg                                >  1:61548/1‑151 (MQ=255)
                                                                                                                      ctcCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGaa                              <  2:1417806/148‑1 (MQ=255)
                                                                                                                      ctcCACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGcc                           >  2:1698906/1‑151 (MQ=255)
                                                                                                                        ccACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGccc                          >  1:958743/1‑150 (MQ=255)
                                                                                                                         cACCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCa                         >  1:903650/1‑150 (MQ=255)
                                                                                                                          aCCCTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAAc                       >  2:1156935/1‑151 (MQ=255)
                                                                                                                           cccTGCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACg                      >  1:1422090/1‑151 (MQ=255)
                                                                                                                               gCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCa                   <  2:1039975/150‑1 (MQ=255)
                                                                                                                               gCCCTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCa                   <  1:1369847/150‑1 (MQ=255)
                                                                                                                                 ccTCGGCGGGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACgg                <  2:1248859/151‑1 (MQ=255)
                                                                                                                                    cggcggGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGGTAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTa              >  2:90481/1‑150 (MQ=255)
                                                                                                                                    cggcggGCGTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTa              <  1:190952/150‑1 (MQ=255)
                                                                                                                                    cggcggGCGTTGGAGCGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGt               <  1:1658045/149‑1 (MQ=255)
                                                                                                                                         ggcgTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGa        <  1:310642/151‑1 (MQ=255)
                                                                                                                                          gcgTTGGAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGAc       >  2:1014650/1‑151 (MQ=255)
                                                                                                                                          gcgTTGGAGTGCTGTTGGCGCGGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGAc       <  2:527670/151‑1 (MQ=255)
                                                                                                                                               ggAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGAcgccgc  >  2:40452/1‑151 (MQ=255)
                                                                                                                                                gAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGAcgccgc  >  1:895056/1‑150 (MQ=255)
                                                                                                                                                gAGTGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGAcgccgc  >  2:96753/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAGGTGTTCCAGGAGACGATGAACTCGCAGGTGATCCTGATTATCGCCGCCATCGCCACGGTGTATATCGTGCTGGGTATCCTTTACGAGAGTTACGTACATCCGCTGACGATTCTCTCCACCCTGCCCTCGGCGGGCGTTGGAGCGCTGTTGGCGCTGGAGCTGTTCAATGCCCCGTTCAGCCTAATCGCCCTGATAGGGATCATGCTATTAATCGGCATCGTGAAGAAAAACGCCATTATGATGGTCGATTTTGCGCTTGAAGCCCAACGGCACGGTAACCTGACGCCGC  >  NC_000913/2160910‑2161203

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