Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 59,968 C→T V127I (GTC→ATC)  rluA ← 23S rRNA pseudouridine(746) and tRNA pseudouridine(32) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,9680CT98.2% 177.4 / NA 57V127I (GTC→ATC) rluA23S rRNA pseudouridine(746) and tRNA pseudouridine(32) synthase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (40/16);  total (41/16)
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.94e-01

CCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCAC  >  NC_000913/59724‑60211
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
ccATACGCCGGATGGGTAATCGTCACCATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCtt                                                                                                                                                                                                                                               >  2:100003/1‑250 (MQ=255)
ccATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTg                                                                                                                                                                                                                                               <  2:511462/251‑1 (MQ=255)
  aTACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGtt                                                                                                                                                                                                                                             <  2:77879/251‑1 (MQ=255)
       ccGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGCGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGCACGCCCGGGAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCCCTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTctgg                                                                                                                                                                                                                                        >  2:406560/1‑248 (MQ=255)
         ggATGGGTAATCGTCAAAATCTCTGCACGTACCAACAAACGTGGTGCCATAGCTATCGATTCTGGTGATGAATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCAAAAGAAATTGATGCGAACGCACGGTAATCGGTTTTAACACCACTCTTGCCGAGTTATACGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCg                                                                                                                                                                                                                                      >  2:376688/1‑251 (MQ=255)
                tAATCGTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACCGATTTTCTGTTTCTGGCGGTTTGGc                                                                                                                                                                                                                               >  2:562131/1‑251 (MQ=255)
                    gtccAACATCTCTTGAGGTAGCAACAAACGTGGTGCCATCGCTCGCGCTTCGGGTGATGCATAAAAACGAGCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACGCTTGCCGTGTTATACGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGGTTACCCGGTTCGTAACAGATGTTCTGTTTCGGGCGGTTTGGCCAGt                                                                                                                                                                                                                           <  2:541989/248‑1 (MQ=255)
                     gTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGTGCGGTTTGGCCAGTc                                                                                                                                                                                                                          >  1:552311/1‑251 (MQ=255)
                      tCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCg                                                                                                                                                                                                                         >  1:426752/1‑251 (MQ=255)
                                     gTAGCAACAAACGTGGTGCCATCGATCTCGCTTCTGGGGATGAATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACCCGCAATTGATGCGAACGCCCGGTAATCGGTTTTACCACCACTCTTGCCGTGTTATACACCGCATACTCCACCCCTTCATATTCCGTATGCGCAGTTTTACCCATTTCGTACCAGATTTTCTGTTTCAGGCGGTTTGCCCAGTCGCCAAtaggcggcggc                                                                                                                                                                                                          >  2:18423/1‑241 (MQ=255)
                                      tAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGc                                                                                                                                                                                                             >  1:393192/1‑247 (MQ=255)
                                         caacaaACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTGATCCGCCGCATACTCAACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGTGCGTGTTGTCCAGTCGCAAATCAGCGGCAGATccc                                                                                                                                                                                                      >  2:133960/1‑250 (MQ=255)
                                              aaCGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAACACGATCGCCGAGAATCGGATGCCCCAGCGCCAGCATATGCACACGCAATTGATGCGAAAGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCCGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATccacc                                                                                                                                                                                                    >  2:92931/1‑248 (MQ=255)
                                                       ccATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCg                                                                                                                                                                                        >  2:131385/1‑251 (MQ=255)
                                                            gCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCggtgg                                                                                                                                                                                    >  1:8027/1‑247 (MQ=255)
                                                                                             ccGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCtt                                                                                                                                                  >  1:102783/1‑251 (MQ=255)
                                                                                               gagaATCGGATGACACAGCGCCAGCATATGCACACGCAATTGATGCGCACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACCTACTGCtttt                                                                                                                                                >  2:240352/1‑251 (MQ=255)
                                                                                                         tGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTCTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCAATTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACCGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGGTGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcg                                                                                                                                      >  2:189259/1‑251 (MQ=255)
                                                                                                            cccAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcgctc                                                                                                                                   >  1:262515/1‑251 (MQ=255)
                                                                                                             ccacggcCAGCCTGTGCCCACGCAATCGATGCGAACGGCCGTTAATCTGTTTTAACACCACGCTTGCCGTGTCATCCTCCGCAGACTCCACCACTTCATACTCCTTCAGCGCACGTTTTCCCGTTTCGTAACAGATTTTCTGTGTCGGGCGTTTTGGCCAGTCGCAAATCAGCTGCAGAGCCACCAGATCTTCTGCGTGGTATGTATGCCCCCAGACGCGGGCCACATACTGCTTGGTCGgctgacgctcg                                                                                                                                  <  2:248756/246‑1 (MQ=255)
                                                                                                                ccgcCAGCATATGCACACGCAATTGATGCGAACTCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCTCCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCgg                                                                                                                               <  1:146406/250‑1 (MQ=255)
                                                                                                                  gcCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGa                                                                                                                              >  1:255859/1‑250 (MQ=255)
                                                                                                                    cAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACt                                                                                                                           >  1:205008/1‑251 (MQ=255)
                                                                                                                            gCACACGCAATTGATGCGAACGCCCGGTACTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAAATGGCGtttt                                                                                                                   >  2:47670/1‑251 (MQ=255)
                                                                                                                                acGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTGTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACt                                                                                                               <  1:92931/251‑1 (MQ=255)
                                                                                                                                              aaCGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCt                                                                                                 >  2:33107/1‑251 (MQ=255)
                                                                                                                                                     ggtaagtggTTTTAACACCACTCTTGCCGTGTTATCCTCCGCATACTCCTCCACTGCATATTCCGTCGGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGTCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGATTTGGTCa                                                                                          <  2:8027/244‑1 (MQ=255)
                                                                                                                                                      gTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTc                                                                                           >  2:240359/1‑249 (MQ=255)
                                                                                                                                                         aTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAgcgc                                                                                      >  1:12947/1‑251 (MQ=255)
                                                                                                                                                            ggTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAACTCAGTGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCAc                                                                                   >  2:12193/1‑251 (MQ=255)
                                                                                                                                                              ttttAACACCACTCTTGCCGTGTTATCCGCCGCCTACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACt                                                                                  >  2:300337/1‑250 (MQ=255)
                                                                                                                                                                ttAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCCCTAc                                                                                >  2:254247/1‑250 (MQ=255)
                                                                                                                                                                   acacCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATc                                                                            >  1:380433/1‑251 (MQ=255)
                                                                                                                                                                    caccacTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCa                                                                           <  1:356188/251‑1 (MQ=255)
                                                                                                                                                                      ccacTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGTCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCAcg                                                                         >  1:440516/1‑251 (MQ=255)
                                                                                                                                                                      ccacTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCAcg                                                                         >  1:468333/1‑251 (MQ=255)
                                                                                                                                                                       cacTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCActc                                                                        >  1:340069/1‑249 (MQ=255)
                                                                                                                                                                                cGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCACGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAAATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTAAGCGCCACTACAATCACGCCGCTGGt                                                                 >  1:357180/1‑249 (MQ=255)
                                                                                                                                                                                cGTGTTATCCGCCGAATCATCAACAACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGAGCTCGCGGAAATGGCGTTTTAACTCCCGCTCCTCGGCTTTGGTCAGCGCCACTACAATCCCGCCGCTGGTac                                                               >  2:259751/1‑250 (MQ=255)
                                                                                                                                                                                 gtgttatccgaagcctacgacaaaacTTCATCTTCAGTCTGCGCAGGTTTAACCGCTTCGTAACAAACTTTCTGTTTCGCGCGGTTTGGCCAGCCGCAAATAAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGCTGCACCCAGCCGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACCGGCGTTTTCACTCCCGCTCCGCGGCCTTGGTCCGCGCCACTACTATCACGCCGCTGGTAGc                                                              >  2:240092/25‑251 (MQ=255)
                                                                                                                                                                                        ccgccgCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATcc                                                       >  1:508296/1‑251 (MQ=255)
                                                                                                                                                                                        ccgccgCATACTCCACCACTGCAAATTCCGGCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTTTTTCGGGCGGTGTGGCCAGTCGCAAATCAGCGGCAGAGCCCCCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGTGCCACATACTGCTTTTTCGGCTCGCGGTCGCGGAACGGGCGTTTTAACTCCCGCTCCGCGGCTTTGGGCAGCGCCACTGCAATCACGCCGCTGGTAGCCATATcc                                                       <  2:270849/251‑1 (MQ=255)
                                                                                                                                                                                           ccgcATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGa                                                    <  1:42047/251‑1 (MQ=255)
                                                                                                                                                                                            cgcATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGAc                                                   >  1:262176/1‑251 (MQ=255)
                                                                                                                                                                                                    ccaccaCTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTCCAATCACGCCGCTGGTAGCCATATCCAGACGATGCAcc                                           >  1:506879/1‑251 (MQ=255)
                                                                                                                                                                                                     caccacTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCCTATCCAGACGATGCACCg                                          >  1:286228/1‑251 (MQ=255)
                                                                                                                                                                                                          cTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGGAGCCATATCCAGACGATGCACCGATTct                                     >  2:191092/1‑251 (MQ=255)
                                                                                                                                                                                                                atTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTctgcctgc                               <  1:345459/251‑1 (MQ=255)
                                                                                                                                                                                                                        tGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATc                       >  1:566932/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcgcAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATc                       <  2:508296/250‑1 (MQ=255)
                                                                                                                                                                                                                          cgcAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCAc                     <  1:131385/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggTTTACCAGTTTCGTAACAGATTTTCTGTTTAGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGGGCTCGCGGAACGGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGGTTCTGACTGCGGATAATCCCGCTg                 >  2:497900/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tACCCGTTGCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACCTACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCCCCGATTCTGCCTGCGGATAATCTCGCTGAAt              <  1:133960/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aCCCGTTTCGTAACAGATTTGCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACGGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATgc            <  1:334299/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tttCGCAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCCGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGtc       >  1:121999/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGAACTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGATTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCATGCCGCTGGTAGCCATATCCAGCCGATGCCCCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGtcatc    >  2:110998/1‑251 (MQ=255)
                                                                                                                                                                                                                                              atcAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTACGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACTCTGAATGCGCGTCATCAc  <  1:157949/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGGTGCCATCGCTCTCGCTTCTGGTGATGCATAAAAACGATCGCCGAGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCAC  >  NC_000913/59724‑60211

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