Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,643,721 A→G 5.9% V188A (GTC→GCC)  rlmN ← dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,643,7210AG 5.9% 206.7 / 2.8 68V188A (GTC→GCC) rlmNdual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base G (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

GGCGCGTGCAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGACACCCGCAGGTTGCGGTT  >  NC_000913/2643574‑2643866
                                                                                                                                                   |                                                                                                                                                 
ggCGCGTGCAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAacgac                                                                                                                                                  <  5:54193/149‑1 (MQ=255)
   gcgTGCAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACgtt                                                                                                                                               <  8:95965/149‑1 (MQ=255)
    cgTGCAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGTAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACgttg                                                                                                                                              <  5:426520/149‑1 (MQ=255)
        cAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTc                                                                                                                                           <  4:302351/148‑1 (MQ=255)
        cAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCa                                                                                                                                          <  2:266360/149‑1 (MQ=255)
           ggAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGt                                                                                                                                       >  3:205567/1‑149 (MQ=255)
                       tGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCt                                                                                                                           >  5:306303/1‑149 (MQ=255)
                          aaCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgc                                                                                                                        <  6:286332/149‑1 (MQ=255)
                          aaCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgc                                                                                                                        <  3:356040/149‑1 (MQ=255)
                          aaCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgc                                                                                                                        >  7:83965/1‑149 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGGCGTTGTTCAGGTTGagcagc                                                                                                                               >  3:331145/1‑141 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGGCGTTGTTCAGGTTGagcagc                                                                                                                               <  4:325428/141‑1 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGGCGTTGTTCAGGTTGagcagc                                                                                                                               <  4:331145/141‑1 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGGCGTTGTTCAGGTTGagcagc                                                                                                                               >  3:325428/1‑141 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACgtt                                                                                                                                               >  5:96859/1‑125 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACgtt                                                                                                                                               <  6:96859/125‑1 (MQ=255)
                           aCGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACGACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCCTGATTTCCATCGCCGGAACGAAGTTGTTCAGGTTGAGCAGCGGCTctca                                                                                                                       >  6:310202/1‑146 (MQ=255)
                              tCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgcccat                                                                                                                    <  3:404551/149‑1 (MQ=255)
                              tCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgcccat                                                                                                                    <  1:234320/149‑1 (MQ=255)
                                gATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgcccatgc                                                                                                                  >  7:245594/1‑149 (MQ=255)
                                          cccAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgccc                                                                                                                      <  2:223786/135‑1 (MQ=255)
                                          cccAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgccc                                                                                                                      >  1:223786/1‑135 (MQ=255)
                                               tttATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGt                                                                                                   <  8:3096/149‑1 (MQ=255)
                                                  aTCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTa                                                                                                 <  6:306303/148‑1 (MQ=255)
                                                   tCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGt                                                                                               >  8:454220/1‑149 (MQ=255)
                                                       gcgcCGGAACTACACCGTCAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCCAGCAGGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATc                                                                                           <  5:310202/149‑1 (MQ=255)
                                                         gcCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTa                                                                                                 >  6:207877/1‑141 (MQ=255)
                                                         gcCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTa                                                                                                 <  5:207877/141‑1 (MQ=255)
                                                           cGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGAc                                                                                       >  2:291208/1‑149 (MQ=255)
                                                                   cacCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGAc                                                                                 <  7:454220/147‑1 (MQ=255)
                                                                    acCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGt                                                                                               <  1:326137/132‑1 (MQ=255)
                                                                    acCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGt                                                                                               >  2:326137/1‑132 (MQ=255)
                                                                      cGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGag                                                                                                                                   <  8:345811/94‑1 (MQ=255)
                                                                      cGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGag                                                                                                                                   >  7:345811/1‑94 (MQ=255)
                                                                      cGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGt                                                                             >  3:46394/1‑148 (MQ=255)
                                                                      cGGAAGTGGAAAGCGAGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCAATCGACGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCAATGACCATCATCACTACGTTAGTGATCGTACTCTGAACGGt                                                                             >  2:306842/1‑148 (MQ=255)
                                                                        gAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCcatcat                                                                                                           <  5:239067/116‑1 (MQ=255)
                                                                        gAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCgg                                                                              >  5:155396/1‑145 (MQ=255)
                                                                        gAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCgg                                                                              <  6:155396/145‑1 (MQ=255)
                                                                        gAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGACCAGCGGCTCGCCCATGCCcatcat                                                                                                           >  6:239067/1‑116 (MQ=255)
                                                                             ggAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtg                                                                     >  4:171578/1‑149 (MQ=255)
                                                                                 aGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgctgc                                                                 >  5:184266/1‑149 (MQ=255)
                                                                                           gACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCcgacga                                                         <  5:156927/147‑1 (MQ=255)
                                                                                           gACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCcgacga                                                         >  6:156927/1‑147 (MQ=255)
                                                                                           gACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGAtt                                                       <  3:171578/149‑1 (MQ=255)
                                                                                            aCGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGAttt                                                      <  8:245594/149‑1 (MQ=255)
                                                                                              gTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGTACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTc                                                    >  6:338896/1‑149 (MQ=255)
                                                                                              gTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTc                                                    <  4:155690/149‑1 (MQ=255)
                                                                                               tttAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTa                                                                                                 >  1:189350/1‑103 (MQ=255)
                                                                                               tttAGACAGGCCAAAACAGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTa                                                                                                 <  2:189350/103‑1 (MQ=255)
                                                                                                 tAGACAGGCCAAACCCGACCTCATCAAGCATGATTTCCATCGCCGGAACGACGGTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCatcatc                                                                                                          <  7:327392/92‑1 (MQ=255)
                                                                                                   gACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGAAGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgccg                                                >  5:132547/1‑148 (MQ=255)
                                                                                                        gCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGcc                                          >  2:235363/1‑149 (MQ=255)
                                                                                                          cAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCcatcat                                                                                                           >  2:447789/1‑82 (MQ=255)
                                                                                                          cAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCcatcat                                                                                                           <  1:447789/82‑1 (MQ=255)
                                                                                                          cAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCg                                                                                          >  7:63388/1‑99 (MQ=255)
                                                                                                          cAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCg                                                                                          <  8:63388/99‑1 (MQ=255)
                                                                                                             aaCCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACAcc                                     >  6:154632/1‑149 (MQ=255)
                                                                                                               ccGAAATCATCAAGCATGATTTCCATCGCCGGTACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTg                                   <  5:338896/149‑1 (MQ=255)
                                                                                                                cGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCCGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGCCCGGTGACTTTTGCTGCACCGACGATTTTCGCCGCACGCCACACCTgg                                  >  7:296449/1‑149 (MQ=255)
                                                                                                                 gAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGc                                 <  8:83965/149‑1 (MQ=255)
                                                                                                                       atcaAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGATCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGATGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATaa                           <  6:132547/149‑1 (MQ=255)
                                                                                                                            gCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGTCGCACGCCACACCTGGCCGATAATTTcc                      <  4:46394/149‑1 (MQ=255)
                                                                                                                               tGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGa                    >  6:175334/1‑148 (MQ=255)
                                                                                                                                       aTCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGACACCCGCAg           >  5:439618/1‑149 (MQ=255)
                                                                                                                                           ccGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGACACCCGCAGGTTg       <  6:184266/149‑1 (MQ=255)
                                                                                                                                            cGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgctg                                                                  <  5:339427/89‑1 (MQ=255)
                                                                                                                                            cGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgctg                                                                  >  6:339427/1‑89 (MQ=255)
                                                                                                                                              gAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGACACCCGCAGGTTGCgg    <  4:205567/149‑1 (MQ=255)
                                                                                                                                                acgacgTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGACACCCGCAGGTTGCGGtt  <  1:150124/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGCGTGCAGGGAAATTGCCAGTGCAACGTCGATCATATCGCCCAGTTTATCCAGCGCCGGAACTACACCGGAAGTGGAAAGCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGACACCCGCAGGTTGCGGTT  >  NC_000913/2643574‑2643866

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