Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,269,969	0	A	C	24.6%	71.5 / 15.0	64	intergenic (+367/+247)	yhaC/rnpB	pentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (14/29);  new base C (5/9);  total (19/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACT  >  NC_000913/3269822‑3270112                                                                                                                                                    |                                                                                                                                                gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                                                                                                                <  1:594039/149‑1 (MQ=255)     ataGCGCACTCAACAATCTGCAAGCCGCCGGATCGGTTTATTACGTACATCAGGTAATACTGACCGATAAGCCGCTTGCTTTTGGGTATAGTGTAGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCTCCGAAGCgg                                                                                                                                            <  2:140195/147‑1 (MQ=18)     aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                            <  1:1698149/149‑1 (MQ=255)     aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                            <  1:1515893/149‑1 (MQ=255)     aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAACCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                            <  1:467616/149‑1 (MQ=255)      aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGAAAAAACCCGCCGAACCggg                                                                                                                                           >  1:936917/1‑149 (MQ=255)        gcgcACTCAACAATCTGAAACCCGCAGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                         <  1:1503657/149‑1 (MQ=255)        ccgcACTCAACAATCTGAAACCCGCCGTAACGGTGTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGTGTAAAGTGTCGTGTACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                         <  2:1701343/148‑1 (MQ=18)           cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTGCTGCCCCTCCAAAAGTAAAAACCCGCGTCGGCGGGttttt                                                                                                                                      <  1:1105021/149‑1 (MQ=18)            aCTCAACAATCATAAACCCGCCGGATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttg                                                                                                                                     <  2:446653/149‑2 (MQ=11)             cTCACCAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGAATAGTGTCGTGGACAGTCATTGATCTTTCTGCCCCTCCAAACGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                    <  1:1533113/149‑1 (MQ=255)             cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAACACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTAc                                                                                                                                    <  1:1218406/149‑1 (MQ=18)              tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGACGTGGACAGTCATTCATCTTTCTGCCCCAGCACGAATAAAAACACGCCGAAGCGCGTTTATAcc                                                                                                                                   >  1:1359625/1‑148 (MQ=18)               caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTActt                                                                                                                                  <  2:919905/149‑3 (MQ=255)                aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTActtt                                                                                                                                 <  2:1310812/149‑4 (MQ=17)                aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCATCAGCGGGtttttgctta                                                                                                                                 >  2:1734532/1‑144 (MQ=17)                 acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTGTTACGtta                                                                                                                                >  1:887362/1‑147 (MQ=255)                   aaTCTGAAACCCGCCGGAACGGTTTATTCCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCTGGTTTTTACGTaaaa                                                                                                                              <  2:1866522/149‑1 (MQ=255)                     tCTGAAAGCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGTTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTTAAACa                                                                                                                            <  2:1739067/149‑1 (MQ=255)                      cTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                           <  2:1470466/149‑1 (MQ=255)                      cTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                           <  1:1542022/149‑1 (MQ=255)                       tGAAACCCTCCGGAACTGTCTACTACGTACATCAGTTAAAACCGCCCGATAACCCGCTTTCATATGGGTATAGTGTCGTGGACAGTAATTCATCTTTCTGCCCCTCCAAAAGTAACAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                          <  2:369137/149‑1 (MQ=17)                       tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                          <  1:564255/149‑1 (MQ=255)                        gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                         <  1:897459/149‑1 (MQ=255)                        gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                         <  2:54255/149‑1 (MQ=255)                        gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                         <  2:1123062/149‑1 (MQ=255)                            cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGAGCACACTCATTAATCGTAGTGCCCCGCCAATAGTATAAACCCGCCGCAGCGGGATTTTAAGTAAACCAGGGGaaa                                                                                                                     >  1:1167055/1‑149 (MQ=14)                            cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGAAATGTGTAGAGCACACTGATTAAATAATGTGACCCTACAAAAGAAAAAACCCGCGTAAGCGTGTTTTTAAGTAAAGCAGGGGagg                                                                                                                     >  1:37145/1‑147 (MQ=9)                             ccgccgGAACGGTTTATTACGTAAATCAGGTAAAACTTACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCCTGGACAATAATTCATATTTCTGCCCCTCCAAAAAAAAAAACCCGACGAAGAGGGTTTTTACGTAAAACAGGTGaaaa                                                                                                                    >  2:1768500/1‑148 (MQ=14)                                 cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                                <  1:1727318/149‑1 (MQ=255)                                 cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                                <  1:896828/149‑1 (MQ=255)                                  ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                               <  2:1271736/149‑1 (MQ=255)                                             ttaCGTACATCACTTATAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATACATCTTTCTCCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                                                                                                    <  2:175423/149‑1 (MQ=14)                                                   aCATCAGGTAACACTGACCGATAATCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                              <  1:1390407/149‑1 (MQ=255)                                                     aacAGTTAACACTGACCGTTAAGCCGCTTTCTTTTGGGGATAGTGTCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                            <  1:740669/147‑1 (MQ=11)                                                      tCATGTAAAACTGTCCGATAAGCCGCTTTCTTTTGGTTTTAGTGGCGGCGACAGGACTTCATCTTTCAGCACCTCAAAAAGCACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                           <  1:1307610/149‑1 (MQ=11)                                                              aaCTGACCGATAAGCCGCTTTCTTTTGCGTATCGTCTCGAGGACAGTCATTCATCTATCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtg                                                                                   <  1:542382/149‑1 (MQ=11)                                                                 tttccGATAAGGAGCTTGCATTAGTGTATTTTGTCGAGGCAAGTGTTTAATCTTTCTGCCCCTCCAAAAGTTCAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCg                                                                                <  1:862084/146‑1 (MQ=11)                                                                           gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTAAGCCCCTACAAAAGAAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGAACCATAAGCCGCGTTCTGTAGGGGATAGTGTTGTGGagggaag                                                                      >  1:642992/1‑143 (MQ=11)                                                                                        gcgtttaGTGTCGTCGACAGACATTCAGATTTATCCCACTCAACAATTAAAAACCCGCCGAAACGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                         <  1:1604059/144‑1 (MQ=11)                                                                                               gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                           >  2:1075865/1‑108 (MQ=35)                                                                                                tgtCGTGGCCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                          >  1:1782609/1‑124 (MQ=38)                                                                                                          aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                    >  1:827359/1‑120 (MQ=39)                                                                                                          aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                    <  2:827359/120‑1 (MQ=39)                                                                                                             aattcttGTTTTTTCCCGTCCTAAATTATAAACCCGCCGTACCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                    <  1:1789548/148‑1 (MQ=11)                                                                                                              attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCgggtt                                                                                                  >  1:1242334/1‑82 (MQ=255)                                                                                                               ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattc                                                                   <  2:115604/116‑1 (MQ=17)                                                                                                               ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGtcattc                                                                   >  1:115604/1‑116 (MQ=25)                                                                                                                       acagCCCCACCAAAATAAAAAACCCGCCGAACCGGGTTATTACGTAACACAGGTCAAACTGACCGATCACCCGCTTTCATTTGGGTAAAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCAGCGGCggg                          <  1:845050/146‑1 (MQ=11)                                                                                                                         tcccacTCAAAAATTATAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTATCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                        <  1:587973/144‑1 (MQ=11)                                                                                                                          gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                            >  2:320727/1‑64 (MQ=255)                                                                                                                           ccccTCCAAAAGTATAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttt                      <  1:1793879/149‑1 (MQ=11)                                                                                                                           ccacTCCACAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttt                      <  1:222281/146‑1 (MQ=11)                                                                                                                            cacTCAAAAATTTAAAACCCGCCGTAGCGGGTTATTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                     <  2:1294893/147‑1 (MQ=11)                                                                                                                            cacTCAAAAATTATAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                     <  2:936917/147‑1 (MQ=11)                                                                                                                                 tcaaaGCATAAACCCGCCGCATCGGGTTATTACGTTAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgcttt                <  2:140686/147‑6 (MQ=11)                                                                                                                                  aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGGATAGTGTCGTGGCCAGTCATTCATGATAGTTCACCGCCAACACCAAAAACCCGGCGCAGCGCGTTGTTAAGtga               >  1:676128/1‑147 (MQ=11)                                                                                                                                    aagcttaaaCCCGCCTTATCGGGTTTTTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc             <  1:363557/143‑1 (MQ=11)                                                                                                                                      gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                               >  1:1769183/1‑113 (MQ=38)                                                                                                                                      gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                               <  2:1769183/113‑1 (MQ=14)                                                                                                                                       tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTAGTTGAAAGTCATTCATATTTCTGCCCCAACAAAAGCAAAAAACCGCCGAAGCGGGTTTTTTCGTTAACCAg          >  2:1458654/1‑149 (MQ=11)                                                                                                                                       tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGAACGATAAGCCGCTTACTTTTGGGTATAGAGTAGTGGACAATGATTGATCTTTATTCAACTAAAAAATCAAAAACCCGCCAAAAAGAGATTTTAAGTAAATCAg          >  2:412873/1‑149 (MQ=255)                                                                                                                                         aaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCgg                                                                                                     >  2:1319547/1‑55 (MQ=11)                                                                                                                                         aaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCgg                                                                                                     <  1:1319547/56‑2 (MQ=11)                                                                                                                                               cgccgTACCGGTTTTTTACGTACAACAGGTTAACCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt  <  2:1868284/149‑1 (MQ=11)                                                                                                                                                    |                                                                                                                                                GAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACT  >  NC_000913/3269822‑3270112

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