Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,5960GA41.5% 66.3 / 40.7 66F13F (TTC→TTTymgDuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (32/6);  new base A (0/27);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 9.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAC  >  NC_000913/1222452‑1222739
                                                                                                                                                      |                                                                                                                                                
gCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaac                                                                                                                                             >  2:811706/1‑150 (MQ=255)
 cTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaaca                                                                                                                                            <  2:135732/150‑1 (MQ=255)
      attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAgg                                                                                                                                      <  2:302940/151‑1 (MQ=255)
       ttacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGc                                                                                                                                    >  1:688542/1‑151 (MQ=255)
        tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGcc                                                                                                                                   >  2:726610/1‑151 (MQ=255)
              ttaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGcc                                                                                                                               >  2:447688/1‑149 (MQ=255)
               taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGt                                                                                                                            >  2:91265/1‑151 (MQ=255)
                 ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAg                                                                                                                             <  1:827488/148‑1 (MQ=255)
                     aacatTTTCTGCACCCTGTGCGCTTTCATTTTTTGCCTTTTCGTTGGTTTCTGCGAAAGAAATAAAGGTTTCGGACGGCAACCCTTTTTAAATATTGTAGTCCGGCGCCCAGGTCTTGGGTGCAAAAACAAACAGGCCCGCCAGTAAAg                                                                                                                        <  1:405126/149‑1 (MQ=255)
                     aacaaTTTCTGCCACCTGTGCGTTGTGATTTTTTGCCTTTTCGTTGGCTTCTGCGCAAGTAATATGGTTTTCTGACGGCAAACCTTTTTTAATATTGTAGCCCGGCCCCCGGGTCATGGGGGCAAAAAAAAACAGGCCCGCCAGTAAAg                                                                                                                        <  2:200386/149‑1 (MQ=255)
                        aaTTTCCGCAACCTGTGCGCTATCATTTTTTGCCTTTTGGTGGACTTCGGCGCAATAAATATAGGTTTCTGACGGCAACCTTTTTTAAATATTGTCGTCCCGCGCCCAGGCCATGGGTGCGAAAAAAAACAGCCCCGCCAGTAAAg                                                                                                                        <  2:630509/146‑1 (MQ=255)
                                 accttgtgCGCTGTCAGTTTTTGCCTTTTCGTTGGTTTCTGCGCAAGAAATATGGGTTTCGGACGCCAACCCTTTTTAAATATTGTAGTCCTCCGCCAAGTTAATTGGCGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtca                                                                                                            <  1:807937/145‑1 (MQ=255)
                                 acccTGGGCGCTTTCTGTTTTTGCCTTTTCGTTGGTTTCTGCGAAAGTAATATAGGTTTCGGACGGCAACCTTTTTTAAATTTTGTAGCCCGGCCCCAAGGTTTTTGGGGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAAttttt                                                                                                               <  2:572791/144‑1 (MQ=255)
                                 aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                          >  2:509519/1‑151 (MQ=255)
                                  aacTGTGCGCTGACAGTTTTTGCCATTTCGTTGGCTTCTGAGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTAATTGGAGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                          >  2:537250/3‑150 (MQ=255)
                                   cttgtgggcgctCATTTTTTGACATTCGGTGGGATTCCGCTAAAATAATATATTTTTCTGACGAATAGCCTTTTTTAATGTGGTAGCTTGCCCCCCAATTCTTCGGTTCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtc                                                                                                             <  2:913417/135‑1 (MQ=255)
                                   ccTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcat                                                                                                        >  2:918688/1‑151 (MQ=255)
                                          gcTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                 >  2:142986/1‑151 (MQ=255)
                                           cTGTCATTTTTTGCCTTTTGGTTGGCTTCTGCGAAAGAAATATGGTTTTCTGACGACAACCCTTTTTTAATTTTGTCGCCCGCCCCCCAGTTTTTTGGGGCAAAAAAAACCGGCCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                     <  2:501467/146‑1 (MQ=255)
                                              tCATTTTTTGCCTTTTCGTTGGCTTCTGCGAAAGTAATATAGGTTTCGGACGGCAACCTGTTTTAAATTTTGTAGCCCTGCCCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                              <  1:589924/150‑1 (MQ=255)
                                                 ttttttGCCTTTTCTTTGTTTTCGGCGCAAATAATATAGTTTTCGGACGCCAACCCTTTTTTAATATTGTAGTCCTGCCCCCAGTTTTTTGGAGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                               <  2:651593/145‑1 (MQ=255)
                                                  tttttCCTATTGGGTTGGTTTTTGCGCAAAAAATATAGGTTGCTGCGGAAAACCCGTTTTAATTTTTGTAGTCCTGCGCCCGGGTGGTGGGTCCGAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                          <  1:195498/150‑1 (MQ=255)
                                                  ttttggctttttCGTTGTCTTCGGCACAATAATTATTGTTTTGTGAGGGCAAACTTTTTTTAATATTGCTGCCCGACGCCCAGTTATTTGCTACAAAACAAAACTGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                           <  2:395247/140‑1 (MQ=255)
                                                  ttttaGCCATTTCTTTGGTTTGTGCCCAAATAAGATAGGTTGCTGACGGGAAACTGTTTTTAATTTTGAAGCCCCGCGCCCAGTTTATGGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                           <  2:712602/149‑1 (MQ=255)
                                                   ttttGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                        >  1:368671/1‑151 (MQ=255)
                                                    tttGCCTTTTGGTTGCCTTCTGCGCAAGTAATATAGTTTTCTGACGAAAACCCGTTTTAAATATTGTAGCCCGGCCCCCAGTTCATTGGTGAAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                        <  2:456264/150‑1 (MQ=255)
                                                          aTTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                  >  1:845590/1‑150 (MQ=255)
                                                          aTTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                  >  1:1025845/1‑150 (MQ=255)
                                                          aTTTCGTTGGCTTCTGCGCAAGTAATATAGGGTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                  >  2:740081/1‑150 (MQ=255)
                                                           tttCGTTGGTTTCTGCACAAGTAATATAGTTTTCTGACGAAAACCCTTTTTAAATATTGTAGCCCTGCCCCCAGTTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                  <  1:421330/149‑1 (MQ=255)
                                                           tttCGTTGGCTTCTGCGAAAGTAATATAGTTTTCTGACGGCAAACTTTTTTTAATTTTGTAGTCCTGCGCCCAGTTTTTTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTa                                                                                     <  2:705624/146‑1 (MQ=255)
                                                                ttGGTTTCGGCGAAAGAAATATAGGTTTCTGACGGAAACCCGTTTTAAATATTGTAGCCCGGCGCCAAGGTTATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                            <  2:83428/150‑1 (MQ=255)
                                                                ttGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                            >  1:41743/1‑150 (MQ=255)
                                                                 tGGTTTCTGCAAAAGAAATAGGGTTTTGTGGGGGAAAACCTTTTTTAATTTGTTAGCCCGGCCCCAGGTTTTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                           <  2:482999/150‑1 (MQ=255)
                                                                  ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCTTTTTTTAATATTGTAGTCCGGCCCCCAGGTTTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                           <  1:1000368/149‑1 (MQ=255)
                                                                  ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                          >  1:1118145/1‑150 (MQ=255)
                                                                    cTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                       >  2:16819/1‑151 (MQ=255)
                                                                     ttCTGCGCAATTAATATGGTTTTCTGACGGCAACCCTTTTTTAATTTTGTAGCCCTCCCCCCAGTTTTTTGGTGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGca                                                                         <  2:1054075/148‑1 (MQ=255)
                                                                     ttCTGCGAAAGTAATATAGTTTTCTGCCGGAAACCTTTTTTTAATTTTGTAGCCCTGCCCCCAGGTCTTTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACACAGa                                                                               <  1:794442/142‑1 (MQ=255)
                                                                             aaGAAATATAGGTTCCTGACGCAAACCCTTTTTAAATATTGTAGCCCGGCCCCGGGTTCTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAcc                                                              <  1:191737/151‑1 (MQ=255)
                                                                                             tGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCt                                              >  2:36777/1‑151 (MQ=255)
                                                                                             ggACGGAAACCCTTTTTTAATATTGTAGCCCGGCCCCCAGGTCTTGGGGGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCt                                              <  1:310776/150‑1 (MQ=255)
                                                                                              gACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                             >  1:975191/1‑151 (MQ=255)
                                                                                               aCGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                            >  1:633663/1‑151 (MQ=255)
                                                                                                    aaaCCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCg                                       >  2:398458/1‑151 (MQ=255)
                                                                                                     aaCCTTTTTTAAATTTTTTGTCCCTCCCCCCAGTTTTTTGGTGAAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                        <  1:822691/150‑1 (MQ=255)
                                                                                                     aaCCGTTTTTAAATATTGTAGCCCTGCCCCGAGTAGATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                        <  2:900086/150‑1 (MQ=255)
                                                                                                     aaCCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGc                                      >  1:264360/1‑151 (MQ=255)
                                                                                                        ccGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                    >  1:347149/1‑151 (MQ=255)
                                                                                                           tttttAATGTTGTAGCCCGGCGCCAAGGTCTTGGGGGAAAAAAAAACCAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                    <  1:1212015/148‑1 (MQ=255)
                                                                                                           tttttAATATTGTAGCCCGCCCCCCAGTTCTTTGGTGCAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcccc                                   <  1:1054787/149‑1 (MQ=255)
                                                                                                             tttAATTTTGTAGTCCTGCGCCCAGGTCATTGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                <  1:480567/150‑1 (MQ=255)
                                                                                                                   tatTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                          >  1:1034733/1‑151 (MQ=255)
                                                                                                                     ttgttgTCCCGCCCCCAGTTATTGGTTCCAAAAACAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                         <  2:373786/150‑1 (MQ=255)
                                                                                                                       gTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTg                      >  2:907858/1‑151 (MQ=255)
                                                                                                                        tAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGc                     >  1:215391/1‑151 (MQ=255)
                                                                                                                         aGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTAATCGCCCCTTCGAGGCCATTGCCg                   >  2:879524/1‑151 (MQ=255)
                                                                                                                              cTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                >  2:1237605/1‑151 (MQ=255)
                                                                                                                                  cccccAGGTTTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCa               <  1:91268/147‑1 (MQ=255)
                                                                                                                                    cccAGTTCTTTGGTGCAAAAACAACCAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc           <  1:661419/151‑1 (MQ=255)
                                                                                                                                    cccAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc           >  2:717309/1‑151 (MQ=255)
                                                                                                                                     ccAGGTTTTTTGTGGCAAAAAAAAAAGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa             <  1:1014995/149‑1 (MQ=255)
                                                                                                                                     ccAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACg          >  1:291242/1‑151 (MQ=255)
                                                                                                                                      cAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGc         >  1:571902/1‑151 (MQ=255)
                                                                                                                                       aGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGc         >  1:687282/1‑150 (MQ=255)
                                                                                                                                               tGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc  >  1:795725/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAC  >  NC_000913/1222452‑1222739

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