Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,131 C→T 10.7% noncoding (92/120 nt) rrfE → 5S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1310CT10.7% 206.6 / 12.8 74noncoding (92/120 nt)rrfE5S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (4/4);  total (36/38)
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

AGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTA  >  NC_000913/4213018‑4213247
                                                                                                                 |                                                                                                                     
aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTa                                                                                    <  1:708679/149‑1 (MQ=31)
   ggTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGGGTAGGGAACTGCCAGGCATCAAATTAGaa                                                                                 <  3:1255717/149‑1 (MQ=18)
   gctCTGAGAAAACCGGAGCTGCCTGGCGGCAGTAGCGCGGTGAACCCACCTGACCCCACGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGAGGTCGCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATGAGaa                                                                                 <  3:383187/147‑1 (MQ=11)
            aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaa                                                                                >  2:28154/1‑141 (MQ=34)
            aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaa                                                                                <  1:28154/141‑1 (MQ=34)
            aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCgg                                                                        <  3:1175394/149‑1 (MQ=40)
              aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAGCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGTGAGTAGGGAACTGCCAGGCATCAAATTAGaa                                                                                 >  4:1269984/1‑138 (MQ=18)
              aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaaa                                                                              >  2:524413/1‑141 (MQ=35)
              aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaaa                                                                              <  1:524413/141‑1 (MQ=35)
              aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaa                                                                                 <  3:1269984/138‑1 (MQ=33)
                  aaTTTGTCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtctt                                                                  <  2:399315‑M2/149‑30 (MQ=255)
                   aTTTGTCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttc                                                                 <  4:457692‑M2/149‑31 (MQ=255)
                    tttGTCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttcc                                                                <  1:697432‑M2/149‑32 (MQ=255)
                      tGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc                                                                      <  2:1367833/141‑1 (MQ=255)
                      tGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc                                                                      >  1:1367833/1‑141 (MQ=255)
                         cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactg                                                           <  3:359411‑M2/149‑37 (MQ=255)
                         cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaa                                                                                  >  3:1108890‑M2/1‑113 (MQ=255)
                         cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaa                                                                                  <  4:1108890‑M2/126‑14 (MQ=255)
                           ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgag                                                         >  3:457692‑M2/1‑111 (MQ=255)
                            gcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgagc                                                        <  4:641411‑M2/149‑40 (MQ=255)
                              ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcct                                                      >  4:359411‑M2/1‑108 (MQ=255)
                              ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtctttcgactgagcct                                                      >  3:1331794‑M2/1‑108 (MQ=255)
                              ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcct                                                      >  1:399315‑M2/1‑108 (MQ=255)
                               gcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtctttcgactgagcctt                                                     >  1:456454‑M2/1‑107 (MQ=255)
                               gcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcctt                                                     >  2:697432‑M2/1‑107 (MQ=255)
                                 aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaa                                                                                  <  3:72052‑M2/118‑14 (MQ=255)
                                 aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaa                                                                                  >  4:72052‑M2/1‑105 (MQ=255)
                                   tAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgagccttttgt                                                 >  3:641411‑M2/1‑103 (MQ=255)
                                   tAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagc                                                        >  3:1378499‑M2/1‑103 (MQ=255)
                                   tAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagc                                                        <  4:1378499‑M2/142‑40 (MQ=255)
                                     gcgcGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATaa                                                                 >  4:920169/1‑131 (MQ=255)
                                     gcgcGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATaa                                                                 <  3:920169/131‑1 (MQ=255)
                                              tCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcctttt                                                   <  2:752766‑M2/136‑45 (MQ=255)
                                              tCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcctttt                                                   >  1:752766‑M2/1‑92 (MQ=255)
                                                ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaa                                                                                       <  2:1124283‑M2/98‑9 (MQ=255)
                                                ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaa                                                                                       >  1:1124283‑M2/1‑90 (MQ=255)
                                                 cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcc                                                                            >  4:176106/1‑108 (MQ=38)
                                                 cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcc                                                                            <  3:176106/108‑1 (MQ=38)
                                                  aCCTGACCCCATGCCGAACTCAGAAGTGCAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCg                                                           >  2:896910/1‑123 (MQ=255)
                                                  aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggc                                                                          <  4:123090‑M2/109‑22 (MQ=255)
                                                  aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggc                                                                          >  3:123090‑M2/1‑88 (MQ=255)
                                                  aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCg                                                           <  1:896910/123‑1 (MQ=255)
                                                    cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcct                                                      <  3:327972‑M2/127‑42 (MQ=255)
                                                    cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcct                                                      >  4:327972‑M2/1‑86 (MQ=255)
                                                         cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcct                                                      <  3:1428409‑M2/122‑42 (MQ=255)
                                                         cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcct                                                      >  4:1428409‑M2/1‑81 (MQ=255)
                                                          ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtctttcgactgagc                                                        >  4:429407‑M2/1‑80 (MQ=255)
                                                          ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtctttcgactgagc                                                        <  3:429407‑M2/119‑40 (MQ=255)
                                                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcct                                                      >  1:234009‑M2/1‑79 (MQ=255)
                                                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcct                                                      <  2:234009‑M2/120‑42 (MQ=255)
                                                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                                                       >  3:1930/1‑103 (MQ=255)
                                                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                                                       <  4:1930/103‑1 (MQ=255)
                                                                 gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcc                                                       >  4:168154‑M2/1‑73 (MQ=255)
                                                                 gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcc                                                       <  3:168154‑M2/113‑41 (MQ=255)
                                                                     tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgagaggcccagtctttcgactgagc                                                        <  4:178072‑M2/108‑40 (MQ=255)
                                                                     tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgagaggcccagtctttcgactgagc                                                        >  3:178072‑M2/1‑69 (MQ=255)
                                                                       aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcctttt                                                   <  2:747899‑M2/111‑45 (MQ=255)
                                                                       aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcaaataaaacaaaaggcccagtcttccgactgagcctttt                                                   >  1:747899‑M2/1‑67 (MQ=255)
                                                                           gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAc                                                                             >  2:489497/1‑81 (MQ=39)
                                                                           gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAc                                                                             <  1:489497/81‑1 (MQ=39)
                                                                              aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtctttcgactgagc                                                        <  1:1116491‑M2/99‑40 (MQ=255)
                                                                              aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtctttcgactgagc                                                        >  2:1116491‑M2/1‑60 (MQ=255)
                                                                              aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt            <  3:712340/142‑1 (MQ=255)
                                                                              aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt            >  4:712340/1‑142 (MQ=255)
                                                                                aCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcctt                                                     <  4:1410685‑M2/100‑43 (MQ=255)
                                                                                aCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGgcatcaaataaaacgaaaggcccagtcttccgactgagcctt                                                     >  3:1410685‑M2/1‑58 (MQ=255)
                                                                                     gTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaaa                                                                              <  1:411763/70‑1 (MQ=40)
                                                                                     gTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaaa                                                                              >  2:411763/1‑70 (MQ=40)
                                                                                        gcgcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgagcc                                                       >  1:865828‑M2/1‑50 (MQ=255)
                                                                                        gcgcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgagcc                                                       <  2:865828‑M2/90‑41 (MQ=255)
                                                                                                  tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGta  >  3:333314/1‑132 (MQ=255)
                                                                                                  tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGta  <  4:333314/132‑1 (MQ=255)
                                                                                                           gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGCCTTATGGACCGGGGttt                                                     <  1:763016/72‑1 (MQ=17)
                                                                                                           gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGCCTTATGGACCGGGGttt                                                     >  2:763016/1‑72 (MQ=17)
                                                                                                                cccATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgagcctt                                                     >  2:157606‑M2/1‑26 (MQ=255)
                                                                                                                cccATGCGAGAGTAGGGAACTGCCAGacatcagataaaacaaaaggcccagtcttccgactgagcctt                                                     <  1:157606‑M2/68‑43 (MQ=255)
                                                                                                                 |                                                                                                                     
AGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTA  >  NC_000913/4213018‑4213247

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