Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 929,684 A→T 7.8% W427R (TGG→AGG)  cydD ← glutathione/cysteine ABC transporter export permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913929,6840AT 7.8% 381.4 / 16.8 110W427R (TGG→AGG) cydDglutathione/cysteine ABC transporter export permease/ATPase
Reads supporting (aligned to +/- strand):  ref base A (45/56);  new base T (4/5);  total (49/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

ATCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTT  >  NC_000913/929535‑929833
                                                                                                                                                     |                                                                                                                                                     
aTCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCa                                                                                                                                                       <  2:1159786/150‑1 (MQ=255)
aTCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCa                                                                                                                                                       <  2:151694/150‑1 (MQ=255)
 tCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAg                                                                                                                                                      <  2:1646039/150‑1 (MQ=255)
 tCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAg                                                                                                                                                      >  2:1374094/1‑150 (MQ=255)
 tCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAg                                                                                                                                                      <  2:796616/150‑1 (MQ=255)
   aaCGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGga                                                                                                                                                    <  2:1441040/150‑1 (MQ=255)
   aaCGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGga                                                                                                                                                    <  2:1370670/150‑1 (MQ=255)
     cGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGgaga                                                                                                                                                  <  2:1233787/150‑1 (MQ=255)
        cTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTGGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGGACTTGTGGGTTTTGCCCAACCCAGGAGAGAt                                                                                                                                               >  2:234490/1‑150 (MQ=255)
         ttGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATg                                                                                                                                              >  1:1635532/1‑150 (MQ=255)
              ggaggaGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTAc                                                                                                                                         >  2:420655/1‑150 (MQ=255)
                aggagCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGc                                                                                                                                       >  2:881758/1‑150 (MQ=255)
                         aGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCt                                                                                                                              <  2:1591595/150‑1 (MQ=255)
                          gAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTg                                                                                                                             <  1:634690/150‑1 (MQ=255)
                           aaaCTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTgg                                                                                                                            <  2:513160/150‑1 (MQ=255)
                             aCTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTg                                                                                                                          >  2:1440179/1‑150 (MQ=255)
                              cTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCTGGAGAGATGTTTACGCCATGATTCTGGTGa                                                                                                                         >  2:401135/1‑150 (MQ=255)
                              cTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCTGGAGAGATGTTTACGCCATGATTCTGGTGa                                                                                                                         >  2:1265848/1‑150 (MQ=255)
                              cTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCTGGAGAGATGTTTACGCCATGATTCTGGTGa                                                                                                                         >  2:1687893/1‑150 (MQ=255)
                              cTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCTGGAGAGATGTTTACGCCATGATTCTGGTGa                                                                                                                         >  2:163699/1‑150 (MQ=255)
                                 gCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATaa                                                                                                                      >  1:1303629/1‑150 (MQ=255)
                                    gACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATc                                                                                                                   <  1:1361534/150‑1 (MQ=255)
                                    gACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATc                                                                                                                   <  1:1361491/150‑1 (MQ=255)
                                          ggCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTaa                                                                                                             <  2:1721535/150‑1 (MQ=255)
                                          ggCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTaa                                                                                                             <  2:869500/150‑1 (MQ=255)
                                            cGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAAtt                                                                                                           >  2:64380/1‑150 (MQ=255)
                                              ttATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCt                                                                                                         >  1:448680/1‑150 (MQ=255)
                                              ttATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGAAGTAAGTAATCCCTAAATCTTGATGCCGGTAATAGAGGGTTTTGCCCAACCCAGGAGAGATTTTTACGCCATTATTCTGGTGATAAATCTCGTAATTCt                                                                                                         >  2:1442359/1‑150 (MQ=255)
                                                 tCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATc                                                                                                      >  1:206133/1‑150 (MQ=255)
                                                     gcgcTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGt                                                                                                  <  2:1503291/150‑1 (MQ=255)
                                                      cgcTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGtt                                                                                                 >  1:1030787/1‑150 (MQ=255)
                                                      cgcTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGtt                                                                                                 >  1:1680140/1‑150 (MQ=255)
                                                            ttGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCg                                                                                           <  2:737649/150‑1 (MQ=255)
                                                             tGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGt                                                                                          >  1:1122286/1‑150 (MQ=255)
                                                              gTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTa                                                                                         <  1:1412484/150‑1 (MQ=255)
                                                               tAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAg                                                                                        <  1:951211/150‑1 (MQ=255)
                                                                      tGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATccc                                                                                 <  2:582537/150‑1 (MQ=255)
                                                                       gTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCt                                                                                <  1:1287103/150‑1 (MQ=255)
                                                                        ttCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTg                                                                               >  1:164987/1‑150 (MQ=255)
                                                                        ttCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTg                                                                               <  2:605777/150‑1 (MQ=255)
                                                                         tCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGa                                                                              <  1:1193882/150‑1 (MQ=255)
                                                                           gCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGata                                                                            >  1:937774/1‑150 (MQ=255)
                                                                           gCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGata                                                                            >  1:937820/1‑150 (MQ=255)
                                                                              ggCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATga                                                                         <  1:639392/150‑1 (MQ=255)
                                                                                cATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATgaga                                                                       <  1:471838/150‑1 (MQ=255)
                                                                                 aTCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATgagag                                                                      >  2:1692167/1‑150 (MQ=255)
                                                                                 aTCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATgagag                                                                      >  2:1692165/1‑150 (MQ=255)
                                                                                      gTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATGAATCGCGGACTTCTCTCCCGTGGATTCGTAGCGATCCCTGAGATgagagacagc                                                                 >  2:889611/1‑146 (MQ=255)
                                                                                      gTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAAc                                                                 >  2:889733/1‑150 (MQ=255)
                                                                                      gTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAAc                                                                 >  1:1706181/1‑150 (MQ=255)
                                                                                        cgcgCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCa                                                               <  2:709962/150‑1 (MQ=255)
                                                                                        cgcgCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCa                                                               <  1:1399789/150‑1 (MQ=255)
                                                                                             cAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAg                                                          >  1:1112933/1‑150 (MQ=255)
                                                                                                tagtaCGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAgcgc                                                       <  2:1714726/150‑1 (MQ=255)
                                                                                                tagtaCGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAgcgc                                                       <  2:862685/150‑1 (MQ=255)
                                                                                                   taCGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGtt                                                    <  1:205637/150‑1 (MQ=255)
                                                                                                   taCGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGtt                                                    <  1:205718/150‑1 (MQ=255)
                                                                                                      gTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTcag                                                 <  2:801617/150‑1 (MQ=255)
                                                                                                      gTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTcag                                                 <  1:1692165/150‑1 (MQ=255)
                                                                                                      gTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTcag                                                 <  1:1692167/150‑1 (MQ=255)
                                                                                                      gTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTcag                                                 <  1:295342/150‑1 (MQ=255)
                                                                                                      gTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTcag                                                 <  2:801620/150‑1 (MQ=255)
                                                                                                       ttATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTcagc                                                >  2:1602858/1‑150 (MQ=255)
                                                                                                            ccGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGa                                           <  2:1635532/150‑1 (MQ=255)
                                                                                                             cGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAg                                          >  2:226152/1‑150 (MQ=255)
                                                                                                             cGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAg                                          >  1:1471309/1‑150 (MQ=255)
                                                                                                             cGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAg                                          >  1:1270642/1‑150 (MQ=255)
                                                                                                                aaTGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCTGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCtt                                       <  1:163699/150‑1 (MQ=255)
                                                                                                                aaTGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCTGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCtt                                       <  1:401135/150‑1 (MQ=255)
                                                                                                                 aTGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCttt                                      <  2:1496682/150‑1 (MQ=255)
                                                                                                                    ttGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTAc                                   <  1:1344817/150‑1 (MQ=255)
                                                                                                                    ttGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTAc                                   <  1:1241638/150‑1 (MQ=255)
                                                                                                                    gtgctgcCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTAc                                   <  1:454461/149‑1 (MQ=255)
                                                                                                                         gcCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGaa                              <  1:407142/150‑1 (MQ=255)
                                                                                                                          ccGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAAc                             >  2:558686/1‑150 (MQ=255)
                                                                                                                           cGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAAcc                            >  1:642107/1‑150 (MQ=255)
                                                                                                                           cGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAAcc                            >  1:642170/1‑150 (MQ=255)
                                                                                                                           cGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAAcc                            >  1:751904/1‑150 (MQ=255)
                                                                                                                             gTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGc                          >  2:1000145/1‑150 (MQ=255)
                                                                                                                                  tgtgGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGa                     <  2:1121324/150‑1 (MQ=255)
                                                                                                                                  tgtgGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGa                     <  1:1593027/150‑1 (MQ=255)
                                                                                                                                   gtgGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGac                    <  1:756816/150‑1 (MQ=255)
                                                                                                                                   gtgGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGac                    <  1:452177/150‑1 (MQ=255)
                                                                                                                                    ggggTTTTGCCCTACCCTGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTTATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGtcc                   <  1:234490/149‑1 (MQ=255)
                                                                                                                                       gTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGaccaac                >  1:1716537/1‑150 (MQ=255)
                                                                                                                                       gTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGaccaac                <  1:64380/150‑1 (MQ=255)
                                                                                                                                        ttttGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGaccaacc               >  2:612638/1‑150 (MQ=255)
                                                                                                                                        ttttGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGaccaacc               <  1:1455738/150‑1 (MQ=255)
                                                                                                                                        ttttGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGaccaacc               <  1:819466/150‑1 (MQ=255)
                                                                                                                                          ttGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATCAATCGCGTAATTCTCGCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGAccaaccaa             >  2:239732/1‑150 (MQ=255)
                                                                                                                                          ttGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGAccaaccaa             <  2:1122286/150‑1 (MQ=255)
                                                                                                                                               cAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGTGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCg        >  1:1157910/1‑150 (MQ=255)
                                                                                                                                               cAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCACCACCg        >  1:1652292/1‑150 (MQ=255)
                                                                                                                                               cAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCg        >  1:1478104/1‑150 (MQ=255)
                                                                                                                                               cAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCg        >  1:177485/1‑150 (MQ=255)
                                                                                                                                               cAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCg        >  1:622225/1‑150 (MQ=255)
                                                                                                                                                     aGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGtt  >  2:42763/1‑150 (MQ=255)
                                                                                                                                                     aGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGtt  <  2:1550184/150‑1 (MQ=255)
                                                                                                                                                     aGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGtt  >  2:1673233/1‑150 (MQ=255)
                                                                                                                                                     aGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGAGTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGtt  >  2:1625175/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCAACGCCTTGTGGGAGGAGCGGTAGAAACTCGCTGACCCAGGCGTTATCCAGCGCTGCTTGTAATTCTTGTTCGCTGGCATCAGGTCGCGCCAGTAGTACGTTATCCCGCAATGTTGCTGCCGGTAATTGTGGGTTTTGCCCAACCCAGGAGAGATGTTTACGCCATGATTCTGGTGATAAATCGCGTAATTCTATCCCGTTGATTCGTAGCGATCCCTGATATGAGAGAAAACCAGAAAGCGCGTTCAGCAGTGAGCTTTTACCTGAACCGCTGCGACCAACCAACACCGCACGTT  >  NC_000913/929535‑929833

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