Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,838,583 T→A 5.5% I112I (ATT→ATA yicL → EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,5830TA 5.5% 388.7 / 10.1 110I112I (ATT→ATAyicLEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (53/51);  new base A (3/3);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 110 total aligned reads displayed.

TTCTATTATTAACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTG  >  NC_000913/3838436‑3838729
                                                                                                                                                   |                                                                                                                                                  
ttCTATTATTAACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTa                                                                                                                                                   >  2:466851/1‑149 (MQ=255)
     ttattaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTAtcgtc                                                                                                                                              >  1:286512/1‑149 (MQ=255)
       attaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGc                                                                                                                                            >  3:198630/1‑149 (MQ=255)
       attaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCGCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTTCTCTCACCGACTATTATCGTCGc                                                                                                                                            >  3:345350/1‑149 (MQ=255)
        ttaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGTGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTAtc                                                                                                                                                 >  2:799558/1‑143 (MQ=255)
        ttaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGcc                                                                                                                                           >  3:73493/1‑149 (MQ=255)
        ttaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATAAtcgtcggt                                                                                                                                           <  2:780606/149‑3 (MQ=255)
        ttaACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATAAtcgtcggt                                                                                                                                           >  1:780606/1‑147 (MQ=255)
         taACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCt                                                                                                                                          <  3:276088/149‑1 (MQ=255)
          aaCAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTg                                                                                                                                         >  3:871106/1‑149 (MQ=255)
          aaCAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTg                                                                                                                                         >  4:41449/1‑149 (MQ=255)
             aaTCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGtt                                                                                                                                      >  3:764288/1‑149 (MQ=255)
             aaTCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGtt                                                                                                                                      >  4:74045/1‑149 (MQ=255)
             aaTCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATAAtcgtcggtga                                                                                                                                         <  4:883673/146‑5 (MQ=255)
             aaTCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATAAtcgtcggtga                                                                                                                                         >  3:883673/1‑142 (MQ=255)
             aaTCATAAAGAGGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGtt                                                                                                                                      >  3:727765/1‑149 (MQ=255)
                 aTAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCa                                                                                                                                  >  2:49468/1‑149 (MQ=255)
                  tAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCAc                                                                                                                                 <  2:607023/149‑1 (MQ=255)
                     aGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTgg                                                                                                                              <  1:306846/149‑1 (MQ=255)
                     aGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCGGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTgg                                                                                                                              >  4:811018/1‑149 (MQ=255)
                      gATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGt                                                                                                                             >  3:199899/1‑149 (MQ=255)
                      gATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGt                                                                                                                             <  2:664286/149‑1 (MQ=255)
                      gATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGt                                                                                                                             <  2:512479/149‑1 (MQ=255)
                      gATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGt                                                                                                                             <  4:199899/149‑1 (MQ=255)
                           cATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTa                                                                                                                        >  3:712185/1‑149 (MQ=255)
                             ttAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTaaa                                                                                                                      >  1:282081/1‑149 (MQ=255)
                               aGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATc                                                                                                                    >  4:222219/1‑149 (MQ=255)
                                 ccTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATcgc                                                                                                                  <  2:121459/149‑1 (MQ=255)
                                  ctgctgATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATc                                                                                                                    <  1:142939/146‑1 (MQ=255)
                                  ctgctgATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATc                                                                                                                    >  2:142939/1‑146 (MQ=255)
                                     ctgATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGccc                                                                                                              >  1:206467/1‑149 (MQ=255)
                                      tgATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCg                                                                                                             >  3:357144/1‑149 (MQ=255)
                                        aTTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCTACAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGCCTCTCACTGGTGCGTAAATCGCGCCCggg                                                                                                           >  2:67888/1‑149 (MQ=255)
                                                cGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCAtt                                                                                                       <  3:571241/145‑1 (MQ=255)
                                                cGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCAtt                                                                                                       >  4:571241/1‑145 (MQ=255)
                                                 gTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGt                                                                                                  >  2:532781/1‑149 (MQ=255)
                                                  tGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGtt                                                                                                 <  4:151299/149‑1 (MQ=255)
                                                          gcgcTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTgcgc                                                                                         >  2:434411/1‑149 (MQ=255)
                                                             cTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTAt                                                                                      >  1:552610/1‑149 (MQ=255)
                                                                aCTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTAtttt                                                                                   <  3:41449/149‑1 (MQ=255)
                                                                     aCAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACAt                                                                              >  3:139243/1‑149 (MQ=255)
                                                                        gCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGc                                                                           <  2:353532/149‑1 (MQ=255)
                                                                         cTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCt                                                                          <  2:930217/149‑1 (MQ=255)
                                                                         cTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCt                                                                          >  2:841005/1‑149 (MQ=255)
                                                                         cTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCt                                                                          >  2:216430/1‑149 (MQ=255)
                                                                           cACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTgg                                                                        >  3:262471/1‑149 (MQ=255)
                                                                            aCTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGt                                                                       >  1:854585/1‑149 (MQ=255)
                                                                            aCTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGt                                                                       <  1:434411/149‑1 (MQ=255)
                                                                            aCTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGt                                                                       >  3:274422/1‑149 (MQ=255)
                                                                            aCTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGAGCTGCGATTCCTCTCACCGACGATTATCGTCGCCTGGTTATCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGt                                                                       >  4:568172/1‑149 (MQ=255)
                                                                               tttttGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCgg                                                                    <  4:262471/149‑1 (MQ=255)
                                                                                ttttGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCggg                                                                   >  4:657070/1‑149 (MQ=255)
                                                                                  ttGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCACGGGCATTCTGGTTTTCTGCGCTATTTTGATAACGCTGGTCGGGAc                                                                 <  4:885957/149‑1 (MQ=255)
                                                                                    gCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACtt                                                               <  4:73493/149‑1 (MQ=255)
                                                                                      tAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACtttt                                                             <  1:841156/149‑1 (MQ=255)
                                                                                      tAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACtttt                                                             <  1:238843/149‑1 (MQ=255)
                                                                                       aaCCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACttttt                                                            <  1:470904/149‑1 (MQ=255)
                                                                                       aaCCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCAGCAATTCCTCTCACCGACGATTATCGTCACCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACttttt                                                            >  3:175606/1‑149 (MQ=255)
                                                                                        aCCATCGAAAAATCGAATGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTTGTCGGGGCtttttt                                                           >  2:843833/1‑149 (MQ=255)
                                                                                        aCCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACtttttt                                                           <  4:836327/149‑1 (MQ=255)
                                                                                         ccATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTTGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTa                                                          <  4:175606/149‑1 (MQ=255)
                                                                                         ccATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTa                                                          <  4:404787/149‑1 (MQ=255)
                                                                                           aTCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTAtt                                                        >  2:146454/1‑149 (MQ=255)
                                                                                           aTCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTAtt                                                        <  1:312519/149‑1 (MQ=255)
                                                                                           aTCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTAtt                                                        >  2:148567/1‑149 (MQ=255)
                                                                                            tCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTg                                                       >  4:773626/1‑149 (MQ=255)
                                                                                            tCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTg                                                       <  2:118598/149‑1 (MQ=255)
                                                                                            tCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTg                                                       <  1:350745/149‑1 (MQ=255)
                                                                                              gAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGt                                                     >  1:363433/1‑149 (MQ=255)
                                                                                              gAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGt                                                     <  2:374171/149‑1 (MQ=255)
                                                                                              gAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGt                                                     <  1:306373/149‑1 (MQ=255)
                                                                                               aaaaaTCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTg                                                       <  2:186777/146‑1 (MQ=255)
                                                                                               aaaaaTCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTg                                                       >  1:186777/1‑146 (MQ=255)
                                                                                               aaaaaTCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTg                                                    >  2:362875/1‑149 (MQ=255)
                                                                                                aaaaTCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGa                                                   <  4:500020/149‑1 (MQ=255)
                                                                                                aaaaTCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGa                                                   >  3:923561/1‑149 (MQ=255)
                                                                                                   aTCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTCACATCGCTGGTCGGGACTTTTTTATTGGTGacac                                                <  1:466851/149‑1 (MQ=255)
                                                                                                       aaCGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACgg                                            <  3:97068/149‑1 (MQ=255)
                                                                                                       aaCGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACgg                                            <  1:535803/149‑1 (MQ=255)
                                                                                                          gCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTaa                                         <  3:106099/149‑1 (MQ=255)
                                                                                                            aGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGt                                           <  2:779028/145‑1 (MQ=255)
                                                                                                            aGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGt                                           >  1:779028/1‑145 (MQ=255)
                                                                                                             gCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATcc                                      >  3:817594/1‑149 (MQ=255)
                                                                                                               cACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGa                                    <  4:923561/149‑1 (MQ=255)
                                                                                                                aCGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGAc                                   <  2:798770/149‑1 (MQ=255)
                                                                                                                            ctgcAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGAt                       <  4:357144/149‑1 (MQ=255)
                                                                                                                            ctgcAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGAt                       <  2:661996/149‑1 (MQ=255)
                                                                                                                            ctgcAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGAt                       <  1:146454/149‑1 (MQ=255)
                                                                                                                                  ttCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTcc                 <  4:727765/149‑1 (MQ=255)
                                                                                                                                     ctctCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGc              >  2:7564/1‑149 (MQ=255)
                                                                                                                                       ctcACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCcg            >  2:733874/1‑149 (MQ=255)
                                                                                                                                        tcACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGAc                                   <  3:862625/125‑1 (MQ=255)
                                                                                                                                        tcACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGAc                                   >  4:862625/1‑125 (MQ=255)
                                                                                                                                         aggcgacgaTAATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGAtct                     >  2:904422/4‑138 (MQ=255)
                                                                                                                                         aggcgacgaTAATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGAtct                     <  1:904422/135‑1 (MQ=255)
                                                                                                                                             gacgaTTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGt                                                     <  2:347600/102‑1 (MQ=255)
                                                                                                                                                gaTTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCt   >  3:111896/1‑149 (MQ=255)
                                                                                                                                                 aTTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTg  <  4:139243/149‑1 (MQ=255)
                                                                                                                                                 aTTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTg  >  1:594852/1‑149 (MQ=255)
                                                                                                                                                 aTTATCGTCGCCTGGTGCTCGCTGGTGCGTAAGGCGCACCGGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGTTTGTCGGGACTTTTTTAATGCTGAGACACGGTAATCCGACGTCATTATCGATCGCTCCTGCCTCGTTGTTCTg  <  2:355433/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCTATTATTAACAATCATAAAGATGCCATTAGCCTGCTGATTTTTTCCGTGGTTGGCGCGCTAACTGTACAGCTCACTTTTTTGCTAACCATCGAAAAATCGAACGCAGCCACGGCAACGGTGCTGCAATTCCTCTCACCGACGATTATCGTCGCCTGGTTCTCACTGGTGCGTAAATCGCGCCCGGGCATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTG  >  NC_000913/3838436‑3838729

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