Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,408,537 A→T 5.7% E163V (GAG→GTG)  zntB → Zn(2(+)):H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,408,5370AT 5.7% 272.0 / 5.4 88E163V (GAG→GTG) zntBZn(2(+)):H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base T (3/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGCCA  >  NC_000913/1408390‑1408686
                                                                                                                                                   |                                                                                                                                                     
gggtAATTGGTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagc                                                                                                                                                     <  1:683188/147‑1 (MQ=255)
ggTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagc                                                                                                                                                     <  2:1155575/150‑1 (MQ=255)
ggTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGCGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagc                                                                                                                                                     <  1:1097013/150‑1 (MQ=255)
 gTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagca                                                                                                                                                    <  1:1166465/150‑1 (MQ=255)
 gTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagca                                                                                                                                                    >  1:1280946/1‑150 (MQ=255)
 gTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagca                                                                                                                                                    >  1:1275693/1‑150 (MQ=255)
     aTTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagcagc                                                                                                                                                  <  1:1020248/148‑1 (MQ=255)
     aTTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagcagc                                                                                                                                                  >  2:1020248/1‑148 (MQ=255)
     aTTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTg                                                                                                                                                >  1:869151/1‑150 (MQ=255)
      ttGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTACAGTGAATTTATCGAGCAGCTGc                                                                                                                                               >  2:1089586/1‑150 (MQ=255)
        gTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCAc                                                                                                                                             <  1:1043362/150‑1 (MQ=255)
         tttCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGTCTGGTGGTTGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGGTGCACg                                                                                                                                            <  2:269042/150‑1 (MQ=255)
            cGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATa                                                                                                                                         >  1:1158633/1‑150 (MQ=255)
             gACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATaa                                                                                                                                        >  1:27829/1‑150 (MQ=255)
              aCCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGG‑AGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGTGCAGCTGCTCGATaaat                                                                                                                                      <  2:685232/150‑2 (MQ=255)
               cccGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATaaaa                                                                                                                                      <  1:931838/150‑1 (MQ=255)
               cccGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATaaaa                                                                                                                                      <  1:587796/150‑1 (MQ=255)
               cccGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATaaaa                                                                                                                                      >  2:292198/1‑150 (MQ=255)
                 cGACAACGCAACGGGCGGGCGCTGGACGCTGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGTTTGCGGGTGAGGTGTGGTGGATGTGTGCGATGCGTTGAGGGATCATTCCCGTGATTTTATCGAGCAGCTGCACGATAAAAtt                                                                                                                                    <  1:125241/150‑1 (MQ=255)
                 cGACAACGCAAAGTGATGGCGCTTGAAGATGTGGTGAGCGATCTGGAAGAGGCAACGGGTACGACCGATTTCGGGGGATGTCTGGTGGATGTGTTCGATGCGTTGACCTATCATTACATTGAATTTATCGAGCAGCTTCACGATAAAAtt                                                                                                                                    <  1:130698/150‑1 (MQ=255)
                    cAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATc                                                                                                                                 <  2:869151/150‑1 (MQ=255)
                     aaCGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGGCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGGTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCg                                                                                                                                <  2:1339242/150‑1 (MQ=255)
                      aCGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGa                                                                                                                               >  2:646301/1‑150 (MQ=255)
                       cGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGAc                                                                                                                              <  2:375072/150‑1 (MQ=255)
                       cGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGAc                                                                                                                              >  1:286493/1‑150 (MQ=255)
                         cAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCg                                                                                                                                <  2:1339181/146‑1 (MQ=255)
                         cAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCg                                                                                                                                >  1:1339181/1‑146 (MQ=255)
                         cAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCt                                                                                                                            >  1:1492028/1‑150 (MQ=255)
                         cAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCt                                                                                                                            >  2:500218/1‑150 (MQ=255)
                          aaaGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCtt                                                                                                                           >  2:996506/1‑150 (MQ=255)
                             gTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGTGCAGCTGCTCGATaaattcac                                                                                                                                  >  1:802350/1‑135 (MQ=255)
                             gTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGTGCAGCTGCTCGATaaattcac                                                                                                                                  <  2:802350/140‑6 (MQ=255)
                                 tGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATa                                                                                                                    <  2:1528237/150‑1 (MQ=255)
                                    cgcTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAAtc                                                                                                                 <  1:911103/150‑1 (MQ=255)
                                          aCGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTg                                                                                                           <  2:126567/150‑1 (MQ=255)
                                            gATGTGGTGTGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTTGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTCCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGAt                                                                                                         <  2:432775/150‑1 (MQ=255)
                                             aTGTGGTGAGCGATCTGGAAGAGGGGACGGGTCCGACGGATTTCGGGGGATGGGTGGTGGATGTGTGCGATGCGTTGACCTATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATc                                                                                                        <  2:803598/150‑1 (MQ=255)
                                             aTGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATc                                                                                                        >  1:1122468/1‑150 (MQ=255)
                                                 ggtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGGAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCa                                                                                                    <  2:1112491/150‑1 (MQ=255)
                                                 ggtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCa                                                                                                    >  2:1223633/1‑150 (MQ=255)
                                                  gtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCaa                                                                                                   <  1:616938/150‑1 (MQ=255)
                                                   tgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCaaa                                                                                                  <  1:605981/150‑1 (MQ=255)
                                                    gAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAAt                                                                                                 <  2:903227/150‑1 (MQ=255)
                                                     aGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTTGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAAtt                                                                                                <  1:95340/150‑1 (MQ=255)
                                                        gATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCa                                                                                             >  1:1350852/1‑150 (MQ=255)
                                                           gtGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCAcct                                                                                          <  2:586262/149‑2 (MQ=255)
                                                            tGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACcgc                                                                                         >  1:877312/1‑150 (MQ=255)
                                                             ggAAGAGGGCACCGGTCCGACCCATTCCGGGCGATGCCTGGTGGATGTGTGCGAAGCGATGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTAAAGATAATCTACTTGATCAGCAAATTCCACcgcg                                                                                        >  2:652939/1‑150 (MQ=255)
                                                                agagGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTgg                                                                                     >  1:1505303/1‑150 (MQ=255)
                                                                     gCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGATTTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTcc                                                                                <  2:1275693/150‑1 (MQ=255)
                                                                     gCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTcc                                                                                <  2:1280946/150‑1 (MQ=255)
                                                                      cACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCt                                                                               >  1:347100/1‑150 (MQ=255)
                                                                      cACGGGTCAGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCt                                                                               >  2:1318714/1‑150 (MQ=255)
                                                                        cGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTgg                                                                             <  2:1492028/150‑1 (MQ=255)
                                                                         gggTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGc                                                                            <  1:1407484/150‑1 (MQ=255)
                                                                         gggTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGc                                                                            <  2:1158633/150‑1 (MQ=255)
                                                                          ggTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGAACAGCAAATTCCACCGCGTGGATTCCTGGct                                                                           >  1:1345048/1‑150 (MQ=255)
                                                                                accgaTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTctgctg                                                                     >  2:1401994/1‑150 (MQ=255)
                                                                                accgaTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTctgctg                                                                     <  1:64012/150‑1 (MQ=255)
                                                                                  cgaTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCt                                                                               >  2:729640/1‑138 (MQ=255)
                                                                                  cgaTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCt                                                                               <  1:729640/138‑1 (MQ=255)
                                                                                   gaTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTgcgc                                                                  <  2:27829/150‑1 (MQ=255)
                                                                                       gCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAAc                                                              >  2:722106/1‑150 (MQ=255)
                                                                                        cGGGGGATAGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACa                                                             <  2:442475/150‑1 (MQ=255)
                                                                                        cGGGGGACGGCTCGGGTATGTGTGCGATTCTTTTACCGATCATTCCAGTGAATTTATCGATCAGCTGCACGATAAAATTATCGACCTTGAAGCTAATCTCCTTGATGAGCAAATTCCACCTCGTGGATTCCTGGCTCTGCTGCGCAAACa                                                             <  2:1012600/150‑1 (MQ=255)
                                                                                         gggggATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACaa                                                            <  2:139324/150‑1 (MQ=255)
                                                                                          ggggATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACaat                                                           >  2:835716/1‑150 (MQ=255)
                                                                                              aTGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACaattaat                                                       <  1:424342/150‑1 (MQ=255)
                                                                                               tGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACaattaatt                                                      >  2:1227549/1‑150 (MQ=255)
                                                                                                 tcTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATtgt                                                    <  1:190132/149‑1 (MQ=255)
                                                                                                         aTGTGTGCGAAGCCTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCAAGATAAAATTATCGACCTTGAAGATAATCTACTTGATCAGCAAATTCCACCGCGTGGAATCGAGGCTCTGCTGCGCAAACAATTAATTGTGATGCGta                                            >  2:1412386/1‑149 (MQ=255)
                                                                                                           gtgtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGc                                          <  1:500218/150‑1 (MQ=255)
                                                                                                            tgtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCt                                         >  2:1069374/1‑150 (MQ=255)
                                                                                                             gtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCta                                        >  2:1110596/1‑150 (MQ=255)
                                                                                                               gCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCtata                                      >  1:1450125/1‑150 (MQ=255)
                                                                                                                 gATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATg                                    <  1:495365/150‑1 (MQ=255)
                                                                                                                 gATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATg                                    >  1:1171671/1‑150 (MQ=255)
                                                                                                                  aTGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAAATGTGATGCGTCGCTATATgg                                   >  2:741930/1‑150 (MQ=255)
                                                                                                                     cGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGa                                                  >  2:413240/1‑132 (MQ=255)
                                                                                                                     cGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGa                                                  <  1:413240/132‑1 (MQ=255)
                                                                                                                     cGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCAc                                >  1:415074/1‑150 (MQ=255)
                                                                                                                      gTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTGATTGTGATGCGTCGCTATATGGCaac                               <  2:1505303/150‑3 (MQ=255)
                                                                                                                         gACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCa                            >  2:69472/1‑150 (MQ=255)
                                                                                                                             gATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGt                        >  1:1322643/1‑150 (MQ=255)
                                                                                                                              aTCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCATCAAATTCCACCTCGTGGTTTCCTGTCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCAGCGCAACGTg                       <  2:1318367/150‑1 (MQ=255)
                                                                                                                              aTCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTg                       <  2:415074/150‑1 (MQ=255)
                                                                                                                                cATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGAt                     >  2:1006525/1‑150 (MQ=255)
                                                                                                                                cATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGAt                     >  1:730631/1‑150 (MQ=255)
                                                                                                                                     cgataattttATCGTGCAGCTGCTCGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTa                >  1:1238338/8‑150 (MQ=255)
                                                                                                                                       ataattttATCGTGCAGCTGCTCGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATg              >  1:416852/6‑150 (MQ=255)
                                                                                                                                        tGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGc             <  1:305424/150‑1 (MQ=255)
                                                                                                                                         gAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCt            <  1:1318714/150‑1 (MQ=255)
                                                                                                                                           aTTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCg          <  2:1004499/150‑1 (MQ=255)
                                                                                                                                                  gAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGcc   >  1:189168/1‑150 (MQ=255)
                                                                                                                                                  gAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGcc   <  1:1006525/150‑1 (MQ=255)
                                                                                                                                                  gAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCACAACGTGATGTTTATGCTCGTCTTGcc   >  2:928329/1‑150 (MQ=255)
                                                                                                                                                   agcagcTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGCCa  >  2:864350/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGCCA  >  NC_000913/1408390‑1408686

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