Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,257,513 G→A 5.3% S285S (AGC→AGT psuG ← pseudouridine 5'‑phosphate glycosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,257,5130GA 5.3% 238.9 / 2.2 75S285S (AGC→AGTpsuGpseudouridine 5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base A (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

GGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATC  >  NC_000913/2257365‑2257660
                                                                                                                                                    |                                                                                                                                                   
ggCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAgg                                                                                                                                                     <  2:108834/149‑1 (MQ=255)
ggCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAgg                                                                                                                                                     <  8:46158/149‑1 (MQ=255)
 gCGGCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTGCTGCCCGTCAGAAATTAACCCGCGAGACGATGATATTCTTTGGCCAGTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGc                                                                                                                                                    <  2:41014/149‑1 (MQ=255)
          ccAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCAccgc                                                                                                                                           >  3:162360/1‑149 (MQ=255)
              gTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGt                                                                                                                                       >  8:42201/1‑149 (MQ=255)
              gTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGt                                                                                                                                       <  7:194752/149‑1 (MQ=255)
                     tATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGtt                                                                                                                                  >  5:287012/1‑147 (MQ=255)
                        ccTTGATGAGTACGGGCGCCGGCTTTCTGCCCGTCAGATATTACCCCGCGAGACGCGCATAATCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGATGGATGTGGGAGTTCAGGATGTCGCCGCCGGTCAGTTAAGc                                                                                                                              >  6:315639/2‑148 (MQ=255)
                            gatgaGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAgg                                                                                                                                                     <  1:409285/121‑1 (MQ=255)
                            gatgaGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAgg                                                                                                                                                     >  2:409285/1‑121 (MQ=255)
                             atgaGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAg                                                                                                                                    >  7:175336/1‑137 (MQ=255)
                             atgaGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAg                                                                                                                                    <  8:175336/137‑1 (MQ=255)
                                  tACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACgcgcgc                                                                                                                     >  3:283344/1‑147 (MQ=255)
                                      ggCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCcagcag                                                                                                               <  6:341005/149‑1 (MQ=255)
                                          aCGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGACTGTc                                                                                                                                                <  7:1700/112‑1 (MQ=255)
                                          aCGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGACTGTc                                                                                                                                                >  8:1700/1‑112 (MQ=255)
                                          aCGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGACTGTc                                                                                                                                                <  7:8084/112‑1 (MQ=255)
                                          aCGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGACTGTc                                                                                                                                                >  8:8084/1‑112 (MQ=255)
                                           cGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATg                                                                                                          <  4:162360/149‑1 (MQ=255)
                                           cGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATg                                                                                                          <  1:241371/149‑1 (MQ=255)
                                            ggCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAgg                                                                                                                                                     >  2:21380/1‑105 (MQ=255)
                                             gCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGc                                                                                                        <  7:70048/149‑1 (MQ=255)
                                             gCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGc                                                                                                        <  7:63058/149‑1 (MQ=255)
                                               tttCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCg                                                                                                       >  4:336632/1‑148 (MQ=255)
                                                 tCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTa                                                                                                     <  3:141496/148‑1 (MQ=255)
                                                     cccGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttc                                                                                                <  1:332906/149‑1 (MQ=255)
                                                         tCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTa                                                                                            >  1:41675/1‑149 (MQ=255)
                                                         tCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTa                                                                                            >  1:263027/1‑149 (MQ=255)
                                                         tCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACAGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGAGTACTTTCTTTa                                                                                            >  1:10066/1‑149 (MQ=255)
                                                           aGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTAc                                                                                           <  6:287012/148‑1 (MQ=255)
                                                            gAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTAcc                                                                                          >  1:389790/1‑148 (MQ=255)
                                                                    cccGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTcc                                                                                 <  7:42201/149‑1 (MQ=255)
                                                                    cccGCGAGAAGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTAGAACACCAGCTGGATGTTGGATTTCAGGCTGTCCCCGCCGGTCAGTTCAGCAACTCGCGCCAGCAGGAATGGCGTCCTTTCTTTACAAATGACTca                                                                                 >  7:365966/1‑148 (MQ=255)
                                                                       gcgAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTa                                                                                            >  5:61516/1‑135 (MQ=255)
                                                                       gcgAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTa                                                                                            <  6:61516/135‑1 (MQ=255)
                                                                       gcgAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTcc                                                                                 >  7:94688/1‑146 (MQ=255)
                                                                       gcgAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTcc                                                                                 <  8:94688/146‑1 (MQ=255)
                                                                             cGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttc                                                                        <  4:283344/149‑1 (MQ=255)
                                                                             cGCTGATATTCTTTGGCAAATTCGCTCGCCAGACTGGCGTTGTTGAACACCAGCTGGATGTTGGATTGCAGGCTGTCACCGCCGGTCAGTGCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttc                                                                        <  3:293828/149‑1 (MQ=255)
                                                                              gCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttctt                                                                                              >  2:367512/1‑126 (MQ=255)
                                                                              gCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttctt                                                                                              <  1:367512/126‑1 (MQ=255)
                                                                                tGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttcagc                                                                     >  4:240720/1‑149 (MQ=255)
                                                                                tGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttcagc                                                                     >  8:38083/1‑149 (MQ=255)
                                                                                 gATATTCTTTGGCAATTTCGCTAGCCACAATGGCGTTGAACAACACAAGCTGCATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGTAATCCCGAACTTTCTTTACCAATGACTCCCTGagcttcagct                                                                    >  6:20941/1‑149 (MQ=255)
                                                                                   tatTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttcagctt                                                                   >  6:89582/1‑148 (MQ=255)
                                                                                     ttCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTAcc                                                                                          <  7:242578/123‑1 (MQ=255)
                                                                                     ttCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTAcc                                                                                          >  8:242578/1‑123 (MQ=255)
                                                                                     ttCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAgcttcagcttca                                                                 <  3:336632/148‑1 (MQ=255)
                                                                                       cTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCttcagcttcagc                                                               >  6:382156/1‑148 (MQ=255)
                                                                                        tttGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACCCCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGc                                                                                                        >  4:27202/1‑106 (MQ=255)
                                                                                        tttGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGGCAGTTCAGCAACGCGCGCCAGCAGGAATGGc                                                                                                        <  3:27202/106‑1 (MQ=255)
                                                                                          ttgCAATTTGGCTCGCCAGAATGGAGTTGATTAACACCAGGTGGTTGTTGGATTTAAGGCTGTCAGCGGCGGTCAGTTCAGGAATGCGCGCCAGCAGGACTGGCGTACTTTCTTTACCAATGACGCCCTGAGCTTAAGCTTCAGCTAc                                                            <  8:365966/146‑1 (MQ=255)
                                                                                          tGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTAc                                                            <  5:89582/148‑1 (MQ=255)
                                                                                            gCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGc                                                         <  5:20941/149‑1 (MQ=255)
                                                                                             cAATTTCGCTCGCGAGATTGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGGGCCAGCAGGTTTGGCGTACTTTCTTTACGAATGACTCCCTGAGCTTCAGCTTCAGCTACCGcc                                                        <  6:13390/149‑1 (MQ=255)
                                                                                              aaTTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCGCTGAGCTTCAGCTTCAGCTACCGcc                                                        >  4:90962/1‑148 (MQ=255)
                                                                                              aaTTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGcc                                                        <  3:90962/148‑1 (MQ=255)
                                                                                                    gctcgcCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTgat                                                    >  2:379667/1‑146 (MQ=255)
                                                                                                    gctcgcCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTgat                                                    <  1:379667/146‑1 (MQ=255)
                                                                                                     ctcgcCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTgatcgat                                                <  1:183913/149‑1 (MQ=255)
                                                                                                         ccAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATcgcc                                            <  2:263027/149‑1 (MQ=255)
                                                                                                                  cGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATa                                   <  2:389790/149‑1 (MQ=255)
                                                                                                                    ttgttgAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAgt                                 <  7:410106/149‑1 (MQ=255)
                                                                                                                    ttgttgAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAgt                                 <  5:382156/149‑1 (MQ=255)
                                                                                                                       ttgAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCTTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAgtgtg                              >  2:89800/1‑149 (MQ=255)
                                                                                                                          aaCACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTc                           >  1:266536/1‑149 (MQ=255)
                                                                                                                               cAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttcagc                                                                     >  4:223177/1‑102 (MQ=255)
                                                                                                                               cAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagcttcagc                                                                     <  3:223177/102‑1 (MQ=255)
                                                                                                                               cAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCAtt                                       >  3:37587/1‑132 (MQ=255)
                                                                                                                               cAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCAtt                                       <  4:37587/132‑1 (MQ=255)
                                                                                                                                 gCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATa                    <  3:240720/149‑1 (MQ=255)
                                                                                                                                    ggATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGc                  >  6:204718/1‑148 (MQ=255)
                                                                                                                                         ttGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttctt                                                                                              <  1:14676/67‑1 (MQ=255)
                                                                                                                                         ttGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttctt                                                                                              >  2:14676/1‑67 (MQ=255)
                                                                                                                                         ttGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTgat                                                    <  1:425818/109‑1 (MQ=255)
                                                                                                                                         ttGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTgat                                                    >  2:425818/1‑109 (MQ=255)
                                                                                                                                              tttCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGCCTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATa                                   >  2:352571/1‑121 (MQ=255)
                                                                                                                                              tttCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGCCTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATa                                   <  1:352571/121‑1 (MQ=255)
                                                                                                                                              tttCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCg       >  2:256169/1‑149 (MQ=255)
                                                                                                                                                  aGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACATTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGCTCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCggcaa   >  5:45702/1‑146 (MQ=255)
                                                                                                                                                   ggCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAgcttcagcttca                                                                 >  6:410017/1‑86 (MQ=255)
                                                                                                                                                   ggCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATc  <  2:41675/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATC  >  NC_000913/2257365‑2257660

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