Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,365 C→A 5.9% R31S (CGT→AGT)  rlmB → 23S rRNA 2'‑O‑ribose G2251 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,3650CA 5.9% 270.4 / 4.4 74R31S (CGT→AGT) rlmB23S rRNA 2'‑O‑ribose G2251 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base A (3/2);  total (40/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 = 9.63e-01

TTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGC  >  NC_000913/4409108‑4409604
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAAAGTTTTAAGGAAGTCTTTATTTTAAAAGGGAGTGGAAATAAACgg                                                                                                                                                                                                                                    >  2:793536/1‑271 (MQ=255)
                                  gcTTTGTTCATGCCGTATGCGGCGTGAACGCCTTATCCTGACTACATAATCACGTAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAATGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                                                                  <  1:1022748/272‑1 (MQ=255)
                                  gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                                                                  <  1:896317/272‑1 (MQ=255)
                                  gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                                                                  <  1:793537/272‑1 (MQ=255)
                                                                                      cGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATATCGACGAGAAAAGCGACgg                                                                                                                                              >  1:304188/1‑272 (MQ=255)
                                                                                                                aTGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAAc                                                                                                                                                                           >  2:820012/1‑217 (MQ=255)
                                                                                                                 tGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAAcc                                                                                                                                                                          <  1:820013/217‑1 (MQ=255)
                                                                                                                   tAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                                                                     >  1:493891/1‑188 (MQ=255)
                                                                                                                    aGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCa                                                                                                                                                                                                    <  2:493890/188‑1 (MQ=255)
                                                                                                                     ggCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                                                                                                                                             >  1:811753/1‑242 (MQ=255)
                                                                                                                        cGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGc                                                                                                                                                                                                                          >  2:91260/1‑162 (MQ=255)
                                                                                                                        cGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGc                                                                                                                                           >  1:165213/1‑241 (MQ=255)
                                                                                                                        cGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAt                                                                                                                                    >  1:357726/1‑248 (MQ=255)
                                                                                                                         ggATAAACGTAGTGCATCAGGAAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCt                                                                                                                                                                                                                         <  1:91260/162‑1 (MQ=255)
                                                                                                                          gATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                                                                     >  1:1004497/1‑181 (MQ=255)
                                                                                                                           aTAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCa                                                                                                                                                                                                    <  2:1004496/181‑1 (MQ=255)
                                                                                                                             aaaCGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAACCCGCCAATACCTCGACGCGAAAAGCGAGGGTCCCGTGACACAGGGCATTACCGcccatcgccttacca                                                                                                        >  2:768097/1‑259 (MQ=255)
                                                                                                                                        aTCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                                                                          >  1:586966/1‑162 (MQ=255)
                                                                                                                                        aTCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAg                                                                                                                                                                                                   >  2:362771/1‑169 (MQ=255)
                                                                                                                                         tCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAg                                                                                                                                                                                                         <  2:586965/162‑1 (MQ=255)
                                                                                                                                         tCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                                                                  <  1:362772/169‑1 (MQ=255)
                                                                                                                                                    gTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                                                                     >  2:294841/1‑155 (MQ=255)
                                                                                                                                                     tAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCa                                                                                                                                                                                                    <  1:294842/155‑1 (MQ=255)
                                                                                                                                                                 aCATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTTGCAAACCGCCAATATCTCGACGAGCAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcg                                                                                                                   >  1:682745/1‑224 (MQ=255)
                                                                                                                                                                           gCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGc                                                                                                                                                   >  1:964588/1‑182 (MQ=255)
                                                                                                                                                                            cGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCg                                                                                                                                                  <  2:964587/182‑1 (MQ=255)
                                                                                                                                                                               aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGc                                                                                                                                                                                                                >  1:213484/1‑117 (MQ=255)
                                                                                                                                                                                aTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGcc                                                                                                                                                                                                               <  2:213484/117‑1 (MQ=255)
                                                                                                                                                                                         tattagagaccggTGCATGCCCTGCTTGAGCGCGCACCTGAACGTTTTCTGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCTTTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTcc                                            <  1:1058387/259‑1 (MQ=255)
                                                                                                                                                                                                  gCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGGGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGGGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCt                                  <  1:764775/272‑1 (MQ=255)
                                                                                                                                                                                                              ctgctgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCAGTGAAGATAAACGTCTGTTACCGCTGATTCACGc                                                                                                                                                                                                                >  1:293297/1‑86 (MQ=255)
                                                                                                                                                                                                               tgctgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCAGTGAAGATAAACGTCTGTTACCGCTGATTCACGcc                                                                                                                                                                                                               <  2:293296/86‑1 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTg                                                                                                                                                                                                                        >  1:642046/1‑73 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGAtt                                                                                                                                                                                                                     >  2:958983/1‑76 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGAtt                                                                                                                                                                                                                     >  2:54560/1‑76 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCCCGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGACAAACGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTAAGTATCAGGAAAACGa                                                                                  >  1:1043369/1‑207 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAg                                                                                                                                                                                                   >  1:1006910/1‑94 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAg                                                                                                                                                                                                   >  2:1016407/1‑94 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgacga                                                                                                                                                          >  1:739331/1‑135 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGAc                                                                                                                                                >  2:147032/1‑145 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                                                                        >  1:936335/1‑153 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                                                                   >  1:372944/1‑158 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                                                                   <  2:372943/158‑1 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTAt                                                                                                                         >  1:268521/1‑168 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGt                                                                                                   >  1:837675/1‑190 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGGTCTGATCGCTTCGCCCGATCCACCGTTCCTGCTGATCCTCGACGGCgacaaccatc                >  1:689100/1‑263 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGACCTGACGGATCAGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCTACGGCGTAACCGAt                 >  1:979999/1‑271 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCCGTTACCGCTGATTCACGCCCTTGAGtctc                                                                                                                                                                                                     >  2:1058502/1‑92 (MQ=255)
                                                                                                                                                                                                                   ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCAGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTTAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGCCGAGAAAAGCGACg                                                                                                                                               >  1:302782/1‑146 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGa                                                                                                                                                                                                                       <  2:642045/73‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTc                                                                                                                                                                                                                    <  1:54560/76‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTc                                                                                                                                                                                                                    <  1:958984/76‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                                                                  <  1:1016408/94‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAgg                                                                                                                                                                                                  <  2:1006909/94‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACgag                                                                                                                                                         <  2:739330/135‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACg                                                                                                                                               <  1:147032/145‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTg                                                                                                                                       <  2:936334/153‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATc                                                                                                                        <  2:268520/168‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                                                                  <  2:837674/190‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCAGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACgg                                                                                                                                              <  2:302781/146‑1 (MQ=255)
                                                                                                                                                                                                                    gAGCGCGCCCCCGAACGTTTTCATGAAGTCTTTGTTTTAAAATGCCGTTAAGATAAACGTCTGTGACCGCTGATTCACGCCCTTGAGTCTCa                                                                                                                                                                                                    <  1:1058503/92‑1 (MQ=255)
                                                                                                                                                                                                                            cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGa                                                                                                                                                                                                                       >  2:344668/1‑65 (MQ=255)
                                                                                                                                                                                                                            cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCCATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGGCAGTATCAGGAAAACGATCTGCCGGATCTGATTGCTTCGCTCGATCAACCGTTCCTGCTGAACCTGGACGGCGTAACCggtggcccttta        >  1:534081/1‑260 (MQ=255)
                                                                                                                                                                                                                             ccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGAt                                                                                                                                                                                                                      <  1:344669/65‑1 (MQ=255)
                                                                                                                                                                                                                              cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                                                                        >  2:521115/1‑142 (MQ=255)
                                                                                                                                                                                                                              cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTg                                                                              >  1:949924/1‑200 (MQ=255)
                                                                                                                                                                                                                              cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTg                                                                              <  2:949923/200‑1 (MQ=255)
                                                                                                                                                                                                                              cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGCATCAGGAAAACGATTGGCCGGATCTGAACGGTTCGCCCGAACAACCGTTCCTGCTGAACCCCCACCGCCGTAACCAACCCGcaaacccc      >  2:771556/1‑266 (MQ=255)
                                                                                                                                                                                                                               tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTg                                                                                                                                       <  1:521116/142‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGgc  <  1:893770/272‑1 (MQ=255)
                                                                                                                                                                                                                                         cAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAACCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTa                       >  2:294769/1‑243 (MQ=255)
                                                                                                                                                                                                                                          aGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAACCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTaa                      <  1:294770/243‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGAt                                     >  2:344671/1‑210 (MQ=255)
                                                                                                                                                                                                                                                              ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATc                                    <  1:344672/210‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
TTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGC  >  NC_000913/4409108‑4409604

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