Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,166,156 G→A 5.6% P521P (CCC→CCT ygiS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,1560GA 5.6% 355.4 / 7.6 107P521P (CCC→CCTygiSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (54/47);  new base A (3/3);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

Only 100 of 112 total aligned reads displayed.

AAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCT  >  NC_000913/3166008‑3166301
                                                                                                                                                    |                                                                                                                                                 
aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  2:14693/149‑1 (MQ=255)
aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  2:792051/149‑1 (MQ=255)
aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   >  2:787879/1‑149 (MQ=255)
aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  4:407042/149‑1 (MQ=255)
aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  1:414070/149‑1 (MQ=255)
         atatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATTGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCCCAGGATTATgt                                                                                                                                          <  3:583904/149‑2 (MQ=255)
            ttGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGc                                                                                                                                       <  2:72551/149‑1 (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)
                           cATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACg                                                                                                                        >  4:543714/1‑149 (MQ=255)
                               aTCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAgg                                                                                                                    >  2:146636/1‑149 (MQ=255)
                               aTCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAgg                                                                                                                    <  4:524303/149‑1 (MQ=255)
                                tCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                   >  1:174800/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)
                                                   ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                   >  3:409318/1‑130 (MQ=255)
                                                   ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  1:146636/149‑1 (MQ=255)
                                                   ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACGCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                   <  4:409318/130‑1 (MQ=255)
                                                    gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               >  3:132124/1‑149 (MQ=255)
                                                    gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               <  1:258302/149‑1 (MQ=255)
                                                    gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               >  1:63667/1‑149 (MQ=255)
                                                         cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  4:132124/149‑1 (MQ=255)
                                                         cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  1:171839/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                                                                                                                                         <  1:775987/89‑3 (MQ=255)
                                                                      gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         <  1:209622/89‑3 (MQ=255)
                                                                      gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         >  2:209622/1‑87 (MQ=255)
                                                                      gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGaa                                                                             >  3:314536/1‑149 (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:524043/1‑149 (MQ=255)
                                                                                             gATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                      <  4:792432/149‑1 (MQ=255)
                                                                                             gATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                      <  1:748184/149‑1 (MQ=255)
                                                                                              aTAAAGGCTTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatga                                                     >  3:738450/1‑149 (MQ=255)
                                                                                              aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  4:160653/106‑1 (MQ=255)
                                                                                              aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                >  3:160653/1‑106 (MQ=255)
                                                                                              aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                >  2:343900/1‑106 (MQ=255)
                                                                                              aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  1:343900/106‑1 (MQ=255)
                                                                                              aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatga                                                     >  3:56364/1‑149 (MQ=255)
                                                                                               tAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgat                                                    <  3:336691/149‑1 (MQ=255)
                                                                                               tAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgat                                                    >  1:647897/1‑149 (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)
                                                                                                   ggCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGtt                                                                                                        <  2:136067/93‑1 (MQ=255)
                                                                                                   ggCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACt                                                <  4:741035/149‑1 (MQ=255)
                                                                                                    gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttgg                                                            >  3:398205/1‑136 (MQ=255)
                                                                                                    gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttgg                                                            <  4:398205/136‑1 (MQ=255)
                                                                                                    gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACtt                                               <  1:737196/149‑1 (MQ=255)
                                                                                                      aTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgattgatcgg                                             >  1:283848/1‑142 (MQ=255)
                                                                                                      aTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTct                                             >  1:210438/1‑149 (MQ=255)
                                                                                                         aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                  <  1:574663/125‑1 (MQ=255)
                                                                                                         aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                  >  2:574663/1‑125 (MQ=255)
                                                                                                          aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     >  4:772677/1‑73 (MQ=255)
                                                                                                          aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     <  3:772677/73‑1 (MQ=255)
                                                                                                          aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcct                                         <  4:314536/149‑1 (MQ=255)
                                                                                                                cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   >  3:719059/1‑149 (MQ=255)
                                                                                                                 ttGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATa                                  >  1:604859/1‑149 (MQ=255)
                                                                                                                   gATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                       <  4:5965/110‑1 (MQ=255)
                                                                                                                   gATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                       >  3:5965/1‑110 (MQ=255)
                                                                                                                       taCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGa                                                                                   <  3:401083/94‑1 (MQ=255)
                                                                                                                       taCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGa                                                                                   >  4:401083/1‑94 (MQ=255)
                                                                                                                        aCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCa                           >  1:792475/1‑149 (MQ=255)
                                                                                                                        aCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCa                           <  1:327286/149‑1 (MQ=255)
                                                                                                                         cAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCATCGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCaa                          >  3:67736/1‑148 (MQ=255)
                                                                                                                          aaCTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       <  1:76557/135‑1 (MQ=255)
                                                                                                                          aaCTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       >  2:76557/1‑135 (MQ=255)
                                                                                                                               tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  1:14916/1‑149 (MQ=255)
                                                                                                                               tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  1:293072/1‑149 (MQ=255)
                                                                                                                               tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  2:223497/1‑149 (MQ=255)
                                                                                                                               tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  3:431654/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)
                                                                                                                                        gACATAATCCTGGGGATTATGCGGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   >  1:493156/1‑125 (MQ=255)
                                                                                                                                        gACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   <  2:493156/125‑1 (MQ=255)
                                                                                                                                        gACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCAt           <  2:63667/149‑1 (MQ=255)
                                                                                                                                        gACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCAt           <  2:792475/149‑1 (MQ=255)
                                                                                                                                         aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATc          >  1:363848/1‑149 (MQ=255)
                                                                                                                                         aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATc          >  3:753085/1‑149 (MQ=255)
                                                                                                                                              aTCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGa     >  2:683505/1‑149 (MQ=255)
                                                                                                                                                 cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTGTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCt  >  2:681744/1‑149 (MQ=255)
                                                                                                                                                 cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACg                      >  4:643222/1‑129 (MQ=255)
                                                                                                                                                 cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACg                      <  3:643222/129‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCT  >  NC_000913/3166008‑3166301

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