Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,451,416 C→A 8.1% G90V (GGC→GTC)  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,4160CA 8.1% 437.0 / 23.4 123G90V (GGC→GTC) rplW50S ribosomal subunit protein L23
Reads supporting (aligned to +/- strand):  ref base C (56/57);  new base A (5/5);  total (61/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 = 6.60e-01

Only 100 of 123 total aligned reads displayed.

AAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTGATTTC  >  NC_000913/3451275‑3451561
                                                                                                                                             |                                                                                                                                                 
aaaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCcttctt                                                                                                                                             >  1:925828/1‑148 (MQ=255)
            cTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCcttct                                                                                                                                              <  1:906386/135‑1 (MQ=255)
                tgCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGc                                                                                                                             <  2:925828/148‑1 (MQ=255)
                      tCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCttt                                                                                                                          <  4:114108/145‑1 (MQ=255)
                      tCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCttt                                                                                                                          >  3:114108/1‑145 (MQ=255)
                       cAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAg                                                                                                                              >  4:69333/1‑140 (MQ=255)
                       cAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGATGACGTAAg                                                                                                                              <  3:69333/140‑1 (MQ=255)
                          ggTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAg                                                                                                                   <  2:502473/148‑1 (MQ=255)
                               aCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTcc                                                                                                                     >  1:523627/1‑141 (MQ=255)
                                ccACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                           <  3:383886/118‑1 (MQ=255)
                                    tttAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgac                                                                                                         >  1:84727/1‑148 (MQ=255)
                                      tAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTAcgacga                                                                                                       <  2:84727/148‑1 (MQ=255)
                                              tGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGAt                                                                                                  <  2:413824/145‑1 (MQ=255)
                                              tGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGAt                                                                                                  >  1:413824/1‑145 (MQ=255)
                                               ggCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGc                                                                                              >  3:210047/1‑148 (MQ=255)
                                                gCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                             >  2:140909/1‑132 (MQ=255)
                                                gCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                             <  1:140909/132‑1 (MQ=255)
                                                   acgacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTc                                                                                          <  4:210047/148‑1 (MQ=255)
                                                    cgacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                           <  3:266719/98‑1 (MQ=255)
                                                    cgacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                           >  4:266719/1‑98 (MQ=255)
                                                         cGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                           <  4:1007446/93‑1 (MQ=255)
                                                         cGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                           >  3:1007446/1‑93 (MQ=255)
                                                                     tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCa                                                                                                                                          >  4:551569/1‑82 (MQ=255)
                                                                     tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCa                                                                                                                                          <  3:551569/82‑1 (MQ=255)
                                                                     tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGt                                                                                                                  <  2:604725/106‑1 (MQ=255)
                                                                     tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGt                                                                                                                  >  1:604725/1‑106 (MQ=255)
                                                                      ttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCACCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTcc                                                                                                                     <  1:534710/102‑1 (MQ=255)
                                                                      ttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCACCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTcc                                                                                                                     >  2:534710/1‑102 (MQ=255)
                                                                           tttAACAACTGTCAGTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGCCCTTCTTTCAGGGTGACGTAAGCttt                                                                                                                          >  4:364635/1‑92 (MQ=255)
                                                                           tttAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacga                                                                                                          >  4:906170/1‑108 (MQ=255)
                                                                           tttAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCGTTTTTCCAGTCGCTacga                                                                                                          <  3:906170/108‑1 (MQ=255)
                                                                           tttAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGCCCTTCTTTCAGGGTGACGTAAGCttt                                                                                                                          <  3:364635/92‑1 (MQ=255)
                                                                            ttAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACgacgac                                                                                                      <  3:154260/111‑1 (MQ=255)
                                                                            ttAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACgacgac                                                                                                      >  4:154260/1‑111 (MQ=255)
                                                                               acaacTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGa                                                                                                                                    >  1:764554/1‑78 (MQ=255)
                                                                               acaacTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGa                                                                                                                                    <  2:764554/78‑1 (MQ=255)
                                                                                  acTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTAcgacga                                                                                                       <  1:787950/104‑1 (MQ=255)
                                                                                  acTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTAcgacga                                                                                                       >  2:787950/1‑104 (MQ=255)
                                                                                          tGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgacg                                                                                                        >  3:765305/1‑95 (MQ=255)
                                                                                          tGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgacg                                                                                                        <  4:765305/95‑1 (MQ=255)
                                                                                          tGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCt                                                                      <  2:846138/129‑1 (MQ=255)
                                                                                             aTTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAg                                                                                                                   >  1:1083565/1‑81 (MQ=255)
                                                                                             aTTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAg                                                                                                                   <  2:1083565/81‑1 (MQ=255)
                                                                                               tACACCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGt                                                                                  >  4:846491/1‑112 (MQ=255)
                                                                                                 ctcctcCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCg                                                                                        <  2:489233/104‑1 (MQ=255)
                                                                                                 ctcctcCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCg                                                                                        >  1:489233/1‑104 (MQ=255)
                                                                                                  tcctccGACTTACTCAGCGCCGCCAACGAAGGCCAGAATCTGGActtc                                                                                                                                               >  1:805633/1‑48 (MQ=21)
                                                                                                  tcctccGACTTACTCAGCGCCGCCAACGAAGGCCAGAATCTGGActtc                                                                                                                                               <  2:805633/48‑1 (MQ=21)
                                                                                                   cctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgact                                              <  3:440705/144‑1 (MQ=255)
                                                                                                   cctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgact                                              >  4:440705/1‑144 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCg                                                                                                                <  2:696821/75‑1 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCg                                                                                                                >  1:696821/1‑75 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACt                                                                            <  2:854998/111‑1 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACt                                                                            >  1:854998/1‑111 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCtt                                                                     <  1:268132/118‑1 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgactt                                       >  1:672066/1‑148 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgac                                         <  1:74533/146‑1 (MQ=255)
                                                                                                      ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgac                                         >  2:74533/1‑146 (MQ=255)
                                                                                                       cGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCg                                                                                                                <  4:599807/74‑1 (MQ=255)
                                                                                                       cGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCg                                                                                                                >  3:599807/1‑74 (MQ=255)
                                                                                                        gACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcg                                     >  3:907170/1‑148 (MQ=255)
                                                                                                          cTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgac                                   >  2:603164/1‑148 (MQ=255)
                                                                                                            tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                             <  2:471255/72‑1 (MQ=255)
                                                                                                            tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                             >  1:471255/1‑72 (MQ=255)
                                                                                                            tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACttt                                                                          >  1:348915/1‑107 (MQ=255)
                                                                                                            tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACttt                                                                          <  2:348915/107‑1 (MQ=255)
                                                                                                            tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTc                                                                         <  2:446508/108‑1 (MQ=255)
                                                                                                              cTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCg                                                                                                    >  3:257469/1‑79 (MQ=255)
                                                                                                              cTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCg                                                                                                    <  4:257469/79‑1 (MQ=255)
                                                                                                               tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaa                              >  3:877260/1‑148 (MQ=255)
                                                                                                               tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaa                              <  4:907170/148‑1 (MQ=255)
                                                                                                               tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaa                              <  2:672066/148‑1 (MQ=255)
                                                                                                                  gcgcCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTcc                                                                        >  1:1125918/1‑103 (MQ=255)
                                                                                                                  gcgcCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTcc                                                                        <  2:1125918/103‑1 (MQ=255)
                                                                                                                   cgccgcCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgacg                                                                                                        <  4:403389/70‑1 (MQ=255)
                                                                                                                       gccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGttt                       >  2:931427/1‑147 (MQ=255)
                                                                                                                       gccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGttt                       <  1:931427/147‑1 (MQ=255)
                                                                                                                         gaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaa                              <  1:219606/134‑1 (MQ=255)
                                                                                                                         gaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaa                              >  2:219606/5‑138 (MQ=255)
                                                                                                                           aCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                                  <  2:286320/100‑1 (MQ=255)
                                                                                                                           aCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTAc                                                              <  2:992290/104‑1 (MQ=255)
                                                                                                                           aCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAggg                                                         >  1:383808/1‑109 (MQ=255)
                                                                                                                           aCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAggg                                                         <  2:383808/109‑1 (MQ=255)
                                                                                                                            cGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTa                                                                               >  1:692421/1‑86 (MQ=255)
                                                                                                                            cGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcga                                    >  3:834549/1‑129 (MQ=255)
                                                                                                                              aaGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTaa                                                                              >  3:442226/1‑85 (MQ=255)
                                                                                                                              aaGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTaa                                                                              <  4:442226/85‑1 (MQ=255)
                                                                                                                               aggccAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgactt                                       >  3:913969/4‑123 (MQ=255)
                                                                                                                                ggccAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                             <  2:730776/130‑1 (MQ=255)
                                                                                                                                  ccAGATTCTGGCCTTCTTTTAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc                      <  4:613971/137‑1 (MQ=255)
                                                                                                                                   cAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTa                                                                               >  2:453274/1‑79 (MQ=255)
                                                                                                                                   cAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTa                                                                               <  1:453274/79‑1 (MQ=255)
                                                                                                                                   cAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTaa                                                                              >  4:1105670/1‑80 (MQ=255)
                                                                                                                                   cAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTaa                                                                              <  3:1105670/80‑1 (MQ=255)
                                                                                                                                    aGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgac                                               >  1:439260/1‑110 (MQ=255)
                                                                                                                                           ggCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTg       <  4:272969/143‑1 (MQ=255)
                                                                                                                                           ggCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTg       >  3:272969/1‑143 (MQ=255)
                                                                                                                                           ggCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTGATTTc  <  4:877260/148‑1 (MQ=255)
                                                                                                                                            gCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcg                                     <  4:753576/112‑1 (MQ=255)
                                                                                                                                            gCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcg                                     >  3:753576/1‑112 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
AAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTGATTTC  >  NC_000913/3451275‑3451561

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