Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,453,330 A→G N254D (AAC→GAC)  ytfT → putative sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,453,3300AG99.1% 409.3 / ‑4.1 115N254D (AAC→GAC) ytfTputative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (58/56);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 117 total aligned reads displayed.

TCGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCGG  >  NC_000913/4453184‑4453477
                                                                                                                                                  |                                                                                                                                                   
tCGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCcgaca                                                                                                                                                  >  4:476608/1‑150 (MQ=255)
 cGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCcgacaa                                                                                                                                                 <  1:1048300/150‑1 (MQ=255)
  gAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACg                                                                                                                                               <  2:1364559/151‑1 (MQ=255)
    aGCCGTTGGTATCAACATTCGGGCGGCAAATAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCt                                                                                                                                             >  1:561132/1‑151 (MQ=255)
    aGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCt                                                                                                                                             >  2:999395/1‑151 (MQ=255)
     gCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTg                                                                                                                                            <  2:671394/151‑1 (MQ=255)
      ccGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTgg                                                                                                                                           >  2:1433209/1‑151 (MQ=255)
        gTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGt                                                                                                                                         <  2:1033385/151‑1 (MQ=255)
             tATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATgg                                                                                                                                    >  2:795990/1‑151 (MQ=255)
              aTCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgc                                                                                                                                   >  1:518742/1‑151 (MQ=255)
                cAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctg                                                                                                                                 >  2:1201886/1‑151 (MQ=255)
                cAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctg                                                                                                                                 >  1:965546/1‑151 (MQ=255)
                 aaCATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctgg                                                                                                                                >  4:216334/1‑151 (MQ=255)
                  aCATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctgga                                                                                                                               <  3:112402/151‑1 (MQ=255)
                   cATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctggag                                                                                                                              <  4:6225/151‑1 (MQ=255)
                    aTTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctggagc                                                                                                                             >  1:811840/1‑151 (MQ=255)
                         ggcggcAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGAc                                                                                                                        >  1:562279/1‑151 (MQ=255)
                            ggcAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTTTTGAGCTGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGTGGATATTCGCGGTGCCGTTGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGca                                                                                                                     <  1:1699502/151‑2 (MQ=255)
                              cAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCAt                                                                                                                   <  2:1671635/151‑1 (MQ=255)
                              cAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCAt                                                                                                                   >  2:1062947/1‑151 (MQ=255)
                                   aaTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTcg                                                                                                              <  1:1020262/151‑1 (MQ=255)
                                    aTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTcgc                                                                                                             >  1:920807/1‑151 (MQ=255)
                                    aTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTcgc                                                                                                             <  1:1209303/151‑1 (MQ=255)
                                      gCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTTTGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGGGCTGGACGCCATTCTCGCgg                                                                                                           >  1:266569/1‑151 (MQ=255)
                                      gCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTcgcg                                                                                                            <  1:472005/150‑1 (MQ=255)
                                       ccGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCggt                                                                                                          <  3:478509/151‑1 (MQ=255)
                                       ccGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCgg                                                                                                           >  1:348163/1‑150 (MQ=255)
                                        cGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACATCGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCggt                                                                                                          <  2:1129012/150‑1 (MQ=255)
                                         ggggTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCggtgg                                                                                                        >  2:1569057/1‑151 (MQ=255)
                                              aaaCACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCggtggt                                                                                                       >  1:826306/1‑147 (MQ=255)
                                                acacGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTgg                                                                                                 <  2:1113371/151‑1 (MQ=255)
                                                   cgcgAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGTCATTATCGTGGCGGCGGATATTCGCGGTGTCGATGCCGACAACGCTGGGTTATGGCTGGAGGCTGACGGCCTTCTCGCGGTGGGGATTggggg                                                                                              >  1:26980/1‑151 (MQ=255)
                                                     cgAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGgcggcg                                                                                            >  2:764092/1‑151 (MQ=255)
                                                     cgAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGgcggcg                                                                                            <  2:1281960/151‑1 (MQ=255)
                                                      gAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGcggcgg                                                                                           >  2:1353385/1‑151 (MQ=255)
                                                      gAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGcggcgg                                                                                           <  2:1296242/151‑1 (MQ=255)
                                                      gAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGcggcgg                                                                                           >  3:382257/1‑151 (MQ=255)
                                                      gAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGGGGTGATTGGcggcgg                                                                                           >  2:119004/1‑151 (MQ=255)
                                                        atcatcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGAt                                                                                         >  3:296529/1‑151 (MQ=255)
                                                         tcatcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATc                                                                                        >  2:1655325/1‑151 (MQ=255)
                                                          catcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATc                                                                                        <  1:160894/150‑1 (MQ=255)
                                                             cGTCATGCTCACTTATGTGTTGAGCGAGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTg                                                                                    <  2:1135250/151‑1 (MQ=255)
                                                               tCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGGGATTGGCGGCGGATCGCTg                                                                                    >  4:311743/1‑149 (MQ=255)
                                                                cATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATg                                                                                 <  2:614578/151‑1 (MQ=255)
                                                                cATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATg                                                                                 <  1:789458/151‑1 (MQ=255)
                                                                 aTGCTCACTTAGGGGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGGTATTCGCGGTGCCGAGGCCGACAACGCTGGTGTATGGCTGGAGCTGGGCGGCCTTCTGGCGGTGGTGATTGGCGGCGGAAAGCTGATgg                                                                                >  1:729072/1‑151 (MQ=255)
                                                                  tGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGATATGGCTGGAGCTGGACGCCATTCTCGCGGTGGGGATTGGCGGCGGCGCGCTGATggg                                                                               >  1:1069361/1‑151 (MQ=255)
                                                                   gCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGGGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGGTCGGTGCTggg                                                                               >  1:1019664/1‑150 (MQ=255)
                                                                   gCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGCTCCCTGATggg                                                                               >  1:13113/1‑150 (MQ=255)
                                                                     tCACTTATGGGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCTGTTATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGGGGATCGCTGGTgggcgg                                                                            >  1:1587029/1‑151 (MQ=255)
                                                                      cACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATgggcggg                                                                           >  3:257776/1‑151 (MQ=255)
                                                                          tATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGgcgc                                                                        <  2:192175/150‑1 (MQ=255)
                                                                            tgtgTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCttt                                                                     >  2:128655/1‑151 (MQ=255)
                                                                             gtgtTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCttt                                                                     <  2:1376533/150‑1 (MQ=255)
                                                                                ttGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGGGCTTTAAcc                                                                 >  1:1201153/1‑151 (MQ=255)
                                                                                ttGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTTAcc                                                                 >  1:1510875/1‑151 (MQ=255)
                                                                                ttGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAAcc                                                                 >  3:379963/1‑151 (MQ=255)
                                                                                 tGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAAcc                                                                 <  1:1721337/150‑1 (MQ=255)
                                                                                        gcggtgtgCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGAGAACGTTGGGTTATGGCTGGAGCTGGATGCCATTCTCGCGGTGGTGATTGGCGGCGGATCTCTGATGGGCTGGCTCGTTAACCTGCTAttt                                                         <  1:892214/148‑1 (MQ=255)
                                                                                         cTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACttt                                                        >  1:5115/1‑151 (MQ=255)
                                                                                               gcggcgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCggtgg                                                  <  4:594142/151‑1 (MQ=255)
                                                                                               gcggcgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCggtg                                                   <  2:998231/150‑1 (MQ=255)
                                                                                                cggcgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCggtggt                                                 <  2:216298/151‑1 (MQ=255)
                                                                                                    gATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTgggg                                             >  4:520598/1‑151 (MQ=255)
                                                                                                        gCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGgcgc                                         <  4:552201/151‑1 (MQ=255)
                                                                                                          gggCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCt                                        <  2:1292912/150‑1 (MQ=255)
                                                                                                          gggCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTg                                       <  3:461742/151‑1 (MQ=255)
                                                                                                           ggCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGTCGGATCGCTGATGGGCGGGGGCTTTAACCTGCTCCTTTCGGTGGTGGGGGCGCTGa                                      <  4:53557/151‑1 (MQ=255)
                                                                                                               ttATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGattat                                  >  2:950928/1‑151 (MQ=255)
                                                                                                               ttATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCCCGCGGGGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGTTTCTTTCGGCGGTGGGGGCGCTGttt                                    >  1:235072/1‑149 (MQ=255)
                                                                                                                tATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGattat                                  <  4:394557/150‑1 (MQ=255)
                                                                                                                 aTCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTc                                >  2:1048117/1‑151 (MQ=255)
                                                                                                                   cGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTc                                >  1:684255/1‑149 (MQ=255)
                                                                                                                   cGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCa                               <  2:90878/150‑1 (MQ=255)
                                                                                                                   cGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCa                               <  1:503637/150‑1 (MQ=255)
                                                                                                                    gTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgg                             <  1:996970/151‑1 (MQ=255)
                                                                                                                    gTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgg                             <  3:6878/151‑1 (MQ=255)
                                                                                                                     tGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAggg                            >  2:965723/1‑151 (MQ=255)
                                                                                                                      ggcggcGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgggg                           <  1:1351978/151‑1 (MQ=255)
                                                                                                                      ggcggcGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgggg                           <  1:128125/151‑1 (MQ=255)
                                                                                                                        cggcggATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGAt                         <  3:377872/151‑1 (MQ=255)
                                                                                                                         ggcggATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATg                        <  2:1630022/151‑1 (MQ=255)
                                                                                                                             gATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTCCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAaca                    <  1:125309/151‑1 (MQ=255)
                                                                                                                               tatTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTGCTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAacac                   <  1:1263271/150‑1 (MQ=255)
                                                                                                                                atTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCg                 <  2:175554/151‑1 (MQ=255)
                                                                                                                                atTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCg                 <  4:389091/151‑1 (MQ=255)
                                                                                                                                atTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCg                 <  4:236830/151‑1 (MQ=255)
                                                                                                                                 ttCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCgg                >  1:1364435/1‑151 (MQ=255)
                                                                                                                                 ttCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGGGGGGGCGCTGATTATTCAGGGGATGAACACCgg                >  1:1292303/1‑151 (MQ=255)
                                                                                                                                    gcgGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAt             <  3:44027/151‑1 (MQ=255)
                                                                                                                                      ggTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAttt           <  1:664232/151‑1 (MQ=255)
                                                                                                                                       gTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGCCGCCATTCTCGAGGTGGTGATGGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAttt           <  3:596525/150‑1 (MQ=255)
                                                                                                                                          ccGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCt       >  2:804156/1‑151 (MQ=255)
                                                                                                                                           cGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCtt      >  1:1263396/1‑151 (MQ=255)
                                                                                                                                            gATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCttt     >  2:817281/1‑151 (MQ=255)
                                                                                                                                            gATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCtt      <  1:1568897/150‑1 (MQ=255)
                                                                                                                                             aTGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTc    <  3:425931/151‑1 (MQ=255)
                                                                                                                                             aTGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTc    >  2:1358143/1‑151 (MQ=255)
                                                                                                                                              tGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCg   >  3:279232/1‑151 (MQ=255)
                                                                                                                                               gCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCgg  >  2:1246685/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCGG  >  NC_000913/4453184‑4453477

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