Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,457,648 A→G 5.4% I73T (ATT→ACT)  yjgA ← ribosome‑associated UPF0307 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,6480AG 5.4% 229.5 / 2.9 74I73T (ATT→ACT) yjgAribosome‑associated UPF0307 family protein
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base G (2/2);  total (41/33)
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.03e-01

GCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGACCCAGATAA  >  NC_000913/4457500‑4457795
                                                                                                                                                    |                                                                                                                                                   
gcgcAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGAAGCATATTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTaa                                                                                                                                                     >  6:145436/1‑149 (MQ=255)
  gcAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATa                                                                                                                                                   >  6:187917/1‑149 (MQ=255)
       ttttCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTg                                                                                                                                              <  2:187772/149‑1 (MQ=255)
           ccAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGc                                                                                                                                           >  2:339950/1‑148 (MQ=255)
                           caccacCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCcgcgc                                                                                                                          >  1:207680/1‑149 (MQ=255)
                             ccaccTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGTCGTCTGCGACCTTCCATCTTAATACGCTTGTCGAGCTCAATAGCcgcgcg                                                                                                                         >  8:272132/1‑148 (MQ=255)
                              caccTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGCCCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCcgcgcg                                                                                                                         >  5:19487/1‑147 (MQ=255)
                                 ctgctgGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCcgcgcg                                                                                                                         >  5:339961/1‑144 (MQ=255)
                                  tgctgGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATc                                                                                                                   >  2:153582/1‑149 (MQ=255)
                                    ctgGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCg                                                                                                                 <  6:313713/149‑1 (MQ=255)
                                       gTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGc                                                                                                                <  3:118658/147‑1 (MQ=255)
                                       gTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCAt                                                                                                              >  3:106268/1‑149 (MQ=255)
                                       gTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCAt                                                                                                              <  6:339961/149‑1 (MQ=255)
                                          gtgACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCt                                                                                                            >  6:99092/1‑148 (MQ=255)
                                           tgACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTaa                                                                                                          <  2:253253/149‑1 (MQ=255)
                                              cGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGAGTAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACgg                                                                                                       <  7:272132/149‑1 (MQ=255)
                                               ggTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACgg                                                                                                       >  6:214082/1‑148 (MQ=255)
                                                gTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCGCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGa                                                                                                     <  7:52378/149‑1 (MQ=255)
                                                 tttttCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGAt                                                                                                    <  1:203660/149‑1 (MQ=255)
                                                  ttttCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATc                                                                                                   >  4:1953/1‑149 (MQ=255)
                                                   tttCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATc                                                                                                   <  2:262674/148‑1 (MQ=255)
                                                   tttCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATc                                                                                                   <  5:187917/148‑1 (MQ=255)
                                                       aGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATc                                                                                              <  3:1953/149‑1 (MQ=255)
                                                       aGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCGGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGACCTTATc                                                                                              <  1:40969/149‑1 (MQ=255)
                                                       aGCTTATCCAGCGCCTGACGAATAGGCGCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATc                                                                                              <  7:323853/149‑1 (MQ=255)
                                                           tATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAGTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAgc                                                                                          >  7:231652/1‑149 (MQ=255)
                                                           tATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAGTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAgc                                                                                          >  7:237492/1‑149 (MQ=255)
                                                             tCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAgcgc                                                                                        >  6:303944/1‑149 (MQ=255)
                                                               cAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTGATCCAGCGCGt                                                                                      >  1:292813/1‑149 (MQ=255)
                                                                   gcCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGttttt                                                                                  <  1:115318/149‑1 (MQ=255)
                                                                    ccTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTc                                                                                 >  3:296153/1‑149 (MQ=255)
                                                                     cTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGCTCCGCATCTTACGGGATCTTATCCAGCGCGTTTTTcc                                                                                >  3:99030/1‑149 (MQ=255)
                                                                     cTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTcc                                                                                >  6:347993/1‑149 (MQ=255)
                                                                      tGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTccc                                                                               >  3:348549/1‑149 (MQ=255)
                                                                      tGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTccc                                                                               >  8:323155/1‑149 (MQ=255)
                                                                      tGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTccc                                                                               >  6:270590/1‑149 (MQ=255)
                                                                        aCGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTccc                                                                               <  4:296153/147‑1 (MQ=255)
                                                                            aTAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAg                                                                                                                                <  1:325926/94‑1 (MQ=255)
                                                                            aTAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAg                                                                                                                                >  2:325926/1‑94 (MQ=255)
                                                                                  tctACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAAt                                                                   >  2:163472/1‑149 (MQ=255)
                                                                                   ctACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAAtt                                                                  >  7:32339/1‑149 (MQ=255)
                                                                                           cgcTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCgcgc                                                           >  8:343641/1‑148 (MQ=255)
                                                                                           cgcTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGcc                                                          <  7:343641/149‑1 (MQ=255)
                                                                                            gcTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGcc                                                          >  3:215120/1‑148 (MQ=255)
                                                                                            gcTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGcc                                                          >  3:210214/1‑148 (MQ=255)
                                                                                                  gcagcaTCTTACCGATGAGCTGCAACTGGCGGCGGCGCCCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGt                                                   >  6:381332/1‑149 (MQ=255)
                                                                                                    agcaTCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTGAATAGCCGCGCGCAGATCCGCATCTAACTGTATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGttt                                                 <  6:19487/149‑1 (MQ=255)
                                                                                                     gcaTCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGtttt                                                <  4:348549/149‑1 (MQ=255)
                                                                                                              cGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCg                                       >  1:202269/1‑149 (MQ=255)
                                                                                                              cGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCg                                       >  2:16022/1‑149 (MQ=255)
                                                                                                                 tGAGCTGCAACTGGCGGCGGCGACCGTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCa                                    >  1:52406/1‑149 (MQ=255)
                                                                                                                   aGCTGCAACTGGCGGCGGCGACCTTCCATCTTACTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATc                                  <  4:348368/149‑1 (MQ=255)
                                                                                                                    gCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCa                                 >  5:83330/1‑149 (MQ=255)
                                                                                                                    gCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCa                                 >  7:125222/1‑149 (MQ=255)
                                                                                                                    gCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGAGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCa                                 >  7:181011/1‑149 (MQ=255)
                                                                                                                        cAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGtt                             <  1:163472/149‑1 (MQ=255)
                                                                                                                          aCTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTa                           <  1:339950/149‑1 (MQ=255)
                                                                                                                              gcggcggcgACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAAtt                        >  1:226196/1‑148 (MQ=255)
                                                                                                                               cggcggcgACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAAttt                       >  6:364931/1‑148 (MQ=255)
                                                                                                                               cggcggcgACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTc                      <  2:207680/149‑1 (MQ=255)
                                                                                                                               cggcggcgACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTc                      <  5:381332/149‑1 (MQ=255)
                                                                                                                                    gcgACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACttt                 <  5:214082/149‑1 (MQ=255)
                                                                                                                                    gcgACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACttt                 <  8:125222/149‑1 (MQ=255)
                                                                                                                                      gACCTTCCATCTTAATACGCTGGGCGAGCTCAATTGCCGCGCTCAGATCCGGATCTAAGGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCTCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTa               <  1:111954/149‑1 (MQ=255)
                                                                                                                                      gACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTa               <  5:364931/149‑1 (MQ=255)
                                                                                                                                      gACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTa               <  1:200269/149‑1 (MQ=255)
                                                                                                                                      gACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTa               <  6:83330/149‑1 (MQ=255)
                                                                                                                                      gACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTa               >  3:159396/1‑149 (MQ=255)
                                                                                                                                       aCCTTCCATCTTAGTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTAc              <  8:231652/149‑1 (MQ=255)
                                                                                                                                       aCCTTCCATCTTAGTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTAc              <  8:237492/149‑1 (MQ=255)
                                                                                                                                          ttCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGa           >  5:305587/1‑149 (MQ=255)
                                                                                                                                              aTCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGa           >  6:44658/1‑145 (MQ=255)
                                                                                                                                              aTCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGa           <  5:44658/145‑1 (MQ=255)
                                                                                                                                               tCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGAcc         >  1:168826/1‑146 (MQ=255)
                                                                                                                                               tCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGAcc         >  3:209733/1‑146 (MQ=255)
                                                                                                                                                   aaTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGACCCAGAt    <  4:159396/147‑1 (MQ=255)
                                                                                                                                                    aTACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGACCCAGATaa  <  8:32339/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTTTTCCCCAGATCAACAATTTCCGCGCCAAGGCGTTTTAGCTCCTCGGCATCACGTTTAATTTCACTTTTACTGACCCAGATAA  >  NC_000913/4457500‑4457795

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