Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,741,991 T→C 5.3% I39I (ATT→ATC yfiR → putative periplasmic inhibitor of YfiN activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,741,9910TC 5.3% 261.3 / 3.3 77I39I (ATT→ATCyfiRputative periplasmic inhibitor of YfiN activity
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base C (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCATTTCGAGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTCAACAAGAGGC  >  NC_000913/2741843‑2742134
                                                                                                                                                    |                                                                                                                                               
gCATTTCGAGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTg                                                                                                                                                >  2:497416/1‑150 (MQ=255)
  aTTTCGAGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTc                                                                                                                                              <  2:1368975/150‑1 (MQ=255)
  aTTTCGAGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTc                                                                                                                                              <  2:1368974/150‑1 (MQ=255)
        aGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTac                                                                                                                                        >  2:1380739/1‑150 (MQ=255)
        aGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTac                                                                                                                                        >  2:5402/1‑150 (MQ=255)
        aGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATAGTCAGTTac                                                                                                                                        >  1:830134/1‑150 (MQ=255)
         gCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTaca                                                                                                                                       >  2:812519/1‑150 (MQ=255)
         gCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTaca                                                                                                                                       >  2:812436/1‑150 (MQ=255)
          cGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTAAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTacac                                                                                                                                      >  1:436001/1‑150 (MQ=255)
             aGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCg                                                                                                                                   <  2:1002909/150‑1 (MQ=255)
             aGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCg                                                                                                                                   <  2:1002907/150‑1 (MQ=255)
             aGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCg                                                                                                                                   >  1:978764/1‑150 (MQ=255)
              gATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGc                                                                                                                                  <  2:1091962/150‑1 (MQ=255)
              gATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGc                                                                                                                                  <  1:454569/150‑1 (MQ=255)
               aTGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCt                                                                                                                                 <  1:719622/150‑1 (MQ=255)
                tggtggGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTg                                                                                                                                <  2:442209/150‑1 (MQ=255)
                tggtggGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTg                                                                                                                                <  2:442204/150‑1 (MQ=255)
                  gtggGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGc                                                                                                                              <  1:129107/150‑1 (MQ=255)
                   tggGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGcc                                                                                                                             >  1:919979/1‑150 (MQ=255)
                           ccTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTAt                                                                                                                     >  2:1259543/1‑150 (MQ=255)
                           ccTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTAt                                                                                                                     >  2:1259602/1‑150 (MQ=255)
                               aaTGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAgg                                                                                                                 <  1:259027/150‑1 (MQ=255)
                               aaTGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAgg                                                                                                                 <  1:380812/150‑1 (MQ=255)
                                 tGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGc                                                                                                               <  2:1333541/150‑1 (MQ=255)
                                 tGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGc                                                                                                               <  1:1218051/150‑1 (MQ=255)
                                    gTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGccac                                                                                                            >  1:341397/1‑150 (MQ=255)
                                         tctcACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAg                                                                                                       <  2:859121/150‑1 (MQ=255)
                                           tcACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCCAAGtt                                                                                                     >  1:81965/1‑150 (MQ=255)
                                                 aCTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTAt                                                                                               >  2:644733/1‑150 (MQ=255)
                                                   tGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCt                                                                                             <  2:585495/150‑1 (MQ=255)
                                                     ttCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCttt                                                                                           >  1:1197371/1‑150 (MQ=255)
                                                        ttACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCa                                                                                        <  2:945390/150‑1 (MQ=255)
                                                           cTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATCGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCt                                                                                     >  2:1156069/1‑150 (MQ=255)
                                                           cTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATCGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCt                                                                                     >  2:1124246/1‑150 (MQ=255)
                                                            ttCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCtt                                                                                    >  1:636047/1‑150 (MQ=255)
                                                            ttCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCtt                                                                                    >  2:968890/1‑150 (MQ=255)
                                                             tCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTc                                                                                   >  1:1064007/1‑150 (MQ=255)
                                                             tCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTc                                                                                   <  2:218884/150‑1 (MQ=255)
                                                             tCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTc                                                                                   <  2:816539/150‑1 (MQ=255)
                                                             tCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTc                                                                                   >  1:1243274/1‑150 (MQ=255)
                                                             tCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTc                                                                                   <  2:816528/150‑1 (MQ=255)
                                                               tGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCAc                                                                                 <  2:752167/150‑1 (MQ=255)
                                                                gATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACg                                                                                >  2:870892/1‑150 (MQ=255)
                                                                gATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACg                                                                                >  1:798904/1‑150 (MQ=255)
                                                                    ctcCTGCTGACGGGAGCACCACTACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGtttt                                                                            <  2:65605/150‑1 (MQ=255)
                                                                    ctcCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAACGTTATGTATCTTTTCATCTTCACGtttt                                                                            <  2:1344378/150‑1 (MQ=255)
                                                                    ctcCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGtttt                                                                            <  2:65579/150‑1 (MQ=255)
                                                                     tcCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTa                                                                           <  2:980082/150‑1 (MQ=255)
                                                                     tcCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTa                                                                           >  1:1010106/1‑150 (MQ=255)
                                                                     tcCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTa                                                                           >  1:1010196/1‑150 (MQ=255)
                                                                        tgctgACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTa                                                                        >  1:1298811/1‑150 (MQ=255)
                                                                             aCGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGcc                                                                   >  1:777081/1‑150 (MQ=255)
                                                                               gggAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCt                                                                 >  1:328003/1‑150 (MQ=255)
                                                                                     aCCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAaga                                                           <  1:1011463/150‑1 (MQ=255)
                                                                                     aCCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAaga                                                           >  2:635564/1‑150 (MQ=255)
                                                                                                cTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGa                                                >  1:1188635/1‑150 (MQ=255)
                                                                                                cTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGa                                                <  2:678713/150‑1 (MQ=255)
                                                                                                    ggAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATc                                            <  1:1375042/150‑1 (MQ=255)
                                                                                                           tCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCt                                     <  2:1027336/150‑1 (MQ=255)
                                                                                                           tCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCt                                     <  2:1027325/150‑1 (MQ=255)
                                                                                                            cTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCtt                                    >  2:366191/1‑150 (MQ=255)
                                                                                                             tGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTa                                   >  2:308275/1‑150 (MQ=255)
                                                                                                             tGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTa                                   >  2:308283/1‑150 (MQ=255)
                                                                                                                  ttGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGc                              <  1:782844/150‑1 (MQ=255)
                                                                                                                    gCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCt                            <  1:616929/150‑1 (MQ=255)
                                                                                                                    gCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCt                            <  2:1038041/150‑1 (MQ=255)
                                                                                                                      gAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGt                          >  2:169891/1‑150 (MQ=255)
                                                                                                                      gAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGt                          <  2:215691/150‑1 (MQ=255)
                                                                                                                          aaCGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGtcatc                      <  2:1400990/150‑1 (MQ=255)
                                                                                                                               gCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTca                 <  2:207527/150‑1 (MQ=255)
                                                                                                                               gCGGATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTca                 <  2:207442/150‑1 (MQ=255)
                                                                                                                                    tgatgGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTc            >  1:23088/1‑150 (MQ=255)
                                                                                                                                      atgGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTcaa          >  2:1330752/1‑150 (MQ=255)
                                                                                                                                           ttCCGGTATCGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTCAACAaga     <  1:1156069/150‑1 (MQ=255)
                                                                                                                                           ttCCGGTATCGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTCAACAaga     <  1:1124246/150‑1 (MQ=255)
                                                                                                                                             ccGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTCAACAAGAgg   >  1:1218768/1‑150 (MQ=255)
                                                                                                                                              cGGTATTGTCAGTTACACCTGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTCAACAAGAGGc  >  2:726766/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCATTTCGAGCGAAGATGGTGGGGATCCCTGAATGCGTTTTTCTCACCGACTGTTCTTACTTCTGATTCTCCTGCTGACGGGAGCACCAATACTCGCTCAGGAGCCTTCTGACGTTGCGAAAAACGTGCGTATGATGGTTTCCGGTATTGTCAGTTACACCCGCTGGCCTGCGCTATCAGGGCCACCAAAGTTATGTATCTTTTCATCTTCACGTTTTAGTACCGCCCTACAAGAGAATGCCGCGACATCTTTGCCTTATCTGCCTGTCATCATTCACACTCAACAAGAGGC  >  NC_000913/2741843‑2742134

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: