Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,1760TC44.4% 113.8 / 60.3 100E304G (GAA→GGA) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (43/12);  new base C (0/44);  total (43/57)
Fisher's exact test for biased strand distribution p-value = 3.74e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGG  >  NC_000913/3601026‑3601326
                                                                                                                                                         |                                                                                                                                                       
tattatTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACtt                                                                                                                                                         >  2:477477/1‑151 (MQ=255)
 attattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTc                                                                                                                                                        >  2:1151164/1‑151 (MQ=255)
 attattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTc                                                                                                                                                        >  2:1402100/1‑151 (MQ=255)
   tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTcc                                                                                                                                                       <  1:1368048/150‑1 (MQ=255)
    attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCg                                                                                                                                                      <  1:1921293/150‑1 (MQ=255)
       cAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAAc                                                                                                                                                 <  2:1852684/151‑1 (MQ=255)
          gCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAAccgc                                                                                                                                              <  2:10653/151‑1 (MQ=255)
           cGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAAccgcc                                                                                                                                             <  1:686185/151‑1 (MQ=255)
           cGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAAccgcc                                                                                                                                             <  2:1885666/151‑1 (MQ=255)
             tAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAAccgccg                                                                                                                                            >  2:57595/1‑150 (MQ=255)
              aaaGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGAt                                                                                                                                          >  2:1458830/1‑151 (MQ=255)
              aaaGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGAt                                                                                                                                          >  2:577431/1‑151 (MQ=255)
                      gTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCg                                                                                                                                  >  2:1454059/1‑151 (MQ=255)
                        aaaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCgc                                                                                                                                 >  2:345533/1‑150 (MQ=255)
                           tgttgtACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCgcagc                                                                                                                              >  1:233480/1‑150 (MQ=255)
                               gccccGTGGCAAGCCCCGCCCCCCCCCAGCCAATCATCGACCATACCAGCAGAAATAGCAGGCATTCATCGAATAATAAGCCATTGATAACAACCTTCGTCCAAAAAACCGGTGCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCa                                                                                                                          <  1:642689/147‑1 (MQ=255)
                                     tGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTc                                                                                                                   >  2:1555015/1‑151 (MQ=255)
                                     ggccagcccACGCCGCCACCCACCCAATCACGGAAAATACCAGAAAAAAGGGCGGGCTTCGATAAAATAATAACCCTTTGATATAAAGCTACGTCCAAAAACCCGGTCCGACTACCGAAACCGCCGATGACAATCGCAGCACCAGAATTTc                                                                                                                   <  2:20897/144‑1 (MQ=255)
                                         aacccccgcccccccccAGCCAATCACCGAGCATCCCAGCAAAAATAGCAGGCTTTCATCAAATGATAACCCATTGATATAAACCTTCGTCCCAAAACCCTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGa                                                                                                                <  2:221493/136‑1 (MQ=255)
                                          acccccGCTCCCCCCCACCCAATCATCGAGCATACCAGCAGAAATAGCAGGCATCCATCGAATGATAAGCCATTGATACAAAACTTCTTCCCAAAACCCTGTCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGa                                                                                                                <  2:1930774/144‑1 (MQ=255)
                                          aGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGAc                                                                                                               >  2:1183162/1‑150 (MQ=255)
                                            ccccgccgcCCCCAACCCAATCACCGAGAATACCAGCAGAAAAAGCAGACATTCAGAAAATGAAAAGCCATTAATATCAACCTTCGTCCCAAAACCCGGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGaaa                                                                                                              <  2:1766854/146‑3 (MQ=255)
                                            cccccctccccccaaGCCAATCACCGACCAACCAAGCACAAAAAGCAGGCATTCATCGAAAGAAAAGCCATTAAAAAAAAACTTCTTCCCAAAAACTTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaa                                                                                                            <  1:140839/137‑1 (MQ=255)
                                            ccACGCTGCCCCCCAGCCAATCATCGAGCATACCACCAAAAATAGCAGGCATTCATCGAATGATAAGCCATTAATATCAACCTTCGTCCCGAAACCCTGTGCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa                                                                                                             <  2:355841/150‑1 (MQ=255)
                                              aCGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt                                                                                                          >  1:1205604/1‑151 (MQ=255)
                                               cGCTCCCCCCCACCAAACCATAGAGCACACAAGCGAAAAAACCGGTCATTCAGAGAATAAAAGCCAATTAATACAAACCTTCCTTCCAAAACCGTGGCCCATTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACa                                                                                                              <  2:341893/146‑1 (MQ=255)
                                                 cccccccccACCAAAACATCGAGCACCCCAACAAAAATGGCAGGCATTCATGGAATAATAACCCTTTGATAACAAACTCCTTCCAAAAACCCTTTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTa                                                                                                        <  1:1589373/144‑1 (MQ=255)
                                                  tGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt                                                                                                      >  2:118532/1‑151 (MQ=255)
                                                    ccccccacccaATCACCGAACACCCCAGCAAAAATAGCAGGCTTTCATAAAATAATACCCCATTAATACAAACCTTCTTCCCAAAACCCTTTCCCCCCCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaa                                                                                                            <  1:448132/140‑1 (MQ=255)
                                                    ccccccacccaACCACCGAGCATACCAGAACAAAAAGCAGGCATCCATCAAATAATAACCCATTGATACAAACCTTCTTCCCAAAACCCTGTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGa                                                                                                    <  1:540824/148‑1 (MQ=255)
                                                    ccccccAGCCAATCATCGACCATCCCAGCAGAAATAGCAGGCTTTAATCGAATAATAACCCATTAATAACAACCTTCGTCCCAAAACCCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTaa                                                                                                       <  2:1152676/145‑1 (MQ=255)
                                                     cACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac                                                                                                   >  2:877079/1‑151 (MQ=255)
                                                     cACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac                                                                                                   >  2:910225/1‑151 (MQ=255)
                                                      cccaaGCAAATCATCGAACACACCAGCAGAAAAAGCAGTCATTCATCAAATAAAAACCAATTAATATAAAACTTCGTCCCAAAACCTTGTCCCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac                                                                                                   <  1:1928737/146‑1 (MQ=255)
                                                        ccAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaac                                                                                                >  2:120777/1‑151 (MQ=255)
                                                            ccAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                                            >  1:801877/1‑151 (MQ=255)
                                                            ccAATCATCGAACACCCCAACAACAATAGCAGGTTTTCATAGAATGATAACCTTTTAATATAAATCTTCTTCCAAAAACCCTGTCCCCCTTCCGAATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                                            <  2:1150828/151‑1 (MQ=255)
                                                              aaTCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcg                                                                                          >  1:1537564/1‑151 (MQ=255)
                                                                  acggacaaTACCAGCAAAAATAGCAGCTTTTCATCGAATAATAACCTATTAACACAAACCTCCGTCCCAAAACCCTGCCCCCCCCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCa                                                                                       <  2:448194/143‑1 (MQ=255)
                                                                  accgaaaaTACCACCAGAAATCGCGGCCTTTCATCAAATAAAAAGCATTTGATATAAACCTCCGTCCCAAAACCCTGTGCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAAt                                                                                             <  2:1020427/137‑1 (MQ=255)
                                                                   tcGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGcc                                                                                        >  2:934523/1‑148 (MQ=255)
                                                                   tcGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGcc                                                                                        >  2:6414/1‑148 (MQ=255)
                                                                    cGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGa                                                                                     >  1:394240/1‑150 (MQ=255)
                                                                    cGACCACACCAACAACAATACAAGACATTCATCAAAAAAAAAGCCATTGATACAAAACTTCTTCCCAAAAACCTTTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg                                                                                      <  1:207899/149‑1 (MQ=255)
                                                                     gAGCAACCCAGAAGAAAAAGCAGGAATTAATAAAATGATACGCTTTTTATATCAACCTTCGTCCCAAAACCCTTTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg                                                                                      <  2:1936372/148‑1 (MQ=255)
                                                                      agcaacccaccaGAAAGAGCAGACATTCATCAAATGAAAACCTATTAATAAAATACTTCTTTCAAAAACCCTGCCCCCCTCACCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg                                                                                      <  1:1449499/141‑1 (MQ=255)
                                                                         aTACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATcc                                                                               >  1:950424/1‑151 (MQ=255)
                                                                             cagcacAAATGGCAGCCATTCATAAAATAAAAACCAATTAATATCAACCTTCTTCCCAAAACCCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCca                                                                             <  2:1513722/149‑1 (MQ=255)
                                                                                 agcaATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTg                                                                       >  2:1900742/1‑151 (MQ=255)
                                                                                 agaaaaaGCAGGCATTCATCAAATGAAAACCCATTAATATAAACCTTCGTCCCAAAACCCTGTCCCATTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATcc                                                                               <  2:938369/137‑1 (MQ=255)
                                                                                  gcaATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                                      >  2:814633/1‑151 (MQ=255)
                                                                                   caATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgcg                                                                     >  1:1671945/1‑151 (MQ=255)
                                                                                    aagaGCAGGCATTCATCAAATGATAACCCATTAATATAAACCTTCGTCCAAAAACCCTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                                      <  1:312141/146‑1 (MQ=255)
                                                                                    aaaagcagcCATTCTCCAAATAATAAGCATTTTATAACAACCTTCTTCCCAAACCCCTGGCCACTTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgcgc                                                                    <  2:1747599/148‑1 (MQ=255)
                                                                                    aaTAGCAGGCTTTCACCAAATAAAAAGCTATTAATATAAACCTTCGTCCCAAAAACCTGTCCCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcag                                                                            <  2:710312/143‑1 (MQ=255)
                                                                                    aaTAGCAGGCATTCATCAAATGAAAAGCATTTAATATCAACCTTCTTCCCAAAACCCTGTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                                      <  2:1809623/149‑1 (MQ=255)
                                                                                      tAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcc                                                                   >  1:1862959/1‑150 (MQ=255)
                                                                                        gCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccacc                                                                >  2:419541/1‑151 (MQ=255)
                                                                                          tggCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt                                                              >  1:461194/2‑151 (MQ=255)
                                                                                           ggaaTTCATCAAATAATAAGCTATTAATACAAACCTTCGTCCCAAAACCCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcca                                                                  <  1:1614212/143‑1 (MQ=255)
                                                                                           ggCATTCATCAAATAAAAAGCCTTTAATACAAAACTTCTTCCCAAAACCCTGGCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATa                                                             <  2:1192079/151‑1 (MQ=255)
                                                                                             cATTCATCAATTGATAACCCTTTGATATAAAACTTCGTTCCAAAACCCTGTGCCTCTTCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACa                                                           <  2:1434752/151‑1 (MQ=255)
                                                                                                      aaTGATAAGCTATTGATAACAAACTTCTTCCCGAAACCCTGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccacca                                                               <  1:5860/138‑1 (MQ=255)
                                                                                                       aTGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcg                                                 >  2:714018/1‑151 (MQ=255)
                                                                                                           aaaaccaaTTAATATAAACCTTCGTCCCAAAACCCTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc                                                <  1:534852/142‑1 (MQ=255)
                                                                                                            aaGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg                                             >  1:417795/1‑150 (MQ=255)
                                                                                                             acccctttaatacCAACCTTCTTCCAAAAACCTTGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg                                             <  2:1475155/141‑1 (MQ=255)
                                                                                                             aatcatttaatatCAATCTTCTTCCAAAAACCCTGTCCCACTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc                                                <  2:453887/138‑1 (MQ=255)
                                                                                                             aGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                           >  2:605404/1‑151 (MQ=255)
                                                                                                              cccATTAATAAAAACCTTCGTCCAAAAACCTTGCGCCACTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgc                                                  <  1:1736807/142‑1 (MQ=255)
                                                                                                                  aTTGATATAAACCTCCGTCCCGAAACCCTGTCCCCCTTCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGa                                         <  1:451702/149‑1 (MQ=255)
                                                                                                                   ttGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                                      >  1:109685/1‑151 (MQ=255)
                                                                                                                    taatacCAACCTTCGTCCAAAAACCCTGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                                      <  2:128214/148‑1 (MQ=255)
                                                                                                                    taataaaaatctccGTCCCAAAACCCTGTCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGa                                         <  1:1137650/134‑1 (MQ=255)
                                                                                                                      atatCAAACTTCTTCCCAAAACCCTGTGCCTCTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATcc                                     <  2:759398/149‑1 (MQ=255)
                                                                                                                       tatCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATc                                  >  2:705427/1‑151 (MQ=255)
                                                                                                                         tCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTg                                >  2:776284/1‑151 (MQ=255)
                                                                                                                           aatcTTCGTCCCGAAACCCTGTCCCCCTCCCGCAGCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTg                                <  2:1422313/146‑1 (MQ=255)
                                                                                                                           aaccTTCGTCCAAAAACCCTGCCCCCCTCCCCCACCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTc                              <  1:1581872/148‑1 (MQ=255)
                                                                                                                           aaaCTTCTTCCAAAAACCCTGTCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCt                                 <  1:875601/148‑1 (MQ=255)
                                                                                                                            aaCTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCg                             >  1:366600/1‑151 (MQ=255)
                                                                                                                            aaCCTCCTTCCCAAAACCCTGCGCCACTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                                      <  1:983739/143‑1 (MQ=255)
                                                                                                                             aCTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGc                            >  2:1866197/1‑151 (MQ=255)
                                                                                                                              cTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACAATACAGGAACGGGCGCAGGATGAATCCATCTGTCGc                            >  1:1074592/1‑150 (MQ=255)
                                                                                                                                  gtcccaaaacCCTGTCCCATTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCa                        <  2:1097349/144‑1 (MQ=255)
                                                                                                                                   ttCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCaa                       >  2:1032967/1‑150 (MQ=255)
                                                                                                                                   tccaaaaaaCCGGCCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCAc                          <  1:1653155/142‑1 (MQ=255)
                                                                                                                                       gAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt                  >  2:155608/1‑151 (MQ=255)
                                                                                                                                       aaaacCCTGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt                  <  1:1383919/150‑1 (MQ=255)
                                                                                                                                        aaacCCTGTCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAg                   <  1:1300426/149‑1 (MQ=255)
                                                                                                                                        aaacCCGGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc                     <  2:767103/147‑1 (MQ=255)
                                                                                                                                        aaaaCCCGGCCCCCTTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCa                    <  2:1083904/149‑1 (MQ=255)
                                                                                                                                         aaaCCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt                >  1:1063908/1‑151 (MQ=255)
                                                                                                                                           aCCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCt              >  1:1844316/1‑151 (MQ=255)
                                                                                                                                                tgtgCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt            >  1:1499773/1‑150 (MQ=255)
                                                                                                                                                  tgCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACg         >  2:572143/1‑151 (MQ=255)
                                                                                                                                                    ccACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGtt       >  2:1771002/1‑151 (MQ=255)
                                                                                                                                                    ccACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTa           <  2:642661/147‑1 (MQ=255)
                                                                                                                                                        ttCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATg   >  1:1937176/1‑151 (MQ=255)
                                                                                                                                                         tCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATgg  >  2:1621752/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAACCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATGG  >  NC_000913/3601026‑3601326

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: