Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,769,777 T→G 7.6% V235G (GTA→GGA)  ydiK → UPF0118 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,7770TG 7.6% 267.0 / 13.4 79V235G (GTA→GGA) ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base G (4/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

AAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGT  >  NC_000913/1769630‑1769923
                                                                                                                                                   |                                                                                                                                                  
aaGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTa                                                                                                                                                   >  3:525285/1‑149 (MQ=255)
 aGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTa                                                                                                                                                   >  2:287435/1‑148 (MQ=255)
          ccATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTa                                                                                                                                          >  2:628130/1‑148 (MQ=255)
           cATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTAcc                                                                                                                                        >  4:162145/1‑149 (MQ=255)
            aTTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCACCGTATCCGGCGTACCt                                                                                                                                       >  4:763538/1‑149 (MQ=255)
               ttGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTAt                                                                                                                                    <  4:163688/149‑1 (MQ=255)
                   aaCCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCaa                                                                                                                                <  3:125408/149‑1 (MQ=255)
                       cGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGtt                                                                                                                            <  4:838051/149‑1 (MQ=255)
                          cTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCt                                                                                                                         <  3:162145/149‑1 (MQ=255)
                               aGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACgg                                                                                                                    >  4:600318/1‑149 (MQ=255)
                                  cGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGgtgt                                                                                                                 <  4:198174/149‑1 (MQ=255)
                                    ttCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGtt                                                                                                                <  1:10944/148‑1 (MQ=255)
                                     tCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTa                                                                                                               >  2:210142/1‑148 (MQ=255)
                                        cgGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATg                                                                                                            <  1:820761/148‑1 (MQ=255)
                                          gagaTGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGAt                                                                                                          <  1:210142/148‑1 (MQ=255)
                                          gagaTGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATc                                                                                                         <  3:600318/149‑1 (MQ=255)
                                          gagaTGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGACCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATc                                                                                                         <  3:703015/149‑1 (MQ=255)
                                               gccgccGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCctct                                                                                                     <  1:509275/148‑1 (MQ=255)
                                                   ccgTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGc                                                                                                <  4:891469/149‑1 (MQ=255)
                                                      gcagcacggcggcACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCtt                                                                                             >  4:402052/8‑149 (MQ=255)
                                                                  cacaGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTgg                                                                                  <  1:287435/148‑1 (MQ=255)
                                                                   acaGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGc                                                                                 <  1:628130/148‑1 (MQ=255)
                                                                    caGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGcc                                                                                >  2:859845/1‑148 (MQ=255)
                                                                      ggCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGgtgt                                                                                                                 >  4:193412/1‑113 (MQ=255)
                                                                      ggCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGgtgt                                                                                                                 <  3:193412/113‑1 (MQ=255)
                                                                      ggCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCg                                                                              >  2:666591/1‑148 (MQ=255)
                                                                        cTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                                                                              <  1:129699/114‑1 (MQ=255)
                                                                        cTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                                                                              >  2:129699/1‑114 (MQ=255)
                                                                           tCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTc                                                                                          >  4:815207/1‑131 (MQ=255)
                                                                           tCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTc                                                                                          <  3:815207/131‑1 (MQ=255)
                                                                           tCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGcc                                                                         >  2:866882/1‑148 (MQ=255)
                                                                            ccGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCgg                                                                       <  3:238232/149‑1 (MQ=255)
                                                                            ccGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCg                                                                        >  1:708799/1‑148 (MQ=255)
                                                                              gcgcgGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTa                                                                     <  4:234354/149‑1 (MQ=255)
                                                                              gcgcgGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTa                                                                     >  3:48942/1‑149 (MQ=255)
                                                                                  ggTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGa                                                                 <  3:159775/149‑1 (MQ=255)
                                                                                          tGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCg                                                         >  4:575377/1‑149 (MQ=255)
                                                                                                   tggtgACGGTGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTgg                                                <  3:763538/149‑1 (MQ=255)
                                                                                                   tggtgACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTgg                                                >  3:153288/1‑149 (MQ=255)
                                                                                                   tggtgACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTgg                                                <  3:851300/149‑1 (MQ=255)
                                                                                                     gtgACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGc                                               <  1:859845/148‑1 (MQ=255)
                                                                                                            cGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                                                                              >  1:396298/1‑78 (MQ=255)
                                                                                                            cGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                                                                              <  2:396298/78‑1 (MQ=255)
                                                                                                             gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCt                                                                                     >  4:858202/1‑102 (MQ=255)
                                                                                                             gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCt                                                                                     <  3:858202/102‑1 (MQ=255)
                                                                                                             gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTactg                                        <  4:269427/147‑1 (MQ=255)
                                                                                                             gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTactg                                        >  3:269427/1‑147 (MQ=255)
                                                                                                              ttAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTactgga                                      >  2:814434/1‑148 (MQ=255)
                                                                                                                 gTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAAt                                                                                                             <  4:404949/74‑1 (MQ=255)
                                                                                                                 gTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAAt                                                                                                             >  3:404949/1‑74 (MQ=255)
                                                                                                                  tACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGc                                 >  4:675805/1‑149 (MQ=255)
                                                                                                                       gCAGTGCTTGGCGGTATCGGCCTCGCCGCATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATac                            >  3:499094/1‑149 (MQ=255)
                                                                                                                           tGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATaccac                         >  1:299687/1‑148 (MQ=255)
                                                                                                                              ttGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGa                               <  4:613869/139‑1 (MQ=255)
                                                                                                                              ttGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGa                               >  3:613869/1‑139 (MQ=255)
                                                                                                                              ttGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTgg                     >  3:851040/1‑149 (MQ=255)
                                                                                                                              ttGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTgg                     >  3:23341/1‑149 (MQ=255)
                                                                                                                                       tCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGt             >  2:463324/1‑148 (MQ=255)
                                                                                                                                       tCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTa            <  3:290111/149‑1 (MQ=255)
                                                                                                                                        cGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTAt           >  4:414062/1‑149 (MQ=255)
                                                                                                                                         ggCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTAt           <  1:666591/148‑1 (MQ=255)
                                                                                                                                         ggCCTCGCCGTATCCGGCGTACCTTATGAAAGGTTGCTAACGGTATCAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCCGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTAtt          >  4:555368/1‑149 (MQ=255)
                                                                                                                                          gCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTg         >  3:132285/1‑149 (MQ=255)
                                                                                                                                           gtacGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATaccacc                        >  1:904233/4‑133 (MQ=255)
                                                                                                                                           gtacGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATaccacc                        <  2:904233/130‑1 (MQ=255)
                                                                                                                                           gtacGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTg         >  1:100463/4‑148 (MQ=255)
                                                                                                                                            tacGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGt                                                                             >  2:789417/3‑79 (MQ=255)
                                                                                                                                            tacGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGt                                                                             <  1:789417/77‑1 (MQ=255)
                                                                                                                                            cTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                                                                                   <  4:403715/73‑1 (MQ=255)
                                                                                                                                            cTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                                                                                   >  3:403715/1‑73 (MQ=255)
                                                                                                                                             tCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGtt       >  1:224339/1‑148 (MQ=255)
                                                                                                                                             tCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGtt       >  2:13894/1‑148 (MQ=255)
                                                                                                                                             tCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTa      >  4:182171/1‑149 (MQ=255)
                                                                                                                                             acGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTa      >  3:824757/2‑149 (MQ=255)
                                                                                                                                               gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGc                                                          >  4:510862/1‑95 (MQ=255)
                                                                                                                                               gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGc                                                          <  3:510862/95‑1 (MQ=255)
                                                                                                                                               gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAgt    <  4:499094/149‑1 (MQ=255)
                                                                                                                                               gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAgt    <  4:525285/149‑1 (MQ=255)
                                                                                                                                                 cGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAgtgt  <  3:575377/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGT  >  NC_000913/1769630‑1769923

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