Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,624,5170GT41.4% 75.5 / 26.6 58E492* (GAA→TAA) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base G (7/27);  new base T (24/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 2.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC  >  NC_000913/2624375‑2624666
                                                                                                                                                |                                                                                                                                                       
aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCa                                                                                                                                                 <  1:55950/151‑1 (MQ=255)
aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCa                                                                                                                                                 <  1:662706/151‑1 (MQ=255)
aGAAAACGGTAAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCa                                                                                                                                                 <  1:547517/151‑1 (MQ=255)
        gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTcc                                                                                                                                       <  1:339463/151‑1 (MQ=255)
        gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTcc                                                                                                                                       <  2:370316/151‑1 (MQ=255)
          gAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCg                                                                                                                                      <  1:86193/150‑1 (MQ=255)
                  gCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACAttt                                                                                                                             <  1:342450/151‑1 (MQ=255)
                    gTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGa                                                                                                                           <  1:305043/151‑1 (MQ=255)
                      tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTAATTTTTTTAAAACCCAT‑A‑CCGccccggggcgttttttt                                                                                                                         >  2:377007/1‑137 (MQ=255)
                      tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTAATTTTTTTAAAAACCAT‑A‑ACGTCCCGTGGCAttt                                                                                                                             >  2:925342/1‑147 (MQ=255)
                       aTGCACACATCGGGGCGGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACTAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGatta                                                                                                                        <  2:898286/151‑1 (MQ=255)
                       aTGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCTCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGatt                                                                                                                         <  1:823548/150‑1 (MQ=255)
                         gCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTTTTTAAAAACCAT‑A‑CCGTCCGGTGGCATTTGattatt                                                                                                                      >  2:972903/1‑151 (MQ=255)
                          cacaccTGGTGACCGGGACCTTTAACGAAAAATCCGCGCGTCTTTATACTGACTATTCGTTGCTGCCCGCCGATGCGCGCATCCCCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGattatt                                                                                                                      <  1:713750/150‑1 (MQ=255)
                                  gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTAATTTATTTGAAAACCAA‑A‑CCGTCCGGGGGCCTTTGATTTTTTTATTGTAt                                                                                                             >  1:835697/1‑151 (MQ=255)
                                  gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTATTTACAACCCAT‑A‑CCGTCCGGGGACATTTGATTATTTTATTGt                                                                                                               >  1:172520/1‑149 (MQ=255)
                                   ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATc                                                                                                            <  1:940156/151‑1 (MQ=255)
                                   ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGT‑TTTTAATTTTTTTTAAACCCCT‑A‑CCCTCCGGGGACATTTGTTTATTTTATgg                                                                                                                >  2:919331/1‑147 (MQ=255)
                                    gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAAACCAC‑A‑CCGTCCGGGGGCATTTGATTATTTTATGGTATcg                                                                                                           >  2:107799/1‑151 (MQ=255)
                                      ccGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATcgcc                                                                                                         <  1:433594/151‑1 (MQ=255)
                                         ggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCa                                                                                                      <  1:546023/151‑1 (MQ=255)
                                           aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCaaa                                                                                                    <  2:654617/151‑1 (MQ=255)
                                           aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAACCCAT‑A‑CCGTCCCGGGGCATTTGATTATTTTATTGTATCGCCGCaa                                                                                                     >  2:493947/1‑150 (MQ=255)
                                             cTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTAATTTTTTTTAAAACCC‑T‑A‑CCGGCCGGGGGGAATTTTTTTTTTAATTGTATTGCCCCcaaa                                                                                                   >  2:386182/1‑150 (MQ=255)
                                              tttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTAATTTTTTTAAAAACCCC‑A‑CCCTCCGGGGGCATTTTATTTTTTTATGGGATCGCCGCcaaa                                                                                                   >  2:894849/1‑149 (MQ=255)
                                                tAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTc                                                                                                <  1:813317/150‑1 (MQ=255)
                                                 aaCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTccc                                                                                              <  2:304142/151‑1 (MQ=255)
                                                   cGAAAAAACCGCGCGTCTTTAGTCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCTCCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgc                                                                                            <  1:792911/151‑1 (MQ=255)
                                                     aaaaacCCCCGCGTCTTTATACTGACTATCCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgcc                                                                                           <  1:220282/150‑1 (MQ=255)
                                                            cgcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTTTTTAAAACCCCT‑A‑CCCTCCGGGGGCACTTTTTTTTTTTATTGTATCTCCCCAAAACcccccccccccctt                                                                                    >  1:322217/1‑146 (MQ=255)
                                                             gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTTTTGAAAACCCCC‑A‑CCCGCCGGGGGCATTTTTTTTTTTTATGGTATCGCCGCAAAACTCCCGCCGCCTCtt                                                                                    >  1:415017/1‑149 (MQ=255)
                                                              cgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAAACCATAA‑CCGCCCGG‑GGCATTTGATTATTTTATGGTATCCCCGCCAAAATCcccccccccttttt                                                                                  >  2:374856/1‑144 (MQ=255)
                                                               gcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTAt                                                                                <  1:95861/151‑1 (MQ=255)
                                                                 gTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTT‑ATTCAATTTATTGAAAAACCAT‑ACCCCCCCGGCGA‑AATTGCTTTTTTTATGGTTTAACCCCCCAACTCCCGCCCCCTACTGttt                                                                                >  2:805673/1‑147 (MQ=255)
                                                                      tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGTGG‑TTTTTATTTTTTGAAAAACCAC‑A‑CCCTCCCGGGACATTTTTTTTTTTTATTGGATTTCCCCAAAAAACCCCCCCCCTAtttttttgaa                                                                            >  1:563421/1‑140 (MQ=255)
                                                                          cTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTT‑TTTACCTTTTTTTAAAACCCCA‑C‑CCCCCCGGGGCAATTTG‑TTTTTTTATGGGTTCGCCCCAAAACTCCCCCCCCCTAttttttttaatgggggg                                                                     >  1:16318/1‑136 (MQ=255)
                                                                               tATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTTTTTAAAACCCAA‑A‑CCGCCCGGGGGCATTTGATTATTTTATGGTATCGCCGCCAAAATCCCCCCCCCCATTGTTTGAAATGGTGGGcc                                                                   >  1:969369/1‑148 (MQ=255)
                                                                                     ttGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGG‑TTTTACTTTTTTTAAAACCCCT‑A‑CCCCCCGGGGGCATTTTATTTTTTTATTGTATCGCCGCAAAACTCCCCCCGCCCCTTTTTTGAAATGGGGGGCCGCGAGGTc                                                           >  1:533587/1‑150 (MQ=255)
                                                                                       gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTAATTTTTTTAAAAACCAA‑A‑CCCTCCGGGGGCATTTGTTTATTTTATGGTATCGCCCCCAAACTCCCCCCCCCCATTGTATGAAAAGGGGGACCGCGGGAACGcc                                                        >  2:624779/1‑151 (MQ=255)
                                                                                         tGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTATTTAAAAACCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTTTATGAAATGGTGGACCGCGAGATCGcc                                                        >  2:897135/1‑149 (MQ=255)
                                                                                          gACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTACTTTTTTTAAAACCCAT‑A‑CCGCCCGGTGGCATTTGATTTTTTTATTGTATTGCCGCCAAACTCCCCCCCCCTATTGTATGAAAATGTGGACCGCGGGGTCGCCCAc                                                     >  1:895130/1‑151 (MQ=255)
                                                                                          gACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCaa                                                      <  1:7312/150‑1 (MQ=255)
                                                                                           aCCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTACTTTATTTAAAACCCAA‑A‑CCGCCCGGTGGCATTTTATTTTTTATTGGAATCGCCGCCAAAATCCCGCCCCCCATTTTATAAAAATGTGGGACCCGAGACCGccca                                                      >  2:104738/1‑147 (MQ=255)
                                                                                             cgccgATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAAACCAT‑A‑CCCTCCGGGGGCATTTTATTTTTTTATGGGATCGCCGCAAAACTCCCCCCGCCCCTTGTTTGAAATGGGGGACCGCGGGATCGCCCAccc                                                   >  2:655918/1‑148 (MQ=255)
                                                                                                   tGCGCGCATCACCAACGAAGTACGGC‑GGGGT‑GTTAACTTTTTTTAAAAACCAT‑A‑CACTGCGGGGGGAATTTATTATTTAATTGGATCTGCCCAAAAAACCCGGCCCCCCTTTTATGAAAaggggggacccggggttg                                                          >  2:886663/1‑125 (MQ=255)
                                                                                                     cgcgcATCACCAACGAAGTACGGC‑GGGTA‑TTTTAATTTTTTTAAAAACCAT‑A‑ACGTCCGGGGGCATTTGTTTTTTTTATAGGATCGCCGCAAAAATCCCGCCCCCTCTTGTATGAAATGGGGGACCGCGAGAACGCCAACGCGCAGCaa                                            >  2:401489/1‑149 (MQ=255)
                                                                                                        gcATCACCAACGAAGTACGGC‑GGGTG‑TTTTATTTTTTTTTAAACACAA‑A‑CCCTCCCGGGGCATTTGGTTATTTAATTGTATTGCCGCAAAACTCCCCCCCCCCATTTTTTTAAAATGGGGAGCGCGGGGGCGCCAACCCCCCGCAAggg                                         >  1:624026/1‑149 (MQ=255)
                                                                                                         cATCACCAACGAAGTACGGC‑GGGTA‑TTTTACTTTTTTTAAAAACCAT‑A‑ACGTCCCGGGGCATTTTATTTTTTTATTGTATCGCCCCAAAAAACCCCCCCCCCCTTTTATGAAAAGGTGGACCCCGAGGTCGCCCACCCCCCGCAAggg                                         >  2:603726/1‑148 (MQ=255)
                                                                                                              ccAACGAAGTACGGC‑GGGGA‑TTTAACTTTATTGTAAAACCAT‑A‑CCGCCCGGGGGCATTTGATTATTTTATTGTATTGCCGCAAAAATCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCCGCAAGGGCTGccc                                   >  1:74359/1‑149 (MQ=255)
                                                                                                               cAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGt                                <  1:932659/151‑1 (MQ=255)
                                                                                                                aaCGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTg                               <  2:698068/151‑1 (MQ=255)
                                                                                                                   gAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTa                            <  2:649177/151‑1 (MQ=255)
                                                                                                                     aGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATc                          <  2:733779/151‑1 (MQ=255)
                                                                                                                              gggTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACAGTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGaa                  <  2:585908/151‑1 (MQ=255)
                                                                                                                                       aaCTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAAt           >  2:312962/1‑150 (MQ=255)
                                                                                                                                        aCTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATa          >  2:71426/1‑150 (MQ=255)
                                                                                                                                          tttATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAAcc       >  2:180087/1‑151 (MQ=255)
                                                                                                                                               tGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGt   <  2:449433/150‑1 (MQ=255)
                                                                                                                                                gAAACCCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTc  <  1:226399/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
AGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC  >  NC_000913/2624375‑2624666

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