Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,166,154 T→C 5.4% Q522R (CAG→CGG) ‡ ygiS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,1540TC 5.4% 362.2 / 7.9 113Q522R (CAG→CGG) ‡ygiSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (54/52);  new base C (3/3);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 116 total aligned reads displayed.

TCAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCT  >  NC_000913/3166006‑3166301
                                                                                                                                                    |                                                                                                                                                   
tCAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCt                                                                                                                                                     <  2:5380/149‑1 (MQ=255)
 cAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTg                                                                                                                                                    <  1:135293/149‑1 (MQ=255)
  aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   >  2:787879/1‑149 (MQ=255)
  aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  2:792051/149‑1 (MQ=255)
  aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  2:14693/149‑1 (MQ=255)
  aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  1:414070/149‑1 (MQ=255)
           atatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATTGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCCCAGGATTATgt                                                                                                                                          <  3:583904/149‑2 (MQ=255)
                   ggAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa                                                                                                                                  >  2:327286/1‑149 (MQ=255)
                      gaTACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTATTTGCTGTAGACATAATCCTGGAGATTATGCAGCGGAAAAcc                                                                                                                               <  3:471852/149‑1 (MQ=255)
                          ctacATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCa                                                                                                                           >  3:696576/1‑149 (MQ=255)
                           tacATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCATCGGAAAACCGCCaa                                                                                                                          >  4:142141/1‑149 (MQ=255)
                           tacATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCGGGGGATTATGCAGCGGAAAACCGCCaa                                                                                                                          >  3:550900/1‑149 (MQ=255)
                            acATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCt                                                                                                                                                     >  2:471698/1‑121 (MQ=255)
                            acATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCt                                                                                                                                                     <  1:471698/121‑1 (MQ=255)
                             cATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACg                                                                                                                        >  4:543714/1‑149 (MQ=255)
                                 aTCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAgg                                                                                                                    <  4:524303/149‑1 (MQ=255)
                                 aTCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAgg                                                                                                                    >  2:146636/1‑149 (MQ=255)
                                  tCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                   >  1:174800/1‑149 (MQ=255)
                                    tGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGtt                                                                                                                 >  2:378938/1‑149 (MQ=255)
                                      ccGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTa                                                                                                               >  2:171839/1‑149 (MQ=255)
                                       cgcgTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTaa                                                                                                              <  2:813805/149‑1 (MQ=255)
                                             gAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGtt                                                                                                        <  4:112492/149‑1 (MQ=255)
                                             gAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATACAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGtt                                                                                                        <  3:479052/149‑1 (MQ=255)
                                                tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGa                                                                                                     <  4:70633/149‑1 (MQ=255)
                                                tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGa                                                                                                     <  4:696576/149‑1 (MQ=255)
                                                    aGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTaa                                                                                                 >  3:749932/1‑149 (MQ=255)
                                                     ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                   >  3:409318/1‑130 (MQ=255)
                                                     ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  1:146636/149‑1 (MQ=255)
                                                     ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACGCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                   <  4:409318/130‑1 (MQ=255)
                                                      gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               <  1:258302/149‑1 (MQ=255)
                                                      gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               >  1:63667/1‑149 (MQ=255)
                                                           cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  1:171839/149‑1 (MQ=255)
                                                           cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  4:132124/149‑1 (MQ=255)
                                                           cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGAGATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  4:550900/149‑1 (MQ=255)
                                                                   aTAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGAt                                                                                  >  1:275355/1‑149 (MQ=255)
                                                                    tAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATa                                                                                 <  1:787879/149‑1 (MQ=255)
                                                                    tAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATa                                                                                 <  4:688156/149‑1 (MQ=255)
                                                                     aGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAg                                                                                <  1:378938/149‑1 (MQ=255)
                                                                        gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         >  2:775987/1‑87 (MQ=255)
                                                                        gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         >  2:209622/1‑87 (MQ=255)
                                                                        gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         <  1:775987/89‑3 (MQ=255)
                                                                        gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGaa                                                                             >  3:314536/1‑149 (MQ=255)
                                                                              ggACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                       <  3:330351/149‑1 (MQ=255)
                                                                              ggACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTAGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                       <  3:543714/149‑1 (MQ=255)
                                                                               gACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCgg                                                                      >  2:674830/1‑149 (MQ=255)
                                                                                    gCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCt                                                                 >  2:261260/1‑149 (MQ=255)
                                                                                        cGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttg                                                             <  3:142141/149‑1 (MQ=255)
                                                                                         gACGGCGATAAAGGCATCAATGTGCCTTGATATACTACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttgg                                                            >  1:100771/1‑149 (MQ=255)
                                                                                           cGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATGAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttggtt                                                          <  1:53236/149‑1 (MQ=255)
                                                                                            ggCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttggttg                                                         >  4:347718/1‑149 (MQ=255)
                                                                                              cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                       <  4:353195/149‑1 (MQ=255)
                                                                                               gATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                      <  4:792432/149‑1 (MQ=255)
                                                                                               gATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                      >  4:524043/1‑149 (MQ=255)
                                                                                                aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  1:343900/106‑1 (MQ=255)
                                                                                                aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  4:160653/106‑1 (MQ=255)
                                                                                                aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                >  3:160653/1‑106 (MQ=255)
                                                                                                aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                >  2:343900/1‑106 (MQ=255)
                                                                                                 tAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgat                                                    >  1:647897/1‑149 (MQ=255)
                                                                                                 tAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgat                                                    <  3:336691/149‑1 (MQ=255)
                                                                                                   aaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCa                                                  <  2:275355/149‑1 (MQ=255)
                                                                                                     ggCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     >  2:90370/1‑80 (MQ=255)
                                                                                                     ggCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     <  1:90370/80‑1 (MQ=255)
                                                                                                     ggCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGtt                                                                                                        >  1:136067/1‑93 (MQ=255)
                                                                                                     ggCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACt                                                <  4:741035/149‑1 (MQ=255)
                                                                                                      gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttgg                                                            <  4:398205/136‑1 (MQ=255)
                                                                                                      gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttgg                                                            >  3:398205/1‑136 (MQ=255)
                                                                                                      gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACtt                                               <  1:737196/149‑1 (MQ=255)
                                                                                                        aTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgattgatcgg                                             >  1:283848/1‑142 (MQ=255)
                                                                                                        aTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTct                                             >  1:210438/1‑149 (MQ=255)
                                                                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                  >  2:574663/1‑125 (MQ=255)
                                                                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                  <  1:574663/125‑1 (MQ=255)
                                                                                                            aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     >  4:772677/1‑73 (MQ=255)
                                                                                                            aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     <  3:772677/73‑1 (MQ=255)
                                                                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       <  1:261260/149‑1 (MQ=255)
                                                                                                                  cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   >  1:763078/1‑149 (MQ=255)
                                                                                                                  cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   >  3:719059/1‑149 (MQ=255)
                                                                                                                   ttGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATa                                  >  1:604859/1‑149 (MQ=255)
                                                                                                                     gATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                       >  3:5965/1‑110 (MQ=255)
                                                                                                                     gATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                       <  4:5965/110‑1 (MQ=255)
                                                                                                                         taCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGa                                                                                   >  4:401083/1‑94 (MQ=255)
                                                                                                                          aCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCa                           >  1:792475/1‑149 (MQ=255)
                                                                                                                           cAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCATCGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCaa                          >  3:67736/1‑148 (MQ=255)
                                                                                                                            aaCTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       >  2:76557/1‑135 (MQ=255)
                                                                                                                            aaCTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       <  1:76557/135‑1 (MQ=255)
                                                                                                                                 tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  1:14916/1‑149 (MQ=255)
                                                                                                                                 tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  3:431654/1‑149 (MQ=255)
                                                                                                                                 tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  1:293072/1‑149 (MQ=255)
                                                                                                                                   tGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCttt                  >  2:116414/1‑149 (MQ=255)
                                                                                                                                      tGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTc               <  2:174800/149‑1 (MQ=255)
                                                                                                                                       gTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCg              >  4:383089/1‑149 (MQ=255)
                                                                                                                                          tgcATAATCCCCAGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCAt           >  2:201883/3‑149 (MQ=255)
                                                                                                                                          gACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCAt           <  2:792475/149‑1 (MQ=255)
                                                                                                                                          gACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCAt           <  2:63667/149‑1 (MQ=255)
                                                                                                                                           aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATc          >  1:363848/1‑149 (MQ=255)
                                                                                                                                           aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATc          >  3:753085/1‑149 (MQ=255)
                                                                                                                                            cATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCa         >  4:381220/1‑149 (MQ=255)
                                                                                                                                                aTCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGa     >  2:683505/1‑149 (MQ=255)
                                                                                                                                                 tCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGAt    >  1:451981/1‑149 (MQ=255)
                                                                                                                                                   cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTGTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCt  >  2:681744/1‑149 (MQ=255)
                                                                                                                                                   cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACg                      <  3:643222/129‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCT  >  NC_000913/3166006‑3166301

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