Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,332,8230AG50.8% 119.3 / 192.4 132pseudogene (1580/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (16/49);  new base G (67/0);  total (83/49)
Fisher's exact test for biased strand distribution p-value = 1.39e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-13

Only 100 of 134 total aligned reads displayed.

GCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCG  >  NC_000913/2332675‑2332970
                                                                                                                                                     |                                                                                                                                                        
 cTATTACTGTATTGTGCGGCAGTACTTAATGGGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGgaga                                                                                                                                                       <  2:236744/150‑1 (MQ=255)
 cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAt                                                                                                                                                      >  2:1642499/1‑151 (MQ=255)
   aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGCGACGGTTAACAAATAGCGGCTCGGTTTATCTGCGGCGGGGGGAttt                                                                                                                                                    >  1:1052590/1‑151 (MQ=255)
   aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATTGGCGGCGGGGGGAttt                                                                                                                                                    >  2:7816/1‑151 (MQ=255)
   aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGttt                                                                                                                                                    >  2:990216/1‑151 (MQ=255)
   aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAttt                                                                                                                                                    >  2:1153812/1‑151 (MQ=255)
   aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAAACCCCGTTCAATGGGGTCCCTTTTGGTGGTGTTGACCCGTACCGCCCCGCCGCCCCGGCCGGGACGGTTTACCAAATACGGGCCGGGTTTTTCGGCGGGGGGGGGGttt                                                                                                                                                    >  1:872833/1‑151 (MQ=255)
    ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGGGATTTa                                                                                                                                                   >  2:44177/1‑151 (MQ=255)
    ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGtttt                                                                                                                                                   >  1:254878/1‑150 (MQ=255)
    ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTa                                                                                                                                                   >  2:692395/1‑151 (MQ=255)
    ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGATCACAAATAGCGGCTCGGCTTATAGGCCGCGGGGAGAgttt                                                                                                                                                   >  1:1488722/1‑150 (MQ=255)
       cTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACg                                                                                                                                                <  1:1620399/151‑1 (MQ=255)
        tGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGc                                                                                                                                               <  2:35351/151‑1 (MQ=255)
         gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGc                                                                                                                                               >  1:970196/1‑150 (MQ=255)
         gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGc                                                                                                                                               >  2:1243398/1‑150 (MQ=255)
         gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTTTCGGCGGCGGGGGGATTTTAccc                                                                                                                                               >  2:68864/1‑148 (MQ=255)
              gtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTAAGCcccaag                                                                                                                                         >  2:516048/1‑149 (MQ=255)
              gtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGCTTATCCGCCGCGGGGGGACTTAACGCCCCAt                                                                                                                                          >  1:549396/1‑150 (MQ=255)
                 cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCATTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAc                                                                                                                                       <  2:281963/150‑1 (MQ=255)
                 cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAc                                                                                                                                       <  1:783471/150‑1 (MQ=255)
                   gCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGATTAACGCcacaatgaca                                                                                                                                     >  1:1044132/1‑144 (MQ=255)
                      gTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAACGCCCCATGAcccccgg                                                                                                                                 >  2:1230404/1‑146 (MQ=255)
                      gTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACCCCCCATGACCCCTgg                                                                                                                                 >  2:1357750/1‑151 (MQ=255)
                        aCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAATTTAAGCCCCATTAGCGCCgggg                                                                                                                               >  1:810912/1‑149 (MQ=255)
                         cTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGATTAAAGCCACAttaaccctcggcg                                                                                                                              >  1:1427603/1‑139 (MQ=255)
                         cTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATTACCGCTGGCGt                                                                                                                              >  1:746874/1‑151 (MQ=255)
                         cTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGt                                                                                                                              <  1:1154918/151‑1 (MQ=255)
                          ttAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAATTTACGCCCCAATACCGCTggggg                                                                                                                              >  1:1721812/1‑147 (MQ=255)
                            aaTGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAAAGCCACATGACCGCTGGGGGCgg                                                                                                                          >  2:1523875/1‑151 (MQ=255)
                              tGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGac                                                                                                                        <  1:103927/151‑1 (MQ=255)
                               gAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCCGGCGTCGGgca                                                                                                                       >  2:1082703/1‑151 (MQ=255)
                               gAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAGTTAACGCCACATGACCGCTGGccgcgggca                                                                                                                        >  1:1652938/1‑143 (MQ=255)
                                  tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGacac                                                                                                                      <  2:485001/149‑1 (MQ=255)
                                   aaTGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGt                                                                                                                   <  1:830158/151‑1 (MQ=255)
                                      gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCCCATGACCGCCGGCGGCGGACACCGtttt                                                                                                                >  1:501250/1‑150 (MQ=255)
                                      gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGATTTACGCCACATGACCGCTGGCGTCGGACACCGttt                                                                                                                 >  1:854477/1‑150 (MQ=255)
                                      gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGGGTCGGGCACCGttt                                                                                                                 >  2:1216967/1‑150 (MQ=255)
                                       cgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTc                                                                                                                <  1:182821/150‑1 (MQ=255)
                                        gcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATACCGGCTCGGGCTATTGGCGGCGGGGGGATTTTAAAGCACATGGCCCCTGGGGTCCGACACCGTTTCg                                                                                                               >  2:741458/1‑150 (MQ=255)
                                         ccAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGc                                                                                                              <  2:624691/150‑1 (MQ=255)
                                                 cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGCGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGACTTTAAGCCCCCTGACCCCTTGGGTCCGGCAccccttccccgcccttc                                                                                                      >  1:1522389/1‑135 (MQ=255)
                                                 cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTTCCTGCCCTcc                                                                                                      >  2:1031130/1‑150 (MQ=255)
                                                     cAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGTATTTAAGCCACATGAGCGCTGGCGGCTGGACACGTTTCGCCGCCCTCCAGCGa                                                                                                 >  1:860356/1‑151 (MQ=255)
                                                     cAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAAAGCCACATGACCGCCGGCGTCGGACACCCGTTCGCTGCCCTcccacca                                                                                                 >  1:1012645/1‑146 (MQ=255)
                                                          aGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAACGCCCCATGACCGCTGGCGGCGGACACCGTTTCGCCGCCCTCCAGCGACACgg                                                                                           >  2:923560/1‑151 (MQ=255)
                                                            gATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTAAGCCACAAGACCGCCGGCGTCGGAAAACCTTTCGGTGACTTCCACCAACACgggg                                                                                         >  2:330482/1‑151 (MQ=255)
                                                              tctcTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACCTTACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCACCGACACGGGGaa                                                                                       >  2:344912/1‑151 (MQ=255)
                                                                 ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATTACCGCTGGCGTCGGGCACCGGTTTGCTGCCTTCCAGCGGCACGGGGgaaag                                                                                    >  1:1287210/1‑149 (MQ=255)
                                                                 ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATGACCGCTGGGGTCGGACACCGTTTCGCTGCCTTCCCGCGACACGGGGaaaag                                                                                    >  2:800521/1‑149 (MQ=255)
                                                                 ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTAAGCCACATGACCGCTGGCGTCGGACACCGTTTTGCTGCCTTCCCGCGACACGGGGGAACg                                                                                    >  2:1039958/1‑151 (MQ=255)
                                                                 ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCCGCCTTCCAGCGACACGGGGAAACg                                                                                    >  2:896694/1‑151 (MQ=255)
                                                                 ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACg                                                                                    <  1:1353274/151‑1 (MQ=255)
                                                                      gtggtggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGATTAACCCCACATGACCGCTGGGGCCGGGCACCGGTTTGCTGCCTTCCCGGGACACGGGGGAAAGTcccg                                                                               >  2:1396764/1‑149 (MQ=255)
                                                                               aCGCGATATGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGc                                                                      <  1:1290316/150‑1 (MQ=255)
                                                                               aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGc                                                                      <  1:1322289/150‑1 (MQ=255)
                                                                                cgcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCAAATGACCGCTGGCGTCGGACAACGTTTCGCTGCCTTCCCGCGACACGGGGAAAAGGCCGGCCTAAAGCaa                                                                    >  2:565497/1‑151 (MQ=255)
                                                                                    aTACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCg                                                                <  2:299260/151‑1 (MQ=255)
                                                                                            cgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGGGATTTTAACCCACATAACCGCTGGGGTCCGGCACCCCTTTGCTGCCCTCCCCCGACAAGGGGGAAAGGCCCGCGTAACGCGAAATGTTAcccg                                                        >  2:1277362/1‑148 (MQ=255)
                                                                                            cgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGAc                                                       <  2:1214826/151‑1 (MQ=255)
                                                                                             gccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACCCCCCATGACCGCCTGGGTCGGGCACCGTTTTGGCGCCTTTCCGCGACAACGGGGAAAGTCCGGCATAACGGCAAACCGTTAcgac                                                        >  1:22150/1‑147 (MQ=255)
                                                                                              ccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGAcc                                                      <  1:931273/150‑1 (MQ=255)
                                                                                              ccgTCACTGGCGGAGACGGTTAACAAATAGCGGCCCGGTTTATCGGCGGCCGGGAGATTTTACGCCACATTACCGCTGGCGTCGGCCACCCGTTTGCTGCCCTCCAGCGACACGGGGGAACGGCCGGCATAACCGAAATCGTTAcccccca                                                     >  2:1138621/1‑146 (MQ=255)
                                                                                                  cACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAAAGCCCCATGACCGCTGGGGGCGGGAACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAAACGTTACCGACCATTGa                                                 >  1:1116024/1‑151 (MQ=255)
                                                                                                    cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                               <  2:518978/151‑1 (MQ=255)
                                                                                                    cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                               <  1:408894/151‑1 (MQ=255)
                                                                                                    cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                               <  1:919711/151‑1 (MQ=255)
                                                                                                    cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                               <  1:717807/151‑1 (MQ=255)
                                                                                                    cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                               <  2:1472511/151‑1 (MQ=255)
                                                                                                    cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                               <  1:1229810/151‑1 (MQ=255)
                                                                                                     tGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACAAGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATag                                              <  1:653183/151‑2 (MQ=255)
                                                                                                          gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGc                                          <  2:797213/150‑1 (MQ=255)
                                                                                                           agaCGGTTAACAAAGAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCCCATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAAttctg                                         <  1:261733/150‑4 (MQ=255)
                                                                                                             aCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTg                                        <  2:1246527/149‑1 (MQ=255)
                                                                                                              cGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCGACCATTTATAATTGCTGAg                                      >  2:886556/1‑150 (MQ=255)
                                                                                                              cGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGGCGGACACCGTTTTGCTGCCTTCCAGCGACACGGGGAAAAGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTTCTGAg                                      >  2:975579/1‑150 (MQ=255)
                                                                                                                gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTTGCGGCGGACACCGCTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCCGCATAACGCAAAGCGGTACCCGACCTTGATTATTTCTTAgcg                                    >  1:723203/1‑150 (MQ=255)
                                                                                                                gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTTCTGAgcg                                    >  2:553907/1‑150 (MQ=255)
                                                                                                                gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTTATAATTGCTGAgcg                                    >  2:152124/1‑150 (MQ=255)
                                                                                                                gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcg                                    <  1:1661982/150‑1 (MQ=255)
                                                                                                                gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcg                                    >  1:700267/1‑150 (MQ=255)
                                                                                                                gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGAATTAACGGCACATTACCGCCCGGGGGGGAAACCCGTTTGGTGGCTTCCAGAGACACGGGGAAAAGGACGGGATCACCCAAAACGTTTACGACCATTGATAAATGCCGAggcg                                   >  1:489337/1‑148 (MQ=255)
                                                                                                                     aCAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCATGACCGCTGGCGTCGGACACCGGTTTGCTGCCTTCCCGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCGACCATTTATAATTTCTGAGCGCGCaa                               >  1:1529720/1‑151 (MQ=255)
                                                                                                                     aCAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGGGTCGGACACCGTTTCGCTGCCTTCCCGCGCACCGGGGGAACGGCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGGGCGGGCaa                               >  2:1195501/1‑151 (MQ=255)
                                                                                                                     aCAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACAAGACCGCTGGCGTCGGGCACCGTTTCGCGGCCTTCCAGCGACACGGGGGAACGTCCGGCCTAACCCAAAACGTTACCGACCATTGATAATTGCTGAGCGCGCaa                               >  1:1717176/1‑151 (MQ=255)
                                                                                                                       aaaTAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCCGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGCCACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTTCTGAGCGCGCaaa                              >  2:1479937/1‑150 (MQ=255)
                                                                                                                        aaTAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACt                            <  2:1634967/151‑1 (MQ=255)
                                                                                                                          tAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTTACGCCCCATGACCGCTGGCGTCCGACACCGTTTCGCTGCCTTTCAGCGACACGGGGAAACGGCCGGCATAACGCAAATCGTTTCCGACCATTGATAATTGCTGAGCGCGCAAACTTa                          >  1:994178/1‑151 (MQ=255)
                                                                                                                               gCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                      <  1:395799/150‑1 (MQ=255)
                                                                                                                               gCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                      <  2:1432909/150‑1 (MQ=255)
                                                                                                                                cTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGc                     <  1:1433309/150‑1 (MQ=255)
                                                                                                                                  cGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcg                  <  1:663431/151‑1 (MQ=255)
                                                                                                                                   ggTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGGTTCGCTGCCTTCCCGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTGCTGAGCGCGCAAACTTAACTGCAcg                  >  1:439756/1‑150 (MQ=255)
                                                                                                                                    gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg                >  1:944410/1‑151 (MQ=255)
                                                                                                                                    gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg                >  1:1236110/1‑151 (MQ=255)
                                                                                                                                            gcggcggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCAttt         <  2:1011212/150‑1 (MQ=255)
                                                                                                                                                  gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc   <  2:437357/150‑1 (MQ=255)
                                                                                                                                                  gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc   <  1:1039411/150‑1 (MQ=255)
                                                                                                                                                  gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc   <  2:39540/150‑1 (MQ=255)
                                                                                                                                                  gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc   <  1:486304/150‑1 (MQ=255)
                                                                                                                                                  gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCg  >  1:870615/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
GCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCG  >  NC_000913/2332675‑2332970

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