Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,958,7710AT41.0% 100.0 / 27.3 61D38E (GAT→GAA) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base A (5/31);  new base T (25/0);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGT  >  NC_000913/2958628‑2958916
                                                                                                                                               |                                                                                                                                                          
aGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGCTGACTTTTTCTC‑ttt                                                                                                                                                    >  2:876136/1‑151 (MQ=255)
 gCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCTCCATTCCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑ttt                                                                                                                                                    <  2:1191449/150‑1 (MQ=255)
 gCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTTTTTTTCTC‑tttt                                                                                                                                                   >  1:41087/1‑151 (MQ=255)
   cccTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAc                                                                                                                                                 <  2:653100/151‑1 (MQ=255)
   cccTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAc                                                                                                                                                 <  1:201387/151‑1 (MQ=255)
    ccTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGGTATTTTTAttcccttt                                                                                                                                                    >  1:343954/1‑141 (MQ=255)
    ccTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGGTATCTTTTTTAA‑tttt                                                                                                                                                   >  2:687006/1‑148 (MQ=255)
       ggTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGa                                                                                                                                             <  1:1193506/150‑1 (MQ=255)
        gTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑g                                                                                                                                          <  1:264787/151‑1 (MQ=255)
           cGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGtt                                                                                                                                       <  2:929643/151‑1 (MQ=255)
             ccTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTTTTAC‑tttttttgatgggttt                                                                                                                                       >  2:612668/1‑139 (MQ=255)
                cGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCGGACGTATGGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCt                                                                                                                                   <  1:985460/150‑1 (MQ=255)
                cGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCt                                                                                                                                   <  1:1021819/150‑1 (MQ=255)
                cGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTg                                                                                                                                  <  1:414356/151‑1 (MQ=255)
                  ggATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑aa                                                                                                                               <  1:880427/151‑1 (MQ=255)
                              gACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTCTTTTTTAC‑TTTTTA‑TGAT‑GGTTTTCTT‑GAT‑CGGCT‑GCCAT‑c                                                                                                                >  2:204237/1‑151 (MQ=255)
                              gACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTTTTTC‑TTTTTA‑GGAT‑GGGTTTCtttatccggctccctc                                                                                                                    >  2:811421/1‑136 (MQ=255)
                                    aCGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGttt                                                                                                          <  2:344995/151‑1 (MQ=255)
                                    aCGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGttt                                                                                                          <  2:956968/151‑1 (MQ=255)
                                         caccacCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTCTTTTTTAC‑TTTTTC‑GGAT‑GGGTTTCTG‑AAA‑CGGGT‑GCCAT‑CCCCttttttcc                                                                                                     >  1:338959/1‑146 (MQ=255)
                                            caccaGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATTTTTTTTAC‑TTTTTC‑GGGT‑GGGTTTCTT‑AAT‑CGGCT‑GCCCT‑ccccctttt                                                                                                        >  1:740024/1‑139 (MQ=255)
                                             accaGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAc                                                                                                 <  2:367979/151‑1 (MQ=255)
                                                 gcgcCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTTTTAC‑TTTTTC‑TGAT‑GGGTTTCTT‑AAT‑AGGCT‑GCCAT‑CCCCGTTCTTCCTTAcct                                                                                               >  1:761289/1‑147 (MQ=255)
                                                  cgcCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTTTTAC‑TTTTAT‑AGAT‑GGGTTTCTT‑AAA‑CGGCT‑GCCCT‑CCCCtttttttctgtctt                                                                                               >  1:682271/1‑137 (MQ=255)
                                                  cgcCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAg                                                                                          <  2:509485/151‑1 (MQ=255)
                                                   gcCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTCTTTTTTCC‑TTTTTCAGGAG‑GGTTTCCGT‑AAT‑‑TGCT‑GCCAT‑CCCCTTTTTTTCTGAC‑tt                                                                                              >  1:639225/1‑147 (MQ=255)
                                                           cGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTCTCTTTTTTTC‑CTTTAT‑GGAT‑GGGTTTCTG‑AAT‑TGTCT‑GCCAC‑CCCCCTTCTTGCTGGC‑TT‑TAggggggcc                                                                                   >  2:792051/1‑144 (MQ=255)
                                                                tAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTCTTTTTCAC‑TTTTTT‑GGAT‑GGGTTCCTT‑AAT‑CGGCT‑GCCAT‑CCCCGTTCTTTCTGACTTT‑AGGGGGGCCCAAAgg                                                                              >  1:190080/1‑150 (MQ=255)
                                                                 aaCTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCa                                                                           <  1:226098/151‑1 (MQ=255)
                                                                  aCTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCa                                                                           <  1:911805/150‑1 (MQ=255)
                                                                     gcctgcGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACt                                                                        <  2:533348/150‑1 (MQ=255)
                                                                           ggATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAacaaca                                                                 <  2:352563/151‑1 (MQ=255)
                                                                           ggATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAacaaca                                                                 <  2:330001/151‑1 (MQ=255)
                                                                                  aaaCCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGCTATATTTTTTCC‑ATTTTC‑GGGT‑GGGTTTCTG‑TAA‑CCGCC‑GCCAT‑CCCCTTTTTGCCTTTC‑TT‑TAGGGGGGCCCAAGGGGGAATAAcaacaa                                                                >  1:846504/1‑145 (MQ=255)
                                                                                    aCCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTCTTTTTCAC‑TTTTTC‑TGGT‑GGGTTTTTT‑AAA‑TTGCT‑TCCCT‑TCCCTTTTTTCCTGCC‑TT‑TAGGGGGCCCCAAAGGCAACCAAAAAAAAAAATTCtt                                                        >  2:613332/1‑151 (MQ=255)
                                                                                    aCCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATATTTTTTAA‑TTTTAT‑GGGT‑GGGTTTCTT‑TAT‑ACGGT‑GCCAT‑CCCCTTTTTTCCTTTC‑TT‑AAGGGGGGCCAAAAGGGAACTAaaaaaaaaaaaactt                                                        >  2:395084/1‑146 (MQ=255)
                                                                                      cAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTg                                                      <  2:204385/151‑1 (MQ=255)
                                                                                        gCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGaa                                                    <  2:909275/151‑1 (MQ=255)
                                                                                           aGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATTTTTTTCAC‑TTTTTT‑GGAT‑GGGTTTCTG‑AAT‑CGGCT‑GCCAT‑CCCCTTTCTTCCTTACTTT‑AGGGGGGCCCAAAGGG‑AACTAaaaaaaaaaaatattttaa                                                    >  1:774755/1‑139 (MQ=255)
                                                                                           aGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCa                                                 <  1:1124179/151‑1 (MQ=255)
                                                                                            gACCACCATACCGTTATCCAGTCGTATAGCCTGATACTTGCGGTTATCTTTATCGC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAG‑CGGCT‑TCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGTTGCGCAAAGGGCAACTAACAACAATGATGCTTTGaa                                                    <  2:430014/147‑1 (MQ=255)
                                                                                            gACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATTTTTTTTAC‑CTTTTT‑GGGT‑GGGTTTCTT‑AAT‑CGGCT‑GCCCTCCCCTTTTTTCCCTGCT‑TT‑AAGGGGGCCCAAA‑GGCAACAAACAACAATAATACttt                                                       >  1:846684/1‑144 (MQ=255)
                                                                                              ccaccaTACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGt                                              <  1:953527/151‑1 (MQ=255)
                                                                                               caccaTACCGTTATCCAGACGTATAGCCTGATACTGGCGGGTTTTTTTTTTAA‑TTTTTC‑GGGTGGGTTTTTTT‑AAT‑‑GGCTGGCAAC‑CCCCTTTTTTCCGAGT‑TT‑AAGGTGGCCCAAAGGG‑AACTAACAAAAATAATTCTTTTTACCaaggg                                             >  1:86911/1‑147 (MQ=255)
                                                                                                  caTACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTTTTTTTTTTAA‑TTTTTT‑CGGT‑GGGTTTCTT‑TAT‑TGCCC‑GCCCT‑CCCCTTTTTTCCTGAC‑TT‑AAGGGGGGCCCAAGGGCCACTAAAAACAAAAATTCTTTTAACCacggg                                             >  2:500189/1‑144 (MQ=255)
                                                                                                     aCCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCg                                        <  2:100831/150‑1 (MQ=255)
                                                                                                        gTTATCCAGACGTATAGCCTGATACTGGCGGTTTTTTTTTTTACTTTTTTC‑GGAG‑GGTTT‑CTT‑AAT‑AGGCC‑GGCCT‑CCCCTTTTTTCCTTTCTTT‑AAGGGGGCCCAAAGGCAACCAAACACCAAAAAAATTTTTAACAAGG‑GCTGCggggg                                    >  2:232567/1‑150 (MQ=255)
                                                                                                          tATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCa                                   <  2:244307/150‑1 (MQ=255)
                                                                                                           aTCCAGACGTATAGCCTGATACTGGCGGTTTTCTTTTTTAC‑TTTTTC‑TGAT‑GGGTTTCTG‑AAT‑CGGCT‑GCCCT‑CCCCTTTTTTCCTTACTTT‑AGGGGGGCCCAAAGGG‑GACTAAAAAAAAAAATGCTTTTTAACAAGGGGTGCggggggtt                                 >  1:503113/1‑147 (MQ=255)
                                                                                                              cAGACGTATAGCCTGATACTGGCGGTTATTTTTTTTAC‑TTTTTT‑GGAA‑GGGTTTCTT‑GAT‑CGGCT‑GCCCT‑CCCCCTTTTTCCCTAC‑TTAAAGGGGGCCCAAAGGG‑AACCAACAAAAAAAAAAGTTTTTAACAGGGGCCGCggggggtttac                              >  1:830486/1‑144 (MQ=255)
                                                                                                                 aCGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGa                           <  2:935045/151‑1 (MQ=255)
                                                                                                                  cGTATAGCCTGATACTGGCGGTTATTTTTATTAA‑TTTTTT‑GGGT‑GGGTTTCTGAAAA‑TGGCG‑GCCTT‑CCCCTTT‑TTTCTGAC‑TT‑AAGGGGGGCCCAAGGGGAAATAACAAAAAAAAAAATTTTAACAAAGTGGTGCGGGGGAtt                                 >  2:3004/1‑144 (MQ=255)
                                                                                                                    tataGCCTGATACTGGCGGTTATTTTTTTTAC‑TTTTTC‑GGGT‑GGGTTTCTT‑TAT‑TGGCT‑GCCCT‑CCCCGTTTTTCCTGTCTTT‑AAGGGGGCCCAAAGGG‑AAATAAAAAAAAAAATAATTTTTAACAGGTGGTGCGGGGGGTTCCCCGAccc                        >  2:1014455/1‑150 (MQ=255)
                                                                                                                       aGCCTGATACTGGCGGTTATCTTTTTTAC‑TTTTTC‑AGCT‑GGGTTTCTT‑AATCCGGGG‑GCCTC‑CCCCTTTTTTCTTACT‑TT‑AAGGGGGCCCAAAGGG‑AAATAACAACAAAAAATCTTTTAAACACGGGGTTCggggggtttac                              >  2:918194/1‑135 (MQ=255)
                                                                                                                        gCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATa                    <  2:809293/151‑1 (MQ=255)
                                                                                                                                 tGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACtt            >  1:139415/1‑150 (MQ=255)
                                                                                                                                  ggCGGTTATCTTTATCAC‑TTTTGC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTGCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTa           <  2:343943/150‑1 (MQ=255)
                                                                                                                                       ttATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTc     >  1:954327/1‑151 (MQ=255)
                                                                                                                                         aTCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCt    <  2:630934/150‑1 (MQ=255)
                                                                                                                                         aTCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCt    <  1:17258/150‑1 (MQ=255)
                                                                                                                                          tCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGt  >  1:837522/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                          
AGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCACCACCAGCGCCGAGAGCGATTTAACTGCCTGCGGATCAGAAACCAGCAAGACCACCATACCGTTATCCAGACGTATAGCCTGATACTGGCGGTTATCTTTATCAC‑TTTTAC‑GGAT‑GGTTTCCTG‑AAT‑CGGCT‑GCCAT‑CCCGTTTCTGCCTGAC‑TT‑AAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGT  >  NC_000913/2958628‑2958916

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