Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 899,657 G→A 6.1% intergenic (+12/+187) rlmC → / ← artJ 23S rRNA m(5)U747 methyltransferase, SAM‑dependent/arginine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913899,6570GA 6.1% 221.9 / 4.0 66intergenic (+12/+187)rlmC/artJ23S rRNA m(5)U747 methyltransferase, SAM‑dependent/arginine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCTTT  >  NC_000913/899515‑899805
                                                                                                                                              |                                                                                                                                                    
gCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGc                                                                                                                                                <  2:1787559/149‑1 (MQ=255)
 cTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCt                                                                                                                                               >  2:1392710/1‑149 (MQ=255)
          ccAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCAtgt                                                                                                                                                     >  1:1752638/1‑134 (MQ=255)
          ccAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCAtgt                                                                                                                                                     <  2:1752638/134‑1 (MQ=255)
           cAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     <  1:1392710/149‑1 (MQ=255)
                 aTGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACaaa                                                                                                                               <  3:1380087/149‑1 (MQ=255)
                    gCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTg                                                                                                                            >  4:860654/1‑149 (MQ=255)
                      gAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGGAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTgcg                                                                                                                          >  3:141850/1‑149 (MQ=255)
                       aaaGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCa                                                                                                                                            <  3:104513/130‑1 (MQ=255)
                       aaaGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCa                                                                                                                                            >  4:104513/1‑130 (MQ=255)
                       aaaGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     >  4:131577/1‑137 (MQ=255)
                       aaaGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     <  3:131577/137‑1 (MQ=255)
                            tatCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgt                                                                                                                    <  2:25560/149‑1 (MQ=255)
                             atCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgtt                                                                                                                   >  4:699014/1‑149 (MQ=255)
                                         ccTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgc                                                                                                       <  3:860654/149‑1 (MQ=255)
                                            gggTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgcggc                                                                                                    >  4:1492328/1‑149 (MQ=255)
                                             ggTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgcggcg                                                                                                   >  2:1202015/1‑149 (MQ=255)
                                                  cGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAAc                                                                                              <  3:699014/149‑1 (MQ=255)
                                                    tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGc                                                                                            <  3:1492328/149‑1 (MQ=255)
                                                     aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           >  3:21261/1‑149 (MQ=255)
                                                     aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           >  3:1458958/1‑149 (MQ=255)
                                                     aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           <  4:1458958/149‑1 (MQ=255)
                                                     aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           <  4:21261/149‑1 (MQ=255)
                                                     aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           <  1:1202015/149‑1 (MQ=255)
                                                      ttGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCTCGCACTATGCAGTGCTGACGCTGCTGGTGGAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCCGATGCGGCGTAAACGACt                                                                                          <  4:141850/149‑1 (MQ=255)
                                                      ttGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                          >  4:939906/1‑149 (MQ=255)
                                                        gAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGa                                                                                                                                          <  4:1350614/99‑1 (MQ=255)
                                                        gAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGa                                                                                                                                          >  3:1350614/1‑99 (MQ=255)
                                                         aaCGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     <  2:411417/103‑1 (MQ=255)
                                                         aaCGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     >  1:411417/1‑103 (MQ=255)
                                                          aCGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                                                      <  3:939906/149‑1 (MQ=255)
                                                                 cAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCttt                                                                                                                                             <  1:1782840/87‑2 (MQ=255)
                                                                 cAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCttt                                                                                                                                             >  2:1782840/1‑86 (MQ=255)
                                                                           atatGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTg                                                                     >  4:1095892/1‑149 (MQ=255)
                                                                            tatGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGc                                                                    >  4:293310/1‑149 (MQ=255)
                                                                               gTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCaaa                                                                 >  3:1608985/1‑149 (MQ=255)
                                                                                    cGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            <  3:1095892/149‑1 (MQ=255)
                                                                                         aCCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTgc                                                                                                                           <  4:1836613/81‑1 (MQ=255)
                                                                                         aCCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTgc                                                                                                                           >  3:1836613/1‑81 (MQ=255)
                                                                                                 cTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       <  2:789448/109‑1 (MQ=255)
                                                                                                 cTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       >  1:789448/1‑109 (MQ=255)
                                                                                                          gCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAgg                                      <  3:293310/149‑1 (MQ=255)
                                                                                                             gACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCt                                   <  4:1608985/149‑1 (MQ=255)
                                                                                                              aCGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgcggc                                                                                                    <  3:415258/83‑1 (MQ=255)
                                                                                                              aCGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgcggc                                                                                                    >  4:415258/1‑83 (MQ=255)
                                                                                                                       gTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                          >  1:827296/1‑84 (MQ=255)
                                                                                                                       gTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                          <  2:827296/84‑1 (MQ=255)
                                                                                                                           aGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            >  1:1066475/1‑110 (MQ=255)
                                                                                                                           aGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            <  2:1066475/110‑1 (MQ=255)
                                                                                                                             cAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGc                            >  1:1507304/1‑140 (MQ=255)
                                                                                                                             cAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGc                            <  2:1507304/140‑1 (MQ=255)
                                                                                                                              aaTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                                                      <  1:90822/81‑1 (MQ=255)
                                                                                                                              aaTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                                                      >  2:90822/1‑81 (MQ=255)
                                                                                                                              aaTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            <  2:358839/107‑1 (MQ=255)
                                                                                                                              aaTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            >  1:358839/1‑107 (MQ=255)
                                                                                                                                tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       >  3:1362555/1‑78 (MQ=255)
                                                                                                                                tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       <  4:1362555/78‑1 (MQ=255)
                                                                                                                                tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATgg                                                >  2:668542/1‑117 (MQ=255)
                                                                                                                                tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATgg                                                <  1:668542/117‑1 (MQ=255)
                                                                                                                                    gaaGCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGAt                                >  3:705723/2‑129 (MQ=255)
                                                                                                                                    gaaGCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGAt                                <  4:705723/128‑1 (MQ=255)
                                                                                                                                     aaGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGa                                 >  2:846416/1‑127 (MQ=255)
                                                                                                                                     aaGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGa                                 <  1:846416/127‑1 (MQ=255)
                                                                                                                                       gCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATg                                                                                                        <  3:744831/54‑1 (MQ=255)
                                                                                                                                       gCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATg                                                                                                        >  4:744831/1‑54 (MQ=255)
                                                                                                                                          gCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAc      >  3:116990/1‑149 (MQ=255)
                                                                                                                                           cATGTGCTGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGAGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAcc     >  4:241109/1‑149 (MQ=255)
                                                                                                                                           cATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAcc     <  3:241109/149‑1 (MQ=255)
                                                                                                                                             tgtgCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCtt   >  1:1806476/1‑149 (MQ=255)
                                                                                                                                              gtgCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCttt  <  4:116990/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCTTT  >  NC_000913/899515‑899805

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