Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,769,780 C→A 7.8% S236Y (TCC→TAC)  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,7800CA 7.8% 265.2 / 10.6 77S236Y (TCC→TAC) ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base A (4/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAG  >  NC_000913/1769640‑1769927
                                                                                                                                            |                                                                                                                                                   
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                                                                                                                  >  2:129699/1‑114 (MQ=255)
                                                              cTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                                                                                  <  1:129699/114‑1 (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                                                                                                                  <  2:396298/78‑1 (MQ=255)
                                                                                                  cGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                                                                                  >  1:396298/1‑78 (MQ=255)
                                                                                                   gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCt                                                                                         <  3:858202/102‑1 (MQ=255)
                                                                                                   gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCt                                                                                         >  4:858202/1‑102 (MQ=255)
                                                                                                   gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTactg                                            >  3:269427/1‑147 (MQ=255)
                                                                                                   gTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTactg                                            <  4:269427/147‑1 (MQ=255)
                                                                                                    ttAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTactgga                                          >  2:814434/1‑148 (MQ=255)
                                                                                                       gTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAAt                                                                                                                 >  3:404949/1‑74 (MQ=255)
                                                                                                       gTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAAt                                                                                                                 <  4:404949/74‑1 (MQ=255)
                                                                                                        tACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGc                                     >  4:675805/1‑149 (MQ=255)
                                                                                                             gCAGTGCTTGGCGGTATCGGCCTCGCCGCATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATac                                >  3:499094/1‑149 (MQ=255)
                                                                                                                 tGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATaccac                             >  1:299687/1‑148 (MQ=255)
                                                                                                                    ttGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGa                                   >  3:613869/1‑139 (MQ=255)
                                                                                                                    ttGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGa                                   <  4:613869/139‑1 (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                                                                                 <  1:789417/77‑1 (MQ=255)
                                                                                                                                  tacGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGt                                                                                 >  2:789417/3‑79 (MQ=255)
                                                                                                                                  cTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                                                                                       >  3:403715/1‑73 (MQ=255)
                                                                                                                                  cTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                                                                                       <  4:403715/73‑1 (MQ=255)
                                                                                                                                   tCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGtt           >  2:13894/1‑148 (MQ=255)
                                                                                                                                   tCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGtt           >  1:224339/1‑148 (MQ=255)
                                                                                                                                   tCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTa          >  4:182171/1‑149 (MQ=255)
                                                                                                                                   acGCCGGATACGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTa          >  3:824757/2‑149 (MQ=255)
                                                                                                                                     gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGc                                                              <  3:510862/95‑1 (MQ=255)
                                                                                                                                     gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGc                                                              >  4:510862/1‑95 (MQ=255)
                                                                                                                                     gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAgt        <  4:499094/149‑1 (MQ=255)
                                                                                                                                     gCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAgt        <  4:525285/149‑1 (MQ=255)
                                                                                                                                       cGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAgtgt      <  3:575377/149‑1 (MQ=255)
                                                                                                                                            ccGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAg  >  2:681034/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAG  >  NC_000913/1769640‑1769927

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: ‑