Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,646,108 T→G 27.8% V242G (GTG→GGG)  rlmJ → 23S rRNA m(6)A2030 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,646,1080TG27.8% 29.2 / 3.3 43V242G (GTG→GGG) rlmJ23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (16/14);  new base G (5/8);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACT  >  NC_000913/3645893‑3646310
                                                                                                                                                                                                                        |                                                                                                                                                                                                           
aaCTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGGGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGGACTTCCCGACAGGGAGCCCCGGGG‑CATGACCGCTTCCGGCCGGATTGTGATTAAAcccccg                                                                                                                                                                                                >  2:16915/1‑229 (MQ=255)
               gtttgtCAGGGGGATAGCAGAAGGTTACAACCTTTTCGCCAGTGGTATTTACGCATGGTGGTATCGGGGGGGTGTGGGTAAGAAAATTAAGCGCATGATCAACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGAAAATTGAACTGGCGGTACTGCCAGACAGCGACCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGaacaa                                                                                                                                                                                 <  2:90009/225‑1 (MQ=255)
                   gtCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCCGGGAAACTGGaaccac                                                                                                                                                                                >  1:555172/1‑226 (MQ=255)
                    tcggggggATAGCGAAAGGTTAAAACCGTTTCGCCATTGGTATTTACGCACGGGGGTATCGGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCGGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGGGATTAACCCGCCGTGGAAACTGGAACAACAGAt                                                                                                                                                                            <  1:479548/224‑1 (MQ=255)
                               gCATAAGTTTAAAAACGTTCCGCAACTGGTATTTACGCACTGTGGTACCCGGTGGTGCTGCGTCAGCAAATTAAGCCCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGtg                                                                                                                                                                  <  1:388511/228‑1 (MQ=255)
                                  gAAGGGTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCACATTAAGCGCATGATCCACGATCTGGAAGAGACCGGCATTCGCAAACTTCTGCAAATTGAACTGGCGGGCCTGCACGACAGCGATTGCCATGG‑CATGACAGATTCCCGTCTGATTTTATTTTACCCCACGCGGGAACTCGaacacaacctgaaaaacggtgcg                                                                                                                                                               >  2:274783/1‑209 (MQ=255)
                                         aCAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCGGACCGCCGTCGCCGGGG‑CATGACCGCTTCCGGCATTATTGGGGTTAACCCGCCCTGGAGACTGGAACAACAGAGGAATAACGTTGTCCCGTgcc                                                                                                                                                        >  1:358676/1‑226 (MQ=255)
                                                tttCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCCATCGCCGTGG‑TATGACCGCCTCCGGCATTATTGTGATTAAACCGCCGGGGAAACTGGAacaaca                                                                                                                                                                               >  1:623889/1‑198 (MQ=255)
                                                        cTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGGGGGAACTGGAACAACAGATAAATAACGGCCTGCCGTGGCTGCACAGCAAACggg                                                                                                                                         >  1:65912/1‑228 (MQ=255)
                                                         tGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACCACAAATGATTAACGTGCTTCCCTGGCTGCACAGCCAACTgtg                                                                                                                                        >  1:244064/1‑226 (MQ=255)
                                                                     aCTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCCCTTCCGGCATGATTGGGATTTACCCGCCGGGGGAACCGGAACAAAAGCTGAATAAATGGGTGCCGGGGCCGCACCACAAACCGGGTCCggcaggacc                                                                                                                            >  2:509243/1‑225 (MQ=255)
                                                                      cTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTTCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGGGATTAACCCGCCGCGGAAACGGGAACAACACAGGAATAAAGTGCTGCCCGGGGTGCGCCGAAAACGGGGTCCGCCGGGccccc                                                                                                                          >  2:647749/1‑228 (MQ=255)
                                                                       tgtgGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTTGGATTAACCCGCCGTGGGAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAACAAACTGGGTCCGGCCGGCACCg                                                                                                                          >  1:552012/1‑228 (MQ=255)
                                                                       tgtgGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGGAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCgg                                                                                                                         >  1:390577/1‑229 (MQ=255)
                                                                       tgtgGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGGGGAAACTGGAACAACAGATGAAAAACGTGCTGCCGTGGCTGGCAAGCCAACCGGTTCCGGCCGGAACCg                                                                                                                          >  1:395805/1‑228 (MQ=255)
                                                                             aTCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGGGATTAACCCGCCGTGGAAAATGGCACAACAGATGAATAAAGTGCTGCCGTGGCCGCCCAGCAAACTGGTTCCGGCAGGCACCAGGCaa                                                                                                                     >  1:222501/1‑226 (MQ=255)
                                                                                 ggtggtGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATTGCCCTGG‑CATGACCTCCTCCGGCATGATTGTCATTACACCGCCCTGGAAACCGGGACAACAGATGAAAAACCTGCTGACGGGGCTGCACAGCCAACCGGTGCCGGCAGGGACCGGGGCCaccacc                                                                                                               >  2:173599/1‑229 (MQ=255)
                                                                                  ttgggtgtggttcaaaaaaTAAAGCGAATAGTACACGATCTGAAAGGGCCCGGTATCCGAAAAATTCTGAAAATGGAACTGGAGGGACTGCTAGAGGGCGACCG‑CGTGGCCATGCCCGCTTCCGGTATGTTGGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCTcc                                                                                                               <  2:247898/212‑1 (MQ=255)
                                                                                   tggtgCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGTCGGGGCTGCCATACAGCGATTGCCGTGG‑CATGACCGCTTCCCGCATGATTTTGATTAACCCGCCCTGAGAACTGGGACAACCAGTAAATAAATACCGGCCCGGGCCGCACAGCCAATTGTGTCCGGCAGGCCCAGGGcccccccccac                                                                                                             >  2:402991/1‑219 (MQ=255)
                                                                                   tggtgCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGGTGAATAACGTGCTGCCCTGGCTGCACAGCAAACTGTTTCCGGCAGCGACCGGGCCAGCCACCgc                                                                                                             >  1:453856/1‑228 (MQ=255)
                                                                                                                 cccACGATCTGGAAGCGCCCGGTAGCCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGCCCGTTTCCGGCATGATTGTGTTTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGcc                                                                              <  1:173599/228‑1 (MQ=255)
                                                                                                                     ctggtggaaaccgcccGGTACTCGAAAATTTTTGAAATTTAATCTGGGGGTACGGCCAGACCGCGTTCGGCGGGG‑AATGACCCTTTCCGGAATGTTTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTg                                                                          <  2:453856/215‑1 (MQ=255)
                                                                                                                       tgttggagggcacggtttctCGCAAATTTTTACAAATTGTCCTGGCGGTACCGCAAGAGACGCCGCCCCGGGC‑CATGCCCTATCCCGGCATGTGTTGGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCt                                                                        <  2:65912/210‑1 (MQ=255)
                                                                                                                           ggAAGCAGCTGGTATGCGCAAAATTCGGCAAATTGAACTGGGGGTACTGCAAGACAGCGCACGCCGTGG‑CATGCCCGCTGCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGc                                                                     <  1:319764/228‑1 (MQ=255)
                                                                                                                           ggAAGCAACCGGTATTCCCAAGATTCTGCAAATTGAACTGGCGGTACTGCCCGACAGCGATCGCCGTGG‑CATAACCGCTTCCGGCCTGATTGTGATTAACCCGCCCTTGAAACTGGAACAACAGATTAATAACGTGCCTCCGTGGCTGCACAGCCAACACGGTCCCGCAAGGCCCCGGCCCCCCCCCGAGAGCTCGA‑TCCAGCCCGAGTAAATGCAGCCCAAGCTGGc                                                                     >  2:92302/1‑228 (MQ=255)
                                                                                                                           gaaggcgcCCGGTATCCGCAAATTTTTGAAAATGTAACTGGGGGTACTGCCAGACAGCGATGGCGTGGG‑CATGACCGCTTCCGGCATAATTGGGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGc                                                                     <  1:402991/224‑1 (MQ=255)
                                                                                                                            gagccgccgggtgtccacAAATTTCGGAAATTGGAACGGGGGGTACTGGCAGACAGCGGCCGCTGTGA‑CAGTCCCGCTCCGGGAAAGATTGGAATAACCCCGCCGTGGAAACGGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCAc                                                                   <  2:215695/212‑1 (MQ=255)
                                                                                                                                  aCCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTATTGCCTGACAGCGCGCGCCTTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTgg                                                                                                                                                                                      <  1:330744/109‑1 (MQ=255)
                                                                                                                                  aCCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CCTGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGCATAAAGTGCTGCCGTGGCCGCACAGCAAACTGGGTCCGTCAGGCCCCGGGGACCCCCCCGCACACTGGA‑TCCGTCCGGGGACATTTCAGCCATTGCGGGccccaaat                                                              >  2:86835/1‑224 (MQ=255)
                                                                                                                                  aCCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAACGATCGGCGTGG‑TAATACCGCCTCCGGCGTTATGGTGATTAAACCGCCGTGGAAACAGGAAACACAGACGAATAAAGGTGTTCCGGCGCCGGCCAGCAAAAAGGGTGCAGCAGGCACCGGGAACGCCACCGTGAGCCGGACTCG‑GCCCGAGAAATTGCGGCCAGTGcgcgcgcccatt                                                              >  2:163731/1‑217 (MQ=255)
                                                                                                                                                     tGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCCCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACCGCAAACCGGTTCCGGCAGGCACCGGGCACg                                                                                                                    >  1:171441/1‑156 (MQ=255)
                                                                                                                                                          tttacactggggtggcggCCAGACAGCGACCGCGGAGG‑AATGACCGCTTCCGGAATGATTGGGATAAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTaa                                     <  2:552011/212‑1 (MQ=255)
                                                                                                                                                               cgtggcggatcgccaaGACGGCTATCGCCGGTG‑CAAGACGTTTGCGGGCATGTTGGTGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAAc                                <  1:183013/214‑1 (MQ=255)
                                                                                                                                                                      gttcTGCCAGACAGCGGTCGCCGTGG‑AATGACCGTTTCCGGCATGGTTGGGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                         <  1:92302/226‑1 (MQ=255)
                                                                                                                                                                      gacccggcgagagcggATCCCCGGGG‑CATGCCCGCTTCGGGATTGTTTGGGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                         <  1:94306/214‑1 (MQ=255)
                                                                                                                                                                      gTACCGCCAGACAGGGATCGCCGTGC‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                         <  1:521346/229‑1 (MQ=255)
                                                                                                                                                                           tcaggACAGCGCTCGCGGGGA‑CAGGCCCGCTCCCGGCATGTTTGTTATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACt                    <  1:411078/225‑1 (MQ=255)
                                                                                                                                                                                 cAGCGATGGCTGGGG‑CATGACCGCTCCGGGCATGATTGTGAATACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTc               <  1:86835/228‑1 (MQ=255)
                                                                                                                                                                                 cAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCCCCGTAAGCTGGA‑TCGTGCCGGAGCAATTGCAGCCCTTGCTTGCACCCATTACGTCCCCGCCTCCAATACCGGAAAACAACGATAAAAACACttgtaca              >  2:487318/1‑224 (MQ=255)
                                                                                                                                                                                  aGCGATCGCCGGGG‑CATGCCCGCTTCCGGCATGATTGTGATTACCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCa              <  1:509243/228‑1 (MQ=255)
                                                                                                                                                                                  aGCGACCGCCGTGG‑CATGACCGTTTCCGGCATGATTGGGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCa              <  1:16915/228‑1 (MQ=255)
                                                                                                                                                                                           cgggggcgagccccagcgcGGGATGAGGGTGAATAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTGCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaaa      <  2:244064/210‑1 (MQ=255)
                                                                                                                                                                                                gcagggccGCTTCGGGCTTGATTGTGATCAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGGGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACt  <  1:308826/221‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                           
AACTGACTATCAAGCGGTGGTCAGCGGGATAGCAGAAGGTTACAAACGTTTCGCCACTGGTATTTACGCACTGTGGTATCCGGTGGTGCTGCGTCAGCAAATTAAGCGCATGATCCACGATCTGGAAGCGACCGGTATTCGCAAAATTCTGCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGG‑CATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGA‑TCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACT  >  NC_000913/3645893‑3646310

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