Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,235,823 G→A A18V (GCG→GTG)  rlmG ← 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,235,8230GA97.8% 333.6 / ‑3.8 92A18V (GCG→GTG) rlmG23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (45/45);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTTATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCGGTTATGAGCAATCT  >  NC_000913/3235678‑3235970
                                                                                                                                                 |                                                                                                                                                   
cTTATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCgg                                                                                                                                                 >  2:292623/1‑150 (MQ=255)
 ttATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCgg                                                                                                                                                 >  1:98193/1‑149 (MQ=255)
  tATGTTCCGCGAGTGCGCAGCTTACCGCACCAAAGGCAACATTCAGAATCAACACCGGGCCGCGGATTTCTGTGTCATCCATTTGTTGAAGAAAATATTCATCCTCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATTCGTCACCg                                                                                                                                                  <  1:983685/147‑1 (MQ=255)
  tATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGaa                                                                                                                                               >  2:941923/1‑150 (MQ=255)
  tATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGaa                                                                                                                                               <  1:295686/150‑1 (MQ=255)
  tATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGaa                                                                                                                                               <  1:699524/150‑1 (MQ=255)
  tATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGaa                                                                                                                                               >  2:509155/1‑150 (MQ=255)
    tGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAAc                                                                                                                                            >  2:450158/1‑151 (MQ=255)
      ttCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACg                                                                                                                                           <  1:278044/150‑1 (MQ=255)
      ttCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACg                                                                                                                                           <  1:776566/150‑1 (MQ=255)
       tCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGt                                                                                                                                          <  1:731329/150‑1 (MQ=255)
        ccGCCAGTGCGCAGCTTAACGCCCCAAAGGCATCATTCAGGAACAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCGCCGGAAAACGTTg                                                                                                                                        <  2:1400929/151‑1 (MQ=255)
        ccGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGtt                                                                                                                                         <  2:379970/150‑1 (MQ=255)
             aGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCCTCCGTCACCGGAAAACGTTGTAgt                                                                                                                                    >  1:264287/1‑150 (MQ=255)
             aGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAgtg                                                                                                                                   <  2:831483/151‑1 (MQ=255)
               tGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTc                                                                                                                                 >  1:1159726/1‑151 (MQ=255)
                gcgcAGCTTAACGCCCCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGCCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCa                                                                                                                                <  1:1005157/151‑1 (MQ=255)
                   cAGCTTATCGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCa                                                                                                                                <  1:73106/148‑1 (MQ=255)
                   cAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTg                                                                                                                             >  1:1231599/1‑151 (MQ=255)
                   cAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTg                                                                                                                             >  2:553981/1‑151 (MQ=255)
                        tAACGCACCAAGGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCAGCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAAc                                                                                                                          <  1:1221777/149‑1 (MQ=255)
                        tAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAAc                                                                                                                          <  1:1152995/149‑1 (MQ=255)
                             cACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAAcc                                                                                                                   >  2:807024/1‑151 (MQ=255)
                                cAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGtt                                                                                                                <  2:369971/151‑1 (MQ=255)
                                cAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGtt                                                                                                                >  2:890125/1‑151 (MQ=255)
                                cAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGtt                                                                                                                >  2:494591/1‑151 (MQ=255)
                                 aaaGGCATCATTCAGGATCAACCCCGGGCCGCGGATTTCTGTGTCGTCCACCTTTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCg                                                                                                                  <  1:1173191/148‑1 (MQ=255)
                                      catcatTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTa                                                                                                           >  2:167936/1‑150 (MQ=255)
                                      catcatTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTa                                                                                                           >  2:1143912/1‑150 (MQ=255)
                                         catTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTg                                                                                                       <  1:1180832/151‑1 (MQ=255)
                                          atTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTg                                                                                                       <  2:1314988/150‑1 (MQ=255)
                                            tCAGGATCAACACCGGGCCGCGGATTTATGTGCCGTCCAACTGTTGCAGCAAATATCCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTAACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGc                                                                                                     >  1:1206993/1‑150 (MQ=255)
                                             cAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTc                                                                                                   >  2:456448/1‑151 (MQ=255)
                                              aGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTc                                                                                                   >  2:482299/1‑150 (MQ=255)
                                              aGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTc                                                                                                   >  2:817528/1‑150 (MQ=255)
                                               ggATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCAt                                                                                                 <  2:669204/151‑1 (MQ=255)
                                                  tCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCg                                                                                               <  1:348179/150‑1 (MQ=255)
                                                   cAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGc                                                                                             >  2:876772/1‑151 (MQ=255)
                                                       acCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctc                                                                                         >  2:914890/1‑151 (MQ=255)
                                                         cGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctcct                                                                                       >  2:475371/1‑151 (MQ=255)
                                                          gggCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTg                                                                                      <  2:38400/151‑1 (MQ=255)
                                                          gggCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAAGATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTg                                                                                      <  2:708336/151‑1 (MQ=255)
                                                             ccGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAAt                                                                                   <  1:1415526/151‑1 (MQ=255)
                                                                  gATGTCTGTGTCGTCGAACTGTTGAAGAAAATATTCAACCGCCGCTCCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTAAACGGAAGCCGTTGTCTAAGTGGCTCATCGGCCCTCCCGAATGATaa                                                                              <  1:1326600/151‑1 (MQ=255)
                                                                     ttCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATaaa                                                                             >  2:841019/1‑149 (MQ=255)
                                                                     ttCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAAt                                                                           >  2:130842/1‑151 (MQ=255)
                                                                        tgtgTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTc                                                                        >  1:1326149/1‑151 (MQ=255)
                                                                         gtgtCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTc                                                                        <  2:1348188/150‑1 (MQ=255)
                                                                            tcgtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGAc                                                                     <  1:438267/150‑1 (MQ=255)
                                                                            tcgccaaaCGGTTGCAGCAAATATTCATCCGCCGCTGCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGa                                                                      <  1:929631/143‑1 (MQ=255)
                                                                               tcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGa                                                                      <  1:742200/146‑1 (MQ=255)
                                                                               tcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTg                                                                 >  1:474551/1‑151 (MQ=255)
                                                                                ccAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGc                                                                >  2:1075234/1‑151 (MQ=255)
                                                                                 cAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGcc                                                               >  2:445798/1‑151 (MQ=255)
                                                                                 cAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGcc                                                               >  1:1064669/1‑151 (MQ=255)
                                                                                     tGTTGCAGCAAATATTCATCCGCCGCTTCCAACGCGGGTAGTGGGGTAACGTCGTCCGTCACCGGGAAAAGATGTAGTGTCAGTGAACGGAAACCGCTGTCTAAGAGGTTACTCGGTCCTCCTGAATTATAAAATTTCTACATTGCCTCTg                                                           >  1:363319/1‑151 (MQ=255)
                                                                                        tGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGaa                                                         >  2:165174/1‑150 (MQ=255)
                                                                                        tGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGaa                                                         <  1:1388523/150‑1 (MQ=255)
                                                                                         gcagcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGaaaa                                                       <  2:594076/151‑1 (MQ=255)
                                                                                           agcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAg                                                      <  1:1388435/150‑1 (MQ=255)
                                                                                             caAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAggg                                                    >  1:331383/1‑150 (MQ=255)
                                                                                              aaaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGt                                                   <  1:423151/150‑1 (MQ=255)
                                                                                              aaaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTg                                                  <  1:241971/151‑1 (MQ=255)
                                                                                               aaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGc                                                 <  1:553933/151‑1 (MQ=255)
                                                                                               aaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGc                                                 <  2:1182759/151‑1 (MQ=255)
                                                                                               aaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGc                                                 >  2:1338642/1‑151 (MQ=255)
                                                                                                  atTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGt                                               <  1:471227/150‑1 (MQ=255)
                                                                                                      aTCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGt                                           >  2:206228/1‑150 (MQ=255)
                                                                                                       tCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGtat                                         >  2:130904/1‑151 (MQ=255)
                                                                                                        ccgccgCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGtata                                        <  1:937617/151‑1 (MQ=255)
                                                                                                           ccgcTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATAccc                                     <  1:208738/151‑1 (MQ=255)
                                                                                                              cTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCttt                                  >  2:1064589/1‑151 (MQ=255)
                                                                                                                tCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCttttt                                >  2:1061888/1‑151 (MQ=255)
                                                                                                                tCCCACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCttt                                  <  2:1193185/149‑1 (MQ=255)
                                                                                                                 cccACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCttttt                                <  1:1407150/150‑1 (MQ=255)
                                                                                                                 cccACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCttttt                                >  2:1003447/1‑150 (MQ=255)
                                                                                                                 cccACGCCTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTCGCTCATCGGCCCTCCTGAATGATAAACTTTCGACATTGCCCCTGAAAAGGGCGCGTGAGTATACCCTTTTtc                               >  1:52944/1‑151 (MQ=255)
                                                                                                                      gCCTGTAGCGGGTTAACATCATCCGTCACCGGAAAACGTTGTAGAGCCAGTGACCGGAAACCGTCGCCTACCTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTCCGTGCGTATACCCTTTTTCTCttt                          >  1:319570/1‑151 (MQ=255)
                                                                                                                        cTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCttttt                        >  2:1147014/1‑151 (MQ=255)
                                                                                                                        cTGTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCttttt                        >  2:1108355/1‑151 (MQ=255)
                                                                                                                          gTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTc                       <  2:552370/150‑1 (MQ=255)
                                                                                                                          gTAGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTc                       <  2:32544/150‑1 (MQ=255)
                                                                                                                           tAGCGGGTTAACGTCATCCGTCACGTGAAAACGTTGTAATGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTGCTAAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATTCCCTGTTTCTCttgtt                        <  1:1071832/148‑1 (MQ=255)
                                                                                                                            aGCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGt                     <  1:450138/150‑1 (MQ=255)
                                                                                                                             gCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTgg                   >  1:601885/1‑151 (MQ=255)
                                                                                                                             gCGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTgg                   <  2:288840/151‑1 (MQ=255)
                                                                                                                              cGGGTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTggg                  <  2:724167/151‑1 (MQ=255)
                                                                                                                                 gTTAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCgg               <  2:528400/151‑1 (MQ=255)
                                                                                                                                   tAACGTCATCCGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCGGtt             >  2:444818/1‑151 (MQ=255)
                                                                                                                                             cGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCGGTTATGAGCAATc   <  1:607141/151‑1 (MQ=255)
                                                                                                                                             cGTCACCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCGGTTATGAGCAATc   >  2:116755/1‑151 (MQ=255)
                                                                                                                                              gTCACCGGAAAACGTTGTAGTGTCAGTGAACGGGAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCGGTTATGAGCAATct  >  1:376023/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTTATGTTCCGCCAGTGCGCAGCTTAACGCACCAAAGGCATCATTCAGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGACATTGCCCCTGAAAAGGGTGCGTGAGTATACCCTTTTTCTCTTTTTCGTGGGCGGTTATGAGCAATCT  >  NC_000913/3235678‑3235970

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