Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 60,254 A→G 6.3% G31G (GGT→GGC rluA ← dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091360,2540AG 6.3% 422.8 / 9.2 128G31G (GGT→GGCrluAdual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (59/60);  new base G (5/3);  total (65/63)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 135 total aligned reads displayed.

CTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGAC  >  NC_000913/60107‑60398
                                                                                                                                                   |                                                                                                                                                
cTCCGCGGCTTTGGTCAGCGCCACTACCATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGGCGACCCGGCACTGACAACaaa                                                                                                                                                  <  2:551743/148‑1 (MQ=255)
          ttGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAGCCgagcggtt                                                                                                                                        >  3:291226/1‑141 (MQ=255)
             gTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgt                                                                                                                                     <  2:1120366/148‑1 (MQ=255)
                aGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACaaa                                                                                                                                                  <  1:1104043/132‑1 (MQ=255)
                 gcgcCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAGCCGAGCGGTttgttg                                                                                                                                   >  4:95064/1‑146 (MQ=255)
                 gcgcCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAGCCGAGCGGTttgttg                                                                                                                                   <  3:95064/146‑1 (MQ=255)
                             tCACGCCGCTCGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAatat                                                                                                                      <  1:1082515/147‑1 (MQ=255)
                             tCACGCCGCTCGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAatat                                                                                                                      >  2:1082515/1‑147 (MQ=255)
                                 gccgcTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCa                                                                                                                 <  4:21478/148‑1 (MQ=255)
                                  ccgcTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                <  3:730501/148‑1 (MQ=255)
                                     cTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCt                                                                                                             >  1:91866/1‑148 (MQ=255)
                                             cATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATcc                                                                                                              <  3:839329/139‑1 (MQ=255)
                                             cATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATcc                                                                                                              >  4:839329/1‑139 (MQ=255)
                                             cATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAgg                                                                                                     <  2:91866/148‑1 (MQ=255)
                                                 tCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATaa                                                                                                 <  2:435462/148‑1 (MQ=255)
                                                  ccAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgt                                                                                                                                     >  4:1108080/1‑111 (MQ=255)
                                                  ccAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgt                                                                                                                                     <  3:1108080/111‑1 (MQ=255)
                                                    aGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATcc                                                                                                              >  4:868039/1‑132 (MQ=255)
                                                    aGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTg                                                                                                            >  3:172864/1‑134 (MQ=255)
                                                    aGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTg                                                                                                            <  4:172864/134‑1 (MQ=255)
                                                    aGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATaacca                                                                                              >  2:169585/1‑148 (MQ=255)
                                                     gACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATaaccaa                                                                                             >  4:250064/1‑148 (MQ=255)
                                                      aCGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttg                                                                                                                                      >  4:791074/1‑106 (MQ=255)
                                                      aCGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttg                                                                                                                                      <  3:791074/106‑1 (MQ=255)
                                                      aCGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACa                                                                                                       <  4:486618/137‑1 (MQ=255)
                                                         aTGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCga                                                                                                                                              <  2:848479/95‑2 (MQ=255)
                                                         aTGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCga                                                                                                                                              >  1:848479/1‑94 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGTCAACAAGCCGAGCGGTttgttgtc                                                                                                                                 >  1:288086/1‑101 (MQ=38)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGTCAACAAGCCGAGCGGTttgttgtc                                                                                                                                 <  2:288086/103‑3 (MQ=38)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCata                                                                                                                          <  4:799717/110‑1 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                >  4:709342/1‑120 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                >  1:576185/1‑120 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                <  3:709342/120‑1 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                <  2:576185/120‑1 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATa                                                                                                         <  4:306640/127‑1 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATa                                                                                                         >  3:306640/1‑127 (MQ=255)
                                                              ccGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTc                                                                                    >  2:169285/1‑148 (MQ=255)
                                                                         tgcGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGAGCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATc                                                                                                               >  1:965826/1‑110 (MQ=255)
                                                                         tgcGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATc                                                                                                               <  2:965826/110‑1 (MQ=255)
                                                                              aTAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAg                                                                    >  3:458679/1‑148 (MQ=255)
                                                                               tAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGaccacc                                                                                                                             <  2:576247/90‑1 (MQ=255)
                                                                               tAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGaccacc                                                                                                                             >  1:576247/1‑90 (MQ=255)
                                                                               tAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGt                                                                   <  3:250064/148‑1 (MQ=255)
                                                                                 aTCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATg                                                                                                                     >  3:635326/1‑96 (MQ=255)
                                                                                   cACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                   >  4:800192/1‑80 (MQ=255)
                                                                                   cACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                   <  3:800192/80‑1 (MQ=255)
                                                                                   cACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGt                                                                      >  1:540543/1‑141 (MQ=255)
                                                                                   cACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTc                                                               <  1:169285/148‑1 (MQ=255)
                                                                                         gAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGa                                                                                                           <  2:751828/98‑1 (MQ=255)
                                                                                         gAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGa                                                                                                           >  1:751828/1‑98 (MQ=255)
                                                                                         gAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGa                                                                          <  3:486125/131‑1 (MQ=255)
                                                                                         gAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGa                                                                          >  4:486125/1‑131 (MQ=255)
                                                                                          aaTGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                >  1:890753/1‑92 (MQ=255)
                                                                                          aaTGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCAt                                                                                                                <  2:890753/92‑1 (MQ=255)
                                                                                              cgcgTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGttt                                                                 >  3:862449/1‑135 (MQ=255)
                                                                                              cgcgTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGttt                                                                 <  4:862449/135‑1 (MQ=255)
                                                                                                     tcaCGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTc                                             >  3:583823/1‑148 (MQ=255)
                                                                                                      caCGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCAtttt                                                   >  1:642750/1‑141 (MQ=255)
                                                                                                      caCGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTcc                                            >  2:155676/1‑148 (MQ=255)
                                                                                                          cTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttg                                                                                                                                      >  4:694062/1‑54 (MQ=255)
                                                                                                          cTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCt                                                                                                                                        <  4:652079/52‑1 (MQ=255)
                                                                                                          cTGTCTTTGTACTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCt                                                                                                                                        >  3:652079/1‑52 (MQ=255)
                                                                                                            gTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTc                                             <  3:960806/141‑1 (MQ=255)
                                                                                                            gTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTc                                             >  4:960806/1‑141 (MQ=255)
                                                                                                             tCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTgg                                                                           <  4:391487/110‑1 (MQ=255)
                                                                                                              cTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                   <  1:613422/53‑1 (MQ=255)
                                                                                                              cTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCttgttg                                                                                                                                   >  2:613422/1‑53 (MQ=255)
                                                                                                                ttGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTg                                                                                                            >  3:476210/1‑74 (MQ=255)
                                                                                                                ttGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTg                                                                                                            <  4:476210/74‑1 (MQ=255)
                                                                                                                ttGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATccc                                                         >  2:644793/1‑125 (MQ=255)
                                                                                                                 tgtgCTCTTCCAGACGACCCGGCACTGACAACAAACCGCACGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGc                                 <  1:155676/148‑1 (MQ=255)
                                                                                                                    gCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGtt                              <  1:169585/148‑1 (MQ=255)
                                                                                                                      tctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTAc                            >  1:1010072/1‑148 (MQ=255)
                                                                                                                       ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGTATCAGGAtg                                                                                                  <  1:85395/77‑2 (MQ=255)
                                                                                                                       ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGTATCAGGAtg                                                                                                  >  2:85395/1‑76 (MQ=255)
                                                                                                                         tCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAc                         >  2:594469/1‑148 (MQ=255)
                                                                                                                           cAGACGACCCGGCACTGACAACAAACCGCTCGGCATGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGCTCCTGCGGTGGCTTGCAGTTTGCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAccg                       >  2:34770/1‑146 (MQ=255)
                                                                                                                            aGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGAt                      <  2:1010072/148‑1 (MQ=255)
                                                                                                                            aGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGAt                      <  1:34770/148‑1 (MQ=255)
                                                                                                                              acgacCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATg                                  >  3:881159/1‑134 (MQ=255)
                                                                                                                               cgacCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCaa                   >  2:796170/1‑148 (MQ=255)
                                                                                                                                 acCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAAc                 >  4:208988/1‑148 (MQ=255)
                                                                                                                                    cGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCa              <  1:796170/148‑1 (MQ=255)
                                                                                                                                          tGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCAt        <  3:208988/148‑1 (MQ=255)
                                                                                                                                          tGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCAt        <  2:1101179/148‑1 (MQ=255)
                                                                                                                                           gACAACAAACCGCTCGGCTTGTTGACCACCATACTATGGTCATCCTGATAc                                                                                                        >  3:1114338/1‑51 (MQ=255)
                                                                                                                                           gACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCg                                <  1:630769/123‑1 (MQ=255)
                                                                                                                                           gACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCg                                >  2:630769/1‑123 (MQ=255)
                                                                                                                                           gACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTAATAc                                                                                                        <  4:1114338/51‑1 (MQ=255)
                                                                                                                                             caacaaGCCGAGCGGTTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATc                                                           >  1:25245/1‑94 (MQ=255)
                                                                                                                                             caacaaGCCGAGCGGTTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATc                                                           <  2:25245/94‑1 (MQ=255)
                                                                                                                                             caacaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCAtttt                                                   >  1:1109061/1‑102 (MQ=255)
                                                                                                                                             caacaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCAtttt                                                   <  2:1109061/102‑1 (MQ=255)
                                                                                                                                               acaaGCCGAGCGGTTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGa   >  2:540501/1‑148 (MQ=255)
                                                                                                                                               acaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCAtttt                                                   >  2:999107/1‑100 (MQ=255)
                                                                                                                                                caaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGAc  >  1:1094024/1‑148 (MQ=255)
                                                                                                                                                 aaaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGa                       <  1:201691/126‑1 (MQ=255)
                                                                                                                                                 aaaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGa                       >  2:201691/1‑126 (MQ=255)
                                                                                                                                                  aaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGtt                                                                  <  2:649034/82‑1 (MQ=255)
                                                                                                                                                  aaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGtt                                                                  >  1:649034/1‑82 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTCCGCGGCTTTGGTCAGCGCCACTACAATCACGCCGCTGGTAGCCATATCCAGACGATGCACCGATTCTGCCTGCGGATAATCACGCTGAATGCGCGTCATCACGCTGTCTTTGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGAC  >  NC_000913/60107‑60398

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