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

GGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACA  >  NC_000913/1222446‑1222740
                                                                                                                                                        |                                                                                                                                                  
ggCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCg                                                                                                                                                    <  1:1119107/151‑1 (MQ=255)
  cAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaa                                                                                                                                                  <  1:499226/151‑1 (MQ=255)
   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaa                                                                                                                                                  >  1:638654/1‑150 (MQ=255)
      gCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaac                                                                                                                                               >  1:354003/1‑150 (MQ=255)
       cTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaaca                                                                                                                                              >  1:604177/1‑150 (MQ=255)
           cattacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑g                                                                                                                                       <  2:735462/151‑1 (MQ=255)
              tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGTCCAGGTCATTGGAGCGAAAACAACCA‑‑GGcc                                                                                                                                    >  2:343271/1‑151 (MQ=255)
               acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGcc                                                                                                                                    <  1:467864/150‑1 (MQ=255)
                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGcc                                                                                                                                >  2:728378/1‑150 (MQ=255)
                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCa                                                                                                                               <  2:290471/151‑1 (MQ=255)
                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCa                                                                                                                               <  2:1012320/151‑1 (MQ=255)
                       ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTa                                                                                                                            >  2:675022/1‑150 (MQ=255)
                             caaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCaa                                                                                                                      >  2:626595/1‑150 (MQ=255)
                               aTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAAtt                                                                                                                   >  1:65216/1‑151 (MQ=255)
                                ttccggaaaCCTGTGCGCTGTCAGTTTTTCCAATTTCTTTGGTTTCTGCGAAAGTAATATGGTTTTCTGCCGGAAACCCTTTTTAAATTTTGAAGTCCCGCGCCCAGGTAATTGGTGCAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAAtt                                                                                                                   <  1:374096/143‑1 (MQ=255)
                                 ttCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAAtttt                                                                                                                 >  2:434735/1‑151 (MQ=255)
                                   cTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAAtttttt                                                                                                               >  2:949116/1‑151 (MQ=255)
                                       aaCCTGTGCGTTGTCATTTTTTGCCATTTCGTTGGTTTCTGCGCAAGAAATATAGGTTTCTGACGGCAACCCTTTTTAAATATTTTAGTCCCGCCCCAAGGTCTTTGGTGAAAAAACAAACG‑‑GGCCCGCCAGTAAAGCAAAttttt                                                                                                                <  1:525121/146‑1 (MQ=255)
                                       aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                           >  2:1114182/1‑151 (MQ=255)
                                         ccggtgCGCTGTCATTTTTTGCTTTTTCGTTGGTTTCTGCGAAAGTAATATAGTTTTCTGACGACAAACCTTTTTAAATTTTGTAGTCCGGCGCCAGGGTCTTTGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTtcatca                                                                                                          <  2:1226515/147‑1 (MQ=255)
                                              gctcTGTCTTTTTTTGCCTTTTTGTTGTTTTCTGCACAAGTAATATAGGTTTCGGACGGAAAACTTTTTTAAATATTGTAGTCCTGCCCCCAGTTCATTGGTGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                      <  2:113275/146‑1 (MQ=255)
                                              gcgtTGTAATTTTTTCCCTTTTCTTTGGCTTCTGCGAAAGAAATATGGTTTTCGGACGGCAACCCTTTTTAAATTTTGTAGTCCGGCCCCCAGTTCATTGGTGAAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                     <  1:655545/150‑1 (MQ=255)
                                                 ggtgtcgTTTTTGCCCATTTCGTTGCTTTCTGCGAAATTAATATAGGTATCAGACGAAAACCCTTTTAAAATATTGAAGCCCGGCGCCATGGTCATTGTAACAAAAAAAAAAA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                     <  2:495878/142‑1 (MQ=255)
                                                  tGTCTGTTTTTGCCTTTTCGTTGGTTTCTGCGAAAGTAATATAGGTTTCGGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                  <  1:128875/149‑1 (MQ=255)
                                                  tGTCATTTTTTGCCTTTTCGTTGGTTTCTGCGCAATTAATATAGGTTTCTGCCGGCAAACCTTTTTTAATTTTGTAGTCCTGCCCCAAGGTCATTGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                   <  1:446402/148‑1 (MQ=255)
                                                    tCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAtt                                                                                              >  2:954853/1‑151 (MQ=255)
                                                       ttttttGCCGTTTGTTTGGTTTTTGCGCAGTAAATTTTTTTTTCTGGCTGCAACCTTTTTTTAATTTTTCGGTCCTGCACCCGGTTCTTGGGCTATAAAAAACATGGGGCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                            <  2:1230395/149‑1 (MQ=255)
                                                        tttttGCTATTTCGTTGGTTTCTGCGAAATTAATATAGGTTTCTGACGGCAAACCTTTTTTAATATTGGAGTCCGCCCCCAAGTTTTTTGGTGAAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                           <  2:478712/150‑1 (MQ=255)
                                                         ttttGCTTTTTCGTTGGCTTCTGCAAAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCAAGGTCTTGGGTGAAAAAACAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAt                                                                                          <  2:147395/150‑1 (MQ=255)
                                                               ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                    >  2:697199/1‑151 (MQ=255)
                                                               ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                    >  2:44400/1‑151 (MQ=255)
                                                                cATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                   >  2:581841/1‑151 (MQ=255)
                                                                  tttCGTTGGTTTCTGCGCAAGTAATATAGTTTTCGGACGGAAAACTGTTTTAAATTTTGTAGCCCTGCGCCCAGTTCATTGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCa                                                                                  <  1:241501/150‑1 (MQ=255)
                                                                  tttCGTTGGCTTCTGCGCAAGTAATATAGGTTTCGGACGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGGGAAAAAACAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcca                                                                                   <  1:432685/149‑2 (MQ=255)
                                                                       ttGGTTTCTGCGAAAGTAATATGGTTTTCTGACGGCAAACCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCTTTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                             <  1:1168255/150‑1 (MQ=255)
                                                                       ttGGCTTCCGCGAAAGTAATATAGTTTTGTGACGGCAACCCTTTTTAAATATTGAAGCCCGGCCCCCAGGTCATTGGTGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                              <  1:697208/149‑1 (MQ=255)
                                                                        tGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                            >  2:551247/1‑150 (MQ=255)
                                                                         gttttCTGCGCAAGTAATATGGGTTGCTGCCGGAAAACTGTTTTAAATTTTGTAGTCCTGCGCCCAGTTTTTTGGTGCAAAAACAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCAACATCATTCCTCATTCAAGTTTACCCAGAATTg                                                                            <  2:1082188/146‑1 (MQ=255)
                                                                         gttttCTGCACAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTGCCCCCAGTTCTTGGGTGAAAAAACAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                            <  1:1010474/146‑1 (MQ=255)
                                                                         ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                           >  2:712158/1‑150 (MQ=255)
                                                                         ggCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                           >  1:774836/1‑150 (MQ=255)
                                                                                 cgcAAGAAATATGGTTTTTTGACGGCAACCTTTTTTAAATATTTTGGTCCGCCCCCCAGTTCATGGGTGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                                  <  1:343285/151‑1 (MQ=255)
                                                                                    aaGTATTAGAGTTTTCGGCGGGCAACCTTTTTTTAATTTTGTAGTCCTGCCCCAAGTTCATTGGAGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                                <  2:639716/150‑1 (MQ=255)
                                                                                    aaGTAATATGGGTTTCGGACGGCAACCCTTTTTTAATATTGTAGCCCTGCGCCAAGTTCTTTGGTGCAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                                 <  1:85605/149‑1 (MQ=255)
                                                                                    aaGTAATATAGGTTTCGGACGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                                 <  1:812605/149‑1 (MQ=255)
                                                                                               gTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACgg                                                    >  1:21645/1‑151 (MQ=255)
                                                                                                tttCTGCCGTCAACCCTTTTTTAATTTAGTAGCCCCGCCCCCAGGCTTTTGGTGAAAAAACAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACg                                                     <  1:1099375/149‑1 (MQ=255)
                                                                                                 ttCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCTGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCCCCCGTACCGCCAATCACGGc                                                   >  2:290092/1‑150 (MQ=255)
                                                                                                     gACGGCAACCTGTTTTAAATTTTGTAGTCCTGCGCCCAGTTCTTGGGTGAAAAAACAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                              <  1:1074459/151‑1 (MQ=255)
                                                                                                      aCGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                             >  1:963961/1‑151 (MQ=255)
                                                                                                         ggaaaCCTTTTTTAATTTTTGTAGCCCTGCCCCCAGTTCTTGGGGGAAAAAAAAACCG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTa                                            <  1:84334/148‑1 (MQ=255)
                                                                                                          cAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAAt                                          >  1:953735/1‑150 (MQ=255)
                                                                                                           aacctttttttAATATTGTACCCCGGCCCCCAGGTCTTTGGTGAAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAc                                                 <  1:1186200/136‑1 (MQ=255)
                                                                                                             aCCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcc                                      >  1:678665/1‑151 (MQ=255)
                                                                                                              ccGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                     >  1:1019000/1‑151 (MQ=255)
                                                                                                                gTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCt                                   >  2:95991/1‑151 (MQ=255)
                                                                                                                 tttttAATTTTGTAGCCCGGCCCCCAGGTCATTGGTGCAAAAAAAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCt                                   <  1:338378/150‑1 (MQ=255)
                                                                                                                 tttttAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                                  >  2:482875/1‑151 (MQ=255)
                                                                                                                  tttaAATATTGTGGTCTTCCCCCCGGGTTATTGGTCAAAAAAAAAACG‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                                  <  1:1233165/150‑1 (MQ=255)
                                                                                                                   tttAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                >  1:1218155/1‑151 (MQ=255)
                                                                                                                     tAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                >  2:6971/1‑149 (MQ=255)
                                                                                                                      aaTATTGTACCCCGCCCCCCAGCTCTTTGGCGAAAAAAAAAAGG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAg                              <  1:728385/150‑1 (MQ=255)
                                                                                                                                  ccTACGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg                    >  2:269143/1‑149 (MQ=255)
                                                                                                                                   cTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                 >  1:657740/1‑151 (MQ=255)
                                                                                                                                     gcgcCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                 >  1:117707/1‑149 (MQ=255)
                                                                                                                                       gcCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa              >  2:646668/1‑150 (MQ=255)
                                                                                                                                         ccAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACg           >  2:968537/1‑151 (MQ=255)
                                                                                                                                              tCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCATGGCACCTTTAATCGCCCCTTCGAGGCCATTCCCGGTCAGTACGACa        >  2:447896/1‑149 (MQ=255)
                                                                                                                                                aTTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg     >  1:895408/1‑150 (MQ=255)
                                                                                                                                                 ttGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc   >  1:640057/1‑151 (MQ=255)
                                                                                                                                                  tGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc   >  1:292142/1‑150 (MQ=255)
                                                                                                                                                   ggTGCAAAAACAACCA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACa  <  1:832921/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACA  >  NC_000913/1222446‑1222740

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