Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,659,470 G→T 7.3% G80V (GGC→GTC) ‡ mdtE → multidrug efflux pump membrane fusion protein MdtE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,659,4700GT 7.3% 223.8 / 6.0 68G80V (GGC→GTC) ‡mdtEmultidrug efflux pump membrane fusion protein MdtE
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base T (3/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

GCCATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGCTGA  >  NC_000913/3659322‑3659619
                                                                                                                                                    |                                                                                                                                                     
gccATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGc                                                                                                                                                      <  1:1398369/150‑1 (MQ=255)
       cGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAg                                                                                                                                               >  2:128253/1‑150 (MQ=255)
       cGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAg                                                                                                                                               <  2:1387532/150‑1 (MQ=255)
        gCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGt                                                                                                                                              <  1:1183294/150‑1 (MQ=255)
           tGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGaa                                                                                                                                           >  1:416235/1‑150 (MQ=255)
            gAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAGttgcg                                                                                                                                          <  1:937546/150‑5 (MQ=255)
                tCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAgg                                                                                                                                      >  2:444853/1‑150 (MQ=255)
                      tcgtcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAg                                                                                                                                               <  2:595425/135‑1 (MQ=255)
                      tcgtcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAg                                                                                                                                               >  1:595425/1‑135 (MQ=255)
                      tcgtcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGAtt                                                                                                                                <  2:1142013/150‑1 (MQ=255)
                         tcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGc                                                                                                                             <  2:982442/150‑1 (MQ=255)
                           acacTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTg                                                                                                                           <  2:1144389/150‑1 (MQ=255)
                            cacTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGt                                                                                                                          <  1:372273/150‑1 (MQ=255)
                               tctcCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATc                                                                                                                       <  1:937146/150‑1 (MQ=255)
                                  cccccGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGa                                                                                                                    >  2:1126206/1‑150 (MQ=255)
                                         ttCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATcc                                                                                                             >  2:607226/1‑150 (MQ=255)
                                         ttCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATcc                                                                                                             <  1:1110157/150‑1 (MQ=255)
                                           cGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGATTTACGTCCCCAGGTGGGCGGTATTATCATTAAACGAAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTg                                                                                                           <  2:400372/150‑1 (MQ=255)
                                           cGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTg                                                                                                           >  2:819886/1‑150 (MQ=255)
                                            ggTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGc                                                                                                          >  1:1187916/1‑150 (MQ=255)
                                               cAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCCAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCATGCAcc                                                                                                       >  1:1356489/1‑150 (MQ=255)
                                                     gtTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACa                                                                                                 >  1:211213/1‑150 (MQ=255)
                                                       tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                               >  2:1552947/1‑150 (MQ=255)
                                                       tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                               >  2:457137/1‑150 (MQ=255)
                                                       tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                               >  1:1072585/1‑150 (MQ=255)
                                                        gAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGc                                                                                              >  1:548134/1‑150 (MQ=255)
                                                        gAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGc                                                                                              >  2:754468/1‑150 (MQ=255)
                                                         aGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGAGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGcc                                                                                             >  2:1348343/1‑150 (MQ=255)
                                                          gCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCg                                                                                            <  1:128253/150‑1 (MQ=255)
                                                             aaTTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGc                                                                                         >  1:1375170/1‑150 (MQ=255)
                                                              aTTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCt                                                                                        <  2:416235/150‑1 (MQ=255)
                                                                tGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTaa                                                                                      <  1:244670/150‑1 (MQ=255)
                                                                  cccGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAAc                                                                                    <  1:1105010/150‑1 (MQ=255)
                                                                  cccGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAAc                                                                                    >  2:314534/1‑150 (MQ=255)
                                                                    cGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTc                                                                                  >  1:442196/1‑150 (MQ=255)
                                                                        aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                              <  1:444853/150‑1 (MQ=255)
                                                                        aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                              <  1:559898/150‑1 (MQ=255)
                                                                        aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                              <  1:882211/150‑1 (MQ=255)
                                                                        aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                              <  1:1126206/150‑1 (MQ=255)
                                                                                  cTTATGAAGTTGCCGAGATACGTCCCGAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGc                                                                    >  2:1387978/1‑150 (MQ=255)
                                                                                  cTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGc                                                                    >  2:1393701/1‑150 (MQ=255)
                                                                                   ttATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCt                                                                   >  2:438238/1‑150 (MQ=255)
                                                                                   ttATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCt                                                                   <  1:1652229/150‑1 (MQ=255)
                                                                                    tATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTg                                                                  >  2:221438/1‑150 (MQ=255)
                                                                                      tGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGc                                                                >  1:1511611/1‑150 (MQ=255)
                                                                                       gAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCg                                                               <  2:1375170/150‑1 (MQ=255)
                                                                                            tGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCggtaga                                                          >  2:889621/1‑145 (MQ=255)
                                                                                            tGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgc                                                          >  2:315881/1‑150 (MQ=255)
                                                                                             gCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgcg                                                         >  1:1035729/1‑150 (MQ=255)
                                                                                               cGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGct                                                       <  1:967524/150‑1 (MQ=255)
                                                                                                    tACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTa                                                  >  2:1300303/1‑150 (MQ=255)
                                                                                                            cAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAACGCGCTCTCTACCGCCAGc                                          >  2:1548130/1‑150 (MQ=255)
                                                                                                             aGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCa                                         <  2:1072585/150‑1 (MQ=255)
                                                                                                              ggTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCaa                                        <  1:1387978/150‑1 (MQ=255)
                                                                                                              ggTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCaa                                        <  1:1393701/150‑1 (MQ=255)
                                                                                                                      gTATTATCATTAAACGCTAGTTTATCGAAGGCGATATAGTGAACGAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGACCTAGCTAAACTCCGCCAAAGGCTCGCTGGCTTAAGCGCTCTCTACCGCCAGCAATGCCCGCa                                <  2:204124/150‑1 (MQ=255)
                                                                                                                      gTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCa                                <  2:109745/150‑1 (MQ=255)
                                                                                                                          tatCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCAc                            >  1:549639/1‑150 (MQ=255)
                                                                                                                              aTTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCttt                        >  1:248523/1‑150 (MQ=255)
                                                                                                                                tAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTaa                      <  2:1187916/150‑1 (MQ=255)
                                                                                                                                  aaCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAcc                    >  2:1275503/1‑150 (MQ=255)
                                                                                                                                      tcaCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCa                >  1:31790/3‑150 (MQ=255)
                                                                                                                                      cAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCa                <  1:315881/150‑1 (MQ=255)
                                                                                                                                       caCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAg               >  2:255802/2‑150 (MQ=255)
                                                                                                                                                gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGc                  >  2:345726/1‑138 (MQ=255)
                                                                                                                                                gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGc                  <  1:345726/138‑1 (MQ=255)
                                                                                                                                                gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTg      >  1:271175/1‑150 (MQ=255)
                                                                                                                                                gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATAAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTg      >  2:324291/1‑150 (MQ=255)
                                                                                                                                                   ggCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGCTg   <  2:1511611/150‑1 (MQ=255)
                                                                                                                                                    gCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAACCTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGCTGa  >  1:659988/1‑150 (MQ=255)
                                                                                                                                                    gCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGCTGa  <  1:819886/150‑1 (MQ=255)
                                                                                                                                                    gCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGCTGa  <  1:152174/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGCTGA  >  NC_000913/3659322‑3659619

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: