Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 899,657 G→A 6.6% 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.6% 296.6 / 8.1 91intergenic (+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 (46/39);  new base A (2/4);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

AGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCTTT  >  NC_000913/899514‑899805
                                                                                                                                               |                                                                                                                                                    
aGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCgg                                                                                                                                                 >  2:84573/1‑149 (MQ=255)
     tAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCa                                                                                                                                            <  2:461971/149‑1 (MQ=255)
     tAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCa                                                                                                                                            <  3:66499/149‑1 (MQ=255)
     tAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCtttt                                                                                                                                            <  2:175781/149‑3 (MQ=255)
     tAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCtttt                                                                                                                                            <  3:638808/149‑3 (MQ=255)
       aCGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGa                                                                                                                                          >  4:626457/1‑149 (MQ=255)
       aCGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGa                                                                                                                                          <  1:622109/149‑1 (MQ=255)
           ccAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGc                                                                                                                                      >  2:27618/1‑149 (MQ=255)
           ccAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGc                                                                                                                                      <  3:708044/149‑1 (MQ=255)
            cAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     <  1:293465/149‑1 (MQ=255)
            cAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCt                                                                                                                                     <  3:90017/149‑1 (MQ=255)
               aCTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGAc                                                                                                                                  <  4:238907/149‑1 (MQ=255)
                      cGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTgc                                                                                                                           <  1:826088/149‑1 (MQ=255)
                        aaaGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTgcgc                                                                                                                         >  4:176427/1‑149 (MQ=255)
                        aaaGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTgcgc                                                                                                                         <  1:634530/149‑1 (MQ=255)
                             tatCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgt                                                                                                                    >  4:280519/1‑149 (MQ=255)
                               tCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGttt                                                                                                                  <  1:681896/149‑1 (MQ=255)
                                ccGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGc                                                                                                                                                   <  2:555483/115‑1 (MQ=255)
                                ccGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGc                                                                                                                                                   >  1:555483/1‑115 (MQ=255)
                                ccGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACaaa                                                                                                                               <  1:597314/135‑1 (MQ=255)
                                ccGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACaaa                                                                                                                               >  2:597314/1‑135 (MQ=255)
                                ccGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTa                                                                                                                 <  4:270516/149‑1 (MQ=255)
                                  gcgAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATg                                                                                                               <  2:322949/149‑1 (MQ=255)
                                  gcgAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATg                                                                                                               >  1:256272/1‑149 (MQ=255)
                                            tGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgt                                                                                                                    <  3:688540/134‑1 (MQ=255)
                                            tGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgt                                                                                                                    >  4:688540/1‑134 (MQ=255)
                                             gggTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgcggc                                                                                                    <  4:391234/149‑1 (MQ=255)
                                               gTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGt                                                                                                  >  2:864611/1‑149 (MQ=255)
                                                tttCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTa                                                                                                 >  3:529508/1‑149 (MQ=255)
                                                     tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGc                                                                                            >  3:597283/1‑149 (MQ=255)
                                                     tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGc                                                                                            >  3:44953/1‑149 (MQ=255)
                                                     tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGc                                                                                            >  2:240382/1‑149 (MQ=255)
                                                     tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCtttt                                                                                                                                            >  1:93799/1‑99 (MQ=255)
                                                     tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCtttt                                                                                                                                            <  2:93799/101‑3 (MQ=255)
                                                      aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           >  3:210148/1‑149 (MQ=255)
                                                      aTTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAACAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTCTGCCGGATGCGGCGTAAACGcc                                                                                           >  1:633417/1‑149 (MQ=255)
                                                       ttGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                          >  4:719382/1‑149 (MQ=255)
                                                        tGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCtt                                                                                         <  4:791093/149‑1 (MQ=255)
                                                        tGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCtt                                                                                         >  1:195280/1‑149 (MQ=255)
                                                          aaCGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       <  3:280519/149‑1 (MQ=255)
                                                          aaCGGGTGCAGCTCTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       >  2:77967/1‑149 (MQ=255)
                                                           aCGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                                                      <  1:848431/149‑1 (MQ=255)
                                                               gTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                  <  1:499025/149‑1 (MQ=255)
                                                               gTGCAGCGTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCATATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                  <  3:549868/149‑1 (MQ=255)
                                                                tGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                 >  1:139759/1‑149 (MQ=255)
                                                                 gCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                >  2:87737/1‑149 (MQ=255)
                                                                 gCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                <  3:429285/149‑1 (MQ=255)
                                                                  cAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                               >  2:808298/1‑149 (MQ=255)
                                                                  cAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                               >  2:390600/1‑149 (MQ=255)
                                                                  cAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                               <  2:194461/149‑1 (MQ=255)
                                                                   aGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                              >  2:477935/1‑149 (MQ=255)
                                                                    gCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                             >  1:726272/1‑149 (MQ=255)
                                                                       tttCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                                          <  3:358254/149‑1 (MQ=255)
                                                                           gATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGt                                                                      <  1:621700/149‑1 (MQ=255)
                                                                           gATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGt                                                                      <  4:529508/149‑1 (MQ=255)
                                                                            atatGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTg                                                                     >  3:692001/1‑149 (MQ=255)
                                                                                     cGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            >  2:618037/1‑149 (MQ=255)
                                                                                     cGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            <  3:684176/149‑1 (MQ=255)
                                                                                     cGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                            <  4:317753/149‑1 (MQ=255)
                                                                                           ccGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACAtt                                                      <  1:164193/149‑1 (MQ=255)
                                                                                             gcgcACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGc                                                    >  3:22039/1‑149 (MQ=255)
                                                                                               gcACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAt                                                  <  3:730492/149‑1 (MQ=255)
                                                                                                   tATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGc                                              >  3:223910/1‑149 (MQ=255)
                                                                                                    aTGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCTTTAt                                                                                       <  2:237811/107‑1 (MQ=255)
                                                                                                    aTGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                       >  1:237811/1‑107 (MQ=255)
                                                                                                      gAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATgg                                                <  4:415242/144‑1 (MQ=255)
                                                                                                      gAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATgg                                                >  3:415242/1‑144 (MQ=255)
                                                                                                      gAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAt                                           <  1:84573/149‑1 (MQ=255)
                                                                                                            cTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGc                                     >  3:382224/1‑149 (MQ=255)
                                                                                                              gACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCt                                   <  4:772458/149‑1 (MQ=255)
                                                                                                              gACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCt                                   >  1:328478/1‑149 (MQ=255)
                                                                                                                   tgctgGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATaa                              >  4:50192/1‑149 (MQ=255)
                                                                                                                   tgctgGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATaa                              >  2:550575/1‑149 (MQ=255)
                                                                                                                   tgctgGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATaa                              <  2:195280/149‑1 (MQ=255)
                                                                                                                   tgctgGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATaa                              >  3:825751/1‑149 (MQ=255)
                                                                                                                          gAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAgc                       <  2:256272/149‑1 (MQ=255)
                                                                                                                           aaGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAgcg                      >  2:96921/1‑149 (MQ=255)
                                                                                                                           aaGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAgcg                      >  1:692109/1‑149 (MQ=255)
                                                                                                                              cAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCAt                   >  2:501961/1‑149 (MQ=255)
                                                                                                                               aaTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATc                  >  1:591215/1‑149 (MQ=255)
                                                                                                                                aTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCa                 >  1:264838/1‑149 (MQ=255)
                                                                                                                                aTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCa                 >  3:476515/1‑149 (MQ=255)
                                                                                                                                 tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           <  1:868826/74‑1 (MQ=255)
                                                                                                                                 tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                           >  2:868826/1‑74 (MQ=255)
                                                                                                                                 tAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAg                >  4:250215/1‑149 (MQ=255)
                                                                                                                                  tctgaaGCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAg                >  1:616436/5‑148 (MQ=255)
                                                                                                                                  tctgaaGCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAg                <  2:616436/144‑1 (MQ=255)
                                                                                                                                    aaaaGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAAGACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa             >  1:260068/1‑149 (MQ=255)
                                                                                                                                     aaaGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa            >  4:792857/1‑149 (MQ=255)
                                                                                                                                      aaGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt           <  3:50192/149‑1 (MQ=255)
                                                                                                                                            cATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAcc     <  1:73896/149‑1 (MQ=255)
                                                                                                                                              tgtgCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCtt   <  1:748268/149‑1 (MQ=255)
                                                                                                                                               gtgCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCttt  <  1:33899/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTACCTTT  >  NC_000913/899514‑899805

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