Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,7270GA41.9% 100.4 / 45.7 74intergenic (‑108/+33)uspF/ompNstress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (38/5);  new base A (1/30);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.91e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGG  >  NC_000913/1435580‑1435872
                                                                                                                                                   |                                                                                                                                                 
aaTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGaa                                                                                                                                                 >  1:118985/1‑150 (MQ=255)
   cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAg                                                                                                                                             <  1:888999/151‑1 (MQ=255)
     gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGcc                                                                                                                                           >  2:1275882/1‑151 (MQ=255)
     gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGcc                                                                                                                                           <  1:202568/151‑1 (MQ=255)
      aaaTATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGCCGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGccc                                                                                                                                          <  1:170408/151‑1 (MQ=255)
        atatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGc                                                                                                                                        <  1:202260/151‑1 (MQ=255)
         tatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                       >  2:294463/1‑151 (MQ=255)
               tAGGGCCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTATGATCCGCATAAAGTGTAACGCTAGATGACACTTTGTGTGTAGTGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAGGAGTATGCAGAAAAAGCCCGCTATaa                                                                                                                                   >  2:614103/1‑149 (MQ=255)
                aGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAAc                                                                                                                                  >  1:531314/1‑149 (MQ=255)
                    aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                                                            >  2:60917/1‑151 (MQ=255)
                     cAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                           >  1:407595/1‑151 (MQ=255)
                     cAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                           >  1:434549/1‑151 (MQ=255)
                      aaGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                           >  1:952488/1‑150 (MQ=255)
                       aaattggTTCTTTTCAAAAACCCTCTTGTTTTATTAGCCGAAAAAATTGAACCGCAAGATAACTTTTTTTGTGAAAGGACGGATTCAACTTTTTTAATTAAGAAAAAAGGAGGAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCag                                                                                                                         <  2:451390/145‑1 (MQ=255)
                             gTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                   >  1:853840/1‑151 (MQ=255)
                                cTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                                                                 >  1:1337877/1‑150 (MQ=255)
                                 ttttaaaaaaaCCTCCTTTTTTAGTATCCGAAAAAAGTTTAACGCCAGATGACTTTTTTTGTGTAAGGACGGAGTTCATATTTTTAATTTAGAAAAAAGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTa                                                                                                                <  1:1030903/143‑1 (MQ=255)
                                  gTTCATAAACCCTCTTTTTTTAGTATCCGAAAAAAGTGTAACCCCAGATGACATTTTTTGTGAAATGAGGGATTTCACTTTTTTAATTTAGATAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                                                               <  2:488972/150‑1 (MQ=255)
                                      aaaaacCCTCTTGTTTAAGTACCCGAATAAAGTGAAACGCAAAAGGACATTTTTTGTGAAATGACGGATTTCAATTTTTTAATTTAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                          <  2:823420/149‑1 (MQ=255)
                                      aaaaacCCTCCTTTTTTATTATCCGAAAAAATTGTAACGCCAGATGACTCTTTTTGTGAAATGAGGGATTTCACTTTTTTATTTAAAATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                                                            <  2:1329759/147‑1 (MQ=255)
                                        aaacCCTCCTTTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACTTTTTTTGTGTAATGACGGAGTTCACTTTTTTATTTTAGAAAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                         <  1:299529/150‑1 (MQ=255)
                                         aaccctcctGTTTTAGTATCCGAAAAAAGTGTACCGCCAGATGACTTTTTTTGTGAAATGAGGGAGTTAATTTTTTTAATTTAAACAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa                                                                                                       <  1:1097143/148‑1 (MQ=255)
                                         aaaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                        >  1:392640/1‑150 (MQ=255)
                                            cctcttTTTTTAGTATCCGAAAAAAGTGTACCGCCAGATGACTCTTTTTGGGAAATGACGGATTTCACTTTTTAAATTTAGAACAAAGGAGGAAAAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa                                                                                                      <  2:1248430/149‑1 (MQ=255)
                                                   ttttAGTATCCACAAAAAGTGTAACTCCAGATGACACTTTTTGAGTAATGACGGAATCCACATTTTTAATTAAAATAAAGGGAGAAAAAAAAATCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCa                                                                                                  <  2:1233427/146‑1 (MQ=255)
                                                    tttAGTTCCCGAAAAAAGTTAAACGCGAGATGACATTTTTTGTGAAATGACGGATTTCACTTTTTAAATTAAGAAAAAGGGAGGAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                                                               <  2:331082/148‑1 (MQ=255)
                                                     ttAGTACCCGAATAAGGTGTACCGCCAGATGACATTTTTGGTGAAATGAGGGATTTCACATTTTTATTTTAAACCAAAGGAGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTa                                                                                            <  1:44790/150‑1 (MQ=255)
                                                      tAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa                                                                                           >  2:1324342/1‑150 (MQ=255)
                                                            ccGCATAAAGTGTAACGCCAGATGACTTTTTTTGTGTAATGAGGGATTTCACTTTTTTAATTTAGATAAAAGGAGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGa                                                                                      <  1:245657/149‑1 (MQ=255)
                                                             cataaaaggtgtAACGCCAATTGACAGTTTTTGAATAAGGAGGGTGTGCACATTTTAATTAAAAATGAAAGGGAAAATAAAAACGAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAg                                                                                         <  2:390286/137‑1 (MQ=255)
                                                              gCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                  >  1:939906/1‑151 (MQ=255)
                                                                   aaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcat                                                                             >  1:1000972/1‑151 (MQ=255)
                                                                         gatgCCGGATGACACGTTTTGGGGAAAGAGGGAGTTCACATTGTTAATTTTGATCGGAGGAGGAAGAATAAAAAAAAAAAACGCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAATGCGACAAGATCATCGGt                                                                         >  2:605547/4‑149 (MQ=255)
                                                                          aCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                         >  1:462987/1‑148 (MQ=255)
                                                                           cGCCAGATGACTTTTTTTGGGAATGGAGGGATTTCAATTTTTTAATTAAAATAAAGGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                                                        <  1:294305/148‑1 (MQ=255)
                                                                            gCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                       >  2:61684/1‑148 (MQ=255)
                                                                             ccAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                       >  1:1136098/1‑147 (MQ=255)
                                                                                gatgcctttttttGTGAAAGGAGGAAGTTCACATTTTTAATTTAGATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg                                                                  <  2:1155510/141‑1 (MQ=255)
                                                                                    acacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                                                            >  2:748091/1‑151 (MQ=255)
                                                                                     cacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                           >  2:267738/1‑151 (MQ=255)
                                                                                      acTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                          >  2:348081/1‑151 (MQ=255)
                                                                                        tttttGTGAAATGACGGATTTCACTTTTTAAATTAAGAAAAAAGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                           <  1:1181471/148‑1 (MQ=255)
                                                                                        tttttGAGTAATGAGGTGGTTCGCATTTTAAATGTAGATCAAAGAAGAAAAAATAAGAAGAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                                                           <  1:1247543/148‑1 (MQ=255)
                                                                                            tgtgTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                                     >  2:52441/1‑150 (MQ=255)
                                                                                                tAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                  >  1:1094259/1‑149 (MQ=255)
                                                                                                 aaggagggATTTCACATTTTTAATTTAAAAAAAAGGAGAAAGAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa                                                    <  1:349470/140‑1 (MQ=255)
                                                                                                 aaTGCCGGAGTTCAATTTTTTAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa                                                    <  2:1313759/146‑1 (MQ=255)
                                                                                                 aaTGACGGAGTTCACATTTTTAATTAAAAAAAAAGGAGGAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAg                                                   <  2:1178312/147‑1 (MQ=255)
                                                                                                       ggAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCAt                                         >  2:1005996/1‑151 (MQ=255)
                                                                                                       ggAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCAt                                         >  2:81098/1‑151 (MQ=255)
                                                                                                         aGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                       >  1:730778/1‑151 (MQ=255)
                                                                                                            tCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc                                     >  2:566156/1‑150 (MQ=255)
                                                                                                             cacaTTTTAAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCa                                    <  1:16366/150‑1 (MQ=255)
                                                                                                              acttttttAATTTAGAAAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                   <  1:250900/147‑1 (MQ=255)
                                                                                                                 tttttAATTAAAAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                                <  2:344411/150‑1 (MQ=255)
                                                                                                                      aaTTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                          >  2:1067525/1‑151 (MQ=255)
                                                                                                                      aaTTAAGACAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                           <  1:188932/150‑1 (MQ=255)
                                                                                                                        ttaaaaaaaaaGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                          <  1:524911/141‑1 (MQ=255)
                                                                                                                        taaaaaacAAAGGAGAAAAAATAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACa                             <  2:1001496/139‑1 (MQ=255)
                                                                                                                         ttAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                        >  1:1127828/1‑150 (MQ=255)
                                                                                                                           aGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt                     >  2:1158754/1‑151 (MQ=255)
                                                                                                                             aTCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                   >  2:364382/1‑151 (MQ=255)
                                                                                                                              tCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                   >  1:1101883/1‑150 (MQ=255)
                                                                                                                              tCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                   >  2:422393/1‑150 (MQ=255)
                                                                                                                                 aaGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                 >  1:405163/1‑149 (MQ=255)
                                                                                                                                   ggaggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc             >  2:89826/1‑151 (MQ=255)
                                                                                                                                   ggagaaAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt              <  1:129588/150‑1 (MQ=255)
                                                                                                                                   ggagaaAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc             <  2:463184/151‑1 (MQ=255)
                                                                                                                                   ggagaaAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc             <  1:60878/151‑1 (MQ=255)
                                                                                                                                    gaggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGATGATTTTATAGTc             >  2:861111/1‑150 (MQ=255)
                                                                                                                                        aaaaaTAAACAAAAAAACCCCGCCATAACAGCGGGCAGGAGGACTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACg         <  2:763780/147‑1 (MQ=255)
                                                                                                                                              aaGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTgg  >  1:667139/1‑151 (MQ=255)
                                                                                                                                              aaGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAGGTAAGAGTCaaggg        >  1:574111/1‑142 (MQ=255)
                                                                                                                                              aaGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTgg  <  1:28471/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGG  >  NC_000913/1435580‑1435872

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