Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG64.1% 51.1 / 78.7 92intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (16/17);  new base G (28/31);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTC  >  NC_000913/1286835‑1287103
                                                                                                                                  |                                                                                                                                          
gCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                                    >  1:236498/1‑139 (MQ=37)
  aGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                  <  2:236498/139‑1 (MQ=255)
              tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCg                                                                                                                      >  1:666129/1‑139 (MQ=255)
              tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCg                                                                                                                      >  2:1507450/1‑139 (MQ=255)
               gCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa                                                                                                                     <  2:666129/139‑1 (MQ=255)
                cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGaa                                                                                                                    <  2:547161/139‑1 (MQ=255)
                              aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTc                                                                                                      <  1:1507450/139‑1 (MQ=255)
                                     gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCAt                                                                                               >  2:371186/1‑139 (MQ=255)
                                     gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCAt                                                                                               <  1:371186/139‑1 (MQ=255)
                                     ggacTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCAt                                                                                               <  2:564462/136‑1 (MQ=0)
                                        cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                            <  2:1498883/139‑1 (MQ=18)
                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGTGCTCGAGCCGAACCTTAATCGAAGCTtc                                                                                                               <  1:4440/118‑1 (MQ=11)
                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGTGCTCGAGCCGAACCTTAATCGAAGCTtc                                                                                                               >  2:4440/1‑118 (MQ=255)
                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  1:800498/1‑139 (MQ=16)
                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  2:628769/1‑139 (MQ=37)
                                          taaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  2:1260787/134‑1 (MQ=11)
                                           ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa                                                                                                                     <  1:664228/111‑1 (MQ=255)
                                           ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa                                                                                                                     >  2:664228/1‑111 (MQ=255)
                                            tGAAAGTGATGGTGGTGGGGGAAGGTTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                        <  2:381835/139‑1 (MQ=12)
                                                  cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTtc                                                                                                               >  1:189827/2‑110 (MQ=255)
                                                  cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTtc                                                                                                               <  2:189827/109‑1 (MQ=255)
                                                   gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                 >  1:1196551/1‑139 (MQ=12)
                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  1:1218224/122‑1 (MQ=12)
                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  2:1218224/1‑122 (MQ=12)
                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                         <  1:903711/139‑1 (MQ=255)
                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                        >  1:653835/1‑139 (MQ=255)
                                                            gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                                        <  2:954995/139‑1 (MQ=16)
                                                             ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAGGCTTCTCATCCTTCCCCGc                                                                                                 >  1:500495/1‑113 (MQ=255)
                                                             ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAGGCTTCTCATCCTTCCCCGc                                                                                                 <  2:500495/113‑1 (MQ=255)
                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                                           >  1:1411962/1‑99 (MQ=255)
                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                                           >  1:1316275/1‑99 (MQ=255)
                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                                           <  2:1411962/99‑1 (MQ=255)
                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                                           <  2:1316275/99‑1 (MQ=255)
                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                   >  1:1369578/1‑139 (MQ=255)
                                                                  aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCt                                                                                                                 <  1:893071/92‑1 (MQ=255)
                                                                  aGGACTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCt                                                                                                                 >  2:893071/1‑92 (MQ=255)
                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                <  1:1340623/139‑1 (MQ=255)
                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCgttggg                                                                <  1:424455/139‑6 (MQ=0)
                                                                       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGa                                                                                                                     >  1:1292603/1‑83 (MQ=255)
                                                                       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGa                                                                                                                     <  2:1292603/83‑1 (MQ=255)
                                                                       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGTAGCTTCTCATcc                                                                                                         <  2:718622/95‑1 (MQ=11)
                                                                       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGTAGCTTCTCATcc                                                                                                         >  1:718622/1‑95 (MQ=11)
                                                                         cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                           >  2:133206/1‑139 (MQ=255)
                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTc                                                          <  2:1369578/139‑1 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCa                                                                                                            <  2:1367067/85‑1 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCa                                                                                                            >  1:1367067/1‑85 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                            <  1:1127064/101‑1 (MQ=9)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                            >  2:1127064/1‑101 (MQ=12)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAACCGAACCTTAATCGAAGCTTCTCAt                                                                                                           >  1:49377/1‑86 (MQ=11)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAACCGAACCTTAATCGAAGCTTCTCAt                                                                                                           <  2:49377/86‑1 (MQ=11)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCa                                                                                                            <  2:766019/85‑1 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCa                                                                                                            >  1:766019/1‑85 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                            >  1:29268/1‑101 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                            <  2:29268/101‑1 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                                           >  2:659418/1‑81 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                                           <  1:659418/81‑1 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                                                 >  2:664176/1‑91 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                                                 <  1:664176/91‑1 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  1:723542/1‑98 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  2:723542/98‑1 (MQ=255)
                                                                                       cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                             <  2:1521742/139‑1 (MQ=255)
                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATc                                                                                                          >  1:877735/1‑69 (MQ=11)
                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTGAATCGAAGCTTCTCATc                                                                                                          <  2:877735/69‑1 (MQ=11)
                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                                                 <  1:1321773/78‑1 (MQ=255)
                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                                                 >  2:1321773/1‑78 (MQ=255)
                                                                                                  gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggt                                 <  2:653835/139‑1 (MQ=11)
                                                                                                   ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggt                                 <  2:361997/139‑1 (MQ=11)
                                                                                                    gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTtc                                                                                                               >  1:811871/1‑60 (MQ=255)
                                                                                                    gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTtc                                                                                                               <  2:811871/60‑1 (MQ=255)
                                                                                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                                                 >  1:1418128/1‑73 (MQ=255)
                                                                                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                                                 <  2:1418128/73‑1 (MQ=255)
                                                                                                         tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                       <  1:642850/95‑1 (MQ=255)
                                                                                                         tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                       >  2:642850/1‑95 (MQ=255)
                                                                                                              gCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGaa                                                                                                                    <  1:1227970/45‑1 (MQ=14)
                                                                                                              gCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGaa                                                                                                                    >  2:1227970/1‑45 (MQ=11)
                                                                                                              gCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCtt                                                                                                                <  1:615089/49‑1 (MQ=14)
                                                                                                              gCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCtt                                                                                                                >  2:615089/1‑49 (MQ=11)
                                                                                                              gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  1:1107031/71‑1 (MQ=255)
                                                                                                              gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  2:1107031/1‑71 (MQ=255)
                                                                                                                    aCGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  2:795352/65‑1 (MQ=11)
                                                                                                                    aCGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  1:795352/1‑65 (MQ=11)
                                                                                                                     cGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  1:365037/1‑64 (MQ=11)
                                                                                                                     cGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  2:365037/64‑1 (MQ=11)
                                                                                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          >  1:500815/1‑64 (MQ=255)
                                                                                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                          <  2:500815/64‑1 (MQ=255)
                                                                                                                         ctctcGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgaaccttcg           >  2:1134431/1‑129 (MQ=11)
                                                                                                                            tcGCTGGCGCTCGAGCCGAACCTTAATCGAAGCtt                                                                                                                >  1:892279/1‑35 (MQ=14)
                                                                                                                               cTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                          <  1:911543/86‑1 (MQ=255)
                                                                                                                               cTGGCGCTCGAGCCGAACCTTAATCGAAGCTTATCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                          >  2:911543/1‑86 (MQ=17)
                                                                                                                               cTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCTGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGtcgcgc     >  2:1261368/1‑139 (MQ=11)
                                                                                                                                 gACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTccc                                                                                                    <  1:1280605/42‑1 (MQ=255)
                                                                                                                                 gACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTccc                                                                                                    >  2:1280605/1‑42 (MQ=255)
                                                                                                                                  aCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc  <  1:658205/139‑1 (MQ=11)
                                                                                                                                  |                                                                                                                                          
GCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTC  >  NC_000913/1286835‑1287103

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