Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,451,404 T→G 5.4% D94A (GAC→GCC)  rplW ← 50S ribosomal subunit protein L23

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,451,4040TG 5.4% 436.2 / 12.2 125D94A (GAC→GCC) rplW50S ribosomal subunit protein L23
Reads supporting (aligned to +/- strand):  ref base T (59/59);  new base G (4/3);  total (63/62)
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.99e-01

Only 100 of 127 total aligned reads displayed.

TTTTTTCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGC  >  NC_000913/3451256‑3451552
                                                                                                                                                    |                                                                                                                                                    
ttttttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGt                                                                                                                                                      <  2:824916/149‑1 (MQ=255)
 tttttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTc                                                                                                                                                     <  2:675501/149‑1 (MQ=255)
  ttttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTcc                                                                                                                                                    <  2:150148/149‑1 (MQ=255)
          cAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTg                                                                                                                                            <  4:790718/149‑1 (MQ=255)
                   aaaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCttt                                                                                                                                   >  2:437875/1‑149 (MQ=255)
                    aaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                  >  4:899034/1‑149 (MQ=255)
                     aaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCa                                                                                                                                 <  2:925547/149‑1 (MQ=255)
                      aGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGATATGTCGGTTTACATTTAACAACTGCTATTGTATTACTCCACCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAg                                                                                                                                <  4:273108/149‑1 (MQ=255)
                      aGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAg                                                                                                                                <  2:756586/149‑1 (MQ=255)
                      aGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAg                                                                                                                                <  4:604173/149‑1 (MQ=255)
                       ggTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAgg                                                                                                                               <  1:383295/149‑1 (MQ=255)
                              ccTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGt                                                                                                                        >  4:138359/1‑149 (MQ=255)
                               cTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTa                                                                                                                       >  4:872288/1‑149 (MQ=255)
                               cTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTa                                                                                                                       >  4:541341/1‑149 (MQ=255)
                                ttGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTaa                                                                                                                      >  3:603855/1‑149 (MQ=255)
                                 tgtgCAGCTTAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATGCTGGCCTTCTTTCAGGGTGACGTAAg                                                                                                                     <  3:638304/149‑1 (MQ=255)
                                  gtgCAGCTCAGGGTTAACCATTTTAACTACGTGGCGACGACTCGGTGATGTCGTTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACGTACTCAGCGCCGCAAACGAAGTCCAGCTTTTGTCCTTCTTTCGGGGTGACGCAAGc                                                                                                                    <  2:324473/149‑1 (MQ=255)
                                        cTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttttt                                                                                                              <  3:769572/149‑1 (MQ=255)
                                           aGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCa                                                                                                           <  3:803631/149‑1 (MQ=255)
                                           aGGGTTAAACAATTTAACTACGTTGCGACGACCCGGAGATGTCGGGTTACCTTTTACAACTTCCATTTTATTTCTTCTCCGACTTACTCAGCGCCGGCAACGGAAGCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCa                                                                                                           <  4:200770/149‑1 (MQ=255)
                                             ggTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGt                                                                                                         >  3:889445/1‑149 (MQ=255)
                                              gTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTc                                                                                                        <  3:755066/149‑1 (MQ=255)
                                              gTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTc                                                                                                        <  3:296027/149‑1 (MQ=255)
                                                 aaCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCa                                                                                                                                 <  3:915571/121‑1 (MQ=255)
                                                     aCTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacga                                                                                                 >  3:378794/1‑149 (MQ=255)
                                                             tACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGAt                                                                                         >  4:399895/1‑149 (MQ=255)
                                                             tACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGAt                                                                                         <  4:344986/149‑1 (MQ=255)
                                                             tACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGAt                                                                                         >  2:655150/1‑149 (MQ=255)
                                                              aCGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgac                                                                                                <  2:604779/141‑1 (MQ=255)
                                                              aCGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATa                                                                                        <  1:715020/149‑1 (MQ=255)
                                                              aCGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTCTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgac                                                                                                >  1:604779/1‑141 (MQ=255)
                                                                  ggCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCt                                                                                    <  4:415470/149‑1 (MQ=255)
                                                                  ggCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCt                                                                                    >  3:276912/1‑149 (MQ=255)
                                                                     gacgacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTc                                                                                 >  2:179458/1‑149 (MQ=255)
                                                                         acCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTg                                                                             >  4:671056/1‑149 (MQ=255)
                                                                           ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                           >  3:766951/1‑149 (MQ=255)
                                                                           ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                           >  4:517175/1‑149 (MQ=255)
                                                                           ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                           <  3:169445/149‑1 (MQ=255)
                                                                             ggAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGt                                                                         >  1:303205/1‑149 (MQ=255)
                                                                                 aTGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTaa                                                                     <  4:617143/149‑1 (MQ=255)
                                                                                  tGTCGTTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAAc                                                                    <  3:455292/149‑1 (MQ=255)
                                                                                  tGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAAc                                                                    >  1:280615/1‑149 (MQ=255)
                                                                                      ggTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTc                                                                <  4:222120/149‑1 (MQ=255)
                                                                                       gTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACTGATACGCTGTCCGTGACGTTTACCTTTcc                                                               >  3:296837/1‑149 (MQ=255)
                                                                                        tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTccc                                                              >  1:757916/1‑149 (MQ=255)
                                                                                         ttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCt                                                             >  2:814999/1‑149 (MQ=255)
                                                                                          tACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCtt                                                            <  3:541341/149‑1 (MQ=255)
                                                                                           aCATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCttt                                                           <  4:889445/149‑1 (MQ=255)
                                                                                             aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                         >  1:74600/1‑149 (MQ=255)
                                                                                                   caacTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCa                                                   >  2:533524/1‑149 (MQ=255)
                                                                                                          cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGc                                                                                                      <  1:387909/91‑1 (MQ=255)
                                                                                                          cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGc                                                                                                      >  2:387909/1‑91 (MQ=255)
                                                                                                          cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACttt                                                                 >  1:10720/1‑128 (MQ=255)
                                                                                                           aTTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTa                                           >  4:595730/1‑149 (MQ=255)
                                                                                                           aTTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTa                                           >  3:798241/1‑149 (MQ=255)
                                                                                                            ttGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTaa                                          >  2:92864/1‑149 (MQ=255)
                                                                                                            ttGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTaa                                          >  4:742904/1‑149 (MQ=255)
                                                                                                            ttGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGCTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTaa                                          <  3:671056/149‑1 (MQ=255)
                                                                                                             tGTATTGCTCCTCCGACTTACTCAGCGCCGACAACGAAGTCCAGATTCTAGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGCCGCTACGTCGACCGATACGCTGTCCGTGACGTTTAACTTTCCCATTTAGTACGAGGGTGTTAAc                                         <  4:651666/149‑1 (MQ=255)
                                                                                                             tGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAc                                         <  3:814860/149‑1 (MQ=255)
                                                                                                                aTTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgac                                      >  3:750300/1‑149 (MQ=255)
                                                                                                                  tACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgactt                                    <  1:110558/149‑1 (MQ=255)
                                                                                                                   aCTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttc                                   <  1:759183/149‑1 (MQ=255)
                                                                                                                     tcctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGc                                                                                                      >  2:722388/1‑80 (MQ=255)
                                                                                                                     tcctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGc                                                                                                      <  1:722388/80‑1 (MQ=255)
                                                                                                                     tcctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcga                                 >  4:345546/1‑149 (MQ=255)
                                                                                                                     tcctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcga                                 <  4:793919/149‑1 (MQ=255)
                                                                                                                      cctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgac                                <  3:345546/149‑1 (MQ=255)
                                                                                                                        tccGACTTACTCAGCGCCGCCAACGAAGTCCAGCTTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGTCGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTCAcgacttcgactt                              >  2:41131/1‑149 (MQ=255)
                                                                                                                        tccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTAcgacga                                                                                              <  4:248136/85‑1 (MQ=255)
                                                                                                                        tccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTAcgacga                                                                                              >  3:248136/1‑85 (MQ=255)
                                                                                                                        tccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgactt                              >  1:555150/1‑149 (MQ=255)
                                                                                                                        tccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgactt                              <  1:411153/149‑1 (MQ=255)
                                                                                                                           gACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcga                           <  2:280615/149‑1 (MQ=255)
                                                                                                                           gACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcga                           <  4:769152/149‑1 (MQ=255)
                                                                                                                             cTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgact                         <  2:827795/149‑1 (MQ=255)
                                                                                                                             cTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgact                         <  4:174716/149‑1 (MQ=255)
                                                                                                                              ttACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgactt                        >  3:61846/1‑149 (MQ=255)
                                                                                                                                  tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                    <  2:495538/149‑1 (MQ=255)
                                                                                                                                  tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                    >  1:545603/1‑149 (MQ=255)
                                                                                                                                     gcgcCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAg                 <  3:91183/149‑1 (MQ=255)
                                                                                                                                      cgccgcCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGt                <  4:296837/149‑1 (MQ=255)
                                                                                                                                         cgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc             <  1:147659/149‑1 (MQ=255)
                                                                                                                                         cgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc             <  1:179458/149‑1 (MQ=255)
                                                                                                                                         cgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc             <  1:672193/149‑1 (MQ=255)
                                                                                                                                         cctgaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc             >  2:566324/8‑149 (MQ=255)
                                                                                                                                          gccAACGAAGTCCAGATTCTGGCCTTCTTTCAAGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt            <  1:285800/149‑1 (MQ=255)
                                                                                                                                           ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg           >  1:778280/1‑149 (MQ=255)
                                                                                                                                            cAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGc          >  3:388430/1‑149 (MQ=255)
                                                                                                                                            cAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGc          <  3:395843/149‑1 (MQ=255)
                                                                                                                                             aaCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGca         >  1:249479/1‑149 (MQ=255)
                                                                                                                                             aaCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGca         >  3:509612/1‑149 (MQ=255)
                                                                                                                                              aCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcac        >  3:862089/1‑149 (MQ=255)
                                                                                                                                                gAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcag      >  2:403587/1‑149 (MQ=255)
                                                                                                                                                gAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttc                             >  1:855532/1‑126 (MQ=255)
                                                                                                                                                gAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttc                             <  2:855532/126‑1 (MQ=255)
                                                                                                                                                 aaGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt            <  2:393524/142‑1 (MQ=255)
                                                                                                                                                 aaGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt            >  1:393524/1‑142 (MQ=255)
                                                                                                                                                   gTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcagcag   >  3:917373/1‑149 (MQ=255)
                                                                                                                                                    tCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACagcagc  >  2:540778/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTTTTCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGC  >  NC_000913/3451256‑3451552

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