Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,161,391 G→A 6.0% P1002P (CCG→CCA mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,161,3910GA 6.0% 206.6 / 2.7 66P1002P (CCG→CCAmdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTATGCAGTCGC  >  NC_000913/2161243‑2161533
                                                                                                                                                    |                                                                                                                                              
cGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCg                                                                                                                                                >  4:66015/1‑149 (MQ=255)
 gATTATGATGACTCCCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCgg                                                                                                                                               >  5:170683/1‑149 (MQ=255)
   ttATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             >  3:22773/1‑149 (MQ=255)
      tgatgaCTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGgtggtg                                                                                                                                          >  7:158636/1‑149 (MQ=255)
      tgatgaCTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGGCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtggt                                                                                                                                           <  3:341356/148‑1 (MQ=255)
       gatgaCTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGgtggtg                                                                                                                                          >  5:115309/1‑148 (MQ=255)
        atgaCTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCAGTGGTGTa                                                                                                                                        >  3:51420/1‑149 (MQ=255)
             tACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             <  1:161375/139‑1 (MQ=255)
             tACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             >  2:161375/1‑139 (MQ=255)
               cccTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCACCCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCttt                                                                                                                                 >  1:21835/1‑149 (MQ=255)
                 cTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCttttt                                                                                                                               <  4:22773/149‑1 (MQ=255)
                        cgcTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             >  5:329843/1‑128 (MQ=255)
                        cgcTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             <  6:329843/128‑1 (MQ=255)
                             tttGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGCCTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCg                                                                                                                   >  8:156329/1‑149 (MQ=255)
                                     gcTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             >  6:99018/1‑115 (MQ=255)
                                     gcTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCggtg                                                                                                                                             <  5:99018/115‑1 (MQ=255)
                                     gcTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCAGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGtt                                                                                                           >  7:152409/1‑149 (MQ=255)
                                       tGCCGCTGGTATTGTCGGGCGGCTCCGGCTCGGCGCTGCGGCCCCCCCTGGGGATACCCATTGTCGGCGGCCTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGACGGTGGTGTATCTCTTCTTCGACCGTCTGCGGCTGCGtttt                                                                                                         >  8:116954/1‑149 (MQ=255)
                                        gCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGttttt                                                                                                        >  1:290706/1‑149 (MQ=255)
                                         ccGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTc                                                                                                       <  2:250264/149‑1 (MQ=255)
                                         ccGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTc                                                                                                       <  6:63731/149‑1 (MQ=255)
                                           gcgggTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCCCCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATCCCCCGCCGGGGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTcgc                                                                                                     <  7:116954/146‑1 (MQ=255)
                                                  tggTCGTGCGGCGACGGCTCGGCGCTGCGGCAACCCCTGGGGATCACCAATGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCGGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAAcc                                                                                              <  2:300416/147‑1 (MQ=255)
                                                    gTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCt                                                                                             <  8:158636/148‑1 (MQ=255)
                                                       gggCGTCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTaaa                                                                                          >  4:175078/1‑148 (MQ=255)
                                                           ggcgACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTaaacaaac                                                                                     <  6:115309/149‑1 (MQ=255)
                                                              gACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACgg                                                                                   >  3:176078/1‑148 (MQ=255)
                                                               aCGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGt                                                                                  >  7:377416/1‑148 (MQ=255)
                                                                  gCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCAGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAAcc                                                                              <  4:51420/149‑1 (MQ=255)
                                                                    tCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTc                                                                                                                       <  1:4634/106‑1 (MQ=255)
                                                                    tCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTc                                                                                                                       >  2:4634/1‑106 (MQ=255)
                                                                    tCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAAcc                                                                              >  5:287847/1‑147 (MQ=255)
                                                                      ggAGCTGCGGCAACCCCTGGGGATCACCTTTGTAGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGTTGGTGTATCTCTTTTTCGACCGTCTGTGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGt                                                                          <  1:404910/149‑1 (MQ=255)
                                                                          cTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAAt                                                                      >  1:46704/1‑149 (MQ=255)
                                                                          cTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAAt                                                                      <  6:170683/149‑1 (MQ=255)
                                                                             cGGCAACCCCTGGGGATCACGATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGAc                                                                   <  1:300376/149‑1 (MQ=255)
                                                                             cGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGa                                                                    <  4:247761/148‑1 (MQ=255)
                                                                                 aaCCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGAt                                                               <  3:175078/149‑1 (MQ=255)
                                                                                   ccccTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAGACCTAAACAAACGGTAACCGAGTAAAAGACAGATCt                                                             >  2:127587/1‑149 (MQ=255)
                                                                                   ccccTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAAAGGTAACCGAGTAAATGACAGATCt                                                             >  1:26713/1‑149 (MQ=255)
                                                                                     ccTGGGGATAACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTc                                                           >  6:333614/1‑149 (MQ=255)
                                                                                           gATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGCCCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGa                                                       >  4:77368/1‑147 (MQ=255)
                                                                                           gATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACa                                                     <  8:377416/149‑1 (MQ=255)
                                                                                                     gTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGtt                                            <  2:290706/148‑1 (MQ=255)
                                                                                                     gTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTg                                           >  3:315734/1‑149 (MQ=255)
                                                                                                           cgACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAatt                                     <  2:268660/148‑1 (MQ=255)
                                                                                                             aCTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgt                                   >  5:312456/1‑149 (MQ=255)
                                                                                                             aCTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgt                                   <  3:66015/149‑1 (MQ=255)
                                                                                                              agcgTAAGGAGCTGGCTCATTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtg                                  >  7:157080/4‑149 (MQ=255)
                                                                                                               tGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtg                                  >  5:185135/1‑148 (MQ=255)
                                                                                                               tGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGCTGTATCTCTTTTTCGACCGTCTGCGGCTGCGGTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTACCCGACAGCACCCCTTGGCAattgag                                  >  2:181421/1‑146 (MQ=255)
                                                                                                                cgTAAGGAGCTGGCTCATTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtgga                                >  5:259698/2‑149 (MQ=255)
                                                                                                                 gTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGAGAGATCTTCCCGACAGCACCCGTTGGCAattgtgg                                 <  5:94540/147‑1 (MQ=255)
                                                                                                                  tAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtggatt                              >  4:264101/1‑149 (MQ=255)
                                                                                                                  tAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtggatt                              <  1:336606/149‑1 (MQ=255)
                                                                                                                   aaTGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtggattg                             >  3:73375/1‑149 (MQ=255)
                                                                                                                   aaTGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtggattg                             >  8:283504/1‑149 (MQ=255)
                                                                                                                     tGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtggattgtg                           <  4:176078/149‑1 (MQ=255)
                                                                                                                            gCTCCTTTCGCGGTATACCACGCCGGTGGTGTATCTCTTTTTCTTCCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCg                    <  1:308953/149‑1 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAGACCTAAACAAACGGTAACCGAGTAAAAGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGc                  <  1:127587/147‑1 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTccc                                                         >  2:298415/1‑108 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTccc                                                         <  1:298415/108‑1 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTccc                                                         <  1:245225/108‑1 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTccc                                                         >  2:245225/1‑108 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGctt                >  7:209299/1‑149 (MQ=255)
                                                                                                                                cTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGctt                >  6:337562/1‑149 (MQ=255)
                                                                                                                                    cGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCttt            <  7:156329/149‑1 (MQ=255)
                                                                                                                                    cGCTGTATACCACGCCAGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGcttctt             <  8:152409/148‑1 (MQ=255)
                                                                                                                                     gcggTATACCACGCCGGGGGTGTATCTCTTTTTCGACCGGCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTa           <  4:73375/146‑1 (MQ=255)
                                                                                                                                      cTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTAt          <  5:333614/149‑1 (MQ=255)
                                                                                                                                      cTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTAt          <  8:157080/149‑1 (MQ=255)
                                                                                                                                         tataCCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTATg         <  7:283504/147‑1 (MQ=255)
                                                                                                                                         tataCCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTATGCa       <  5:141992/149‑1 (MQ=255)
                                                                                                                                              cACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTATGCAGTCGc  >  4:283438/1‑149 (MQ=255)
                                                                                                                                              cACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTATGCAGTCGc  >  4:64446/1‑149 (MQ=255)
                                                                                                                                              cACGCCGATGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTACGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCAACCGTTGGCAATTGTGGATTGTGGCTTTAGGCTTCTTTATGAAGTCGc  >  5:11010/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTATGCAGTCGC  >  NC_000913/2161243‑2161533

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: