Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,166,155 G→C 5.7% Q522E (CAG→GAG) ‡ ygiS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,1550GC 5.7% 375.0 / 11.9 107Q522E (CAG→GAG) ‡ygiSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (54/46);  new base C (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 = 3.14e-01

Only 100 of 113 total aligned reads displayed.

CAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCT  >  NC_000913/3166007‑3166301
                                                                                                                                                    |                                                                                                                                                  
cAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTg                                                                                                                                                    <  1:135293/149‑1 (MQ=255)
 aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  2:14693/149‑1 (MQ=255)
 aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   >  2:787879/1‑149 (MQ=255)
 aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  2:792051/149‑1 (MQ=255)
 aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  1:414070/149‑1 (MQ=255)
 aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                   <  4:407042/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)
                          tacATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCGGGGGATTATGCAGCGGAAAACCGCCaa                                                                                                                          >  3:550900/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)
                                   tGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGtt                                                                                                                 >  2:378938/1‑149 (MQ=255)
                                     ccGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTa                                                                                                               >  2:171839/1‑149 (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:63667/1‑149 (MQ=255)
                                                     gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               >  3:132124/1‑149 (MQ=255)
                                                     gagaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACg                                                                                               <  1:258302/149‑1 (MQ=255)
                                                          cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  4:132124/149‑1 (MQ=255)
                                                          cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGa                                                                                          <  1:171839/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)
                                                                       gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         >  2:775987/1‑87 (MQ=255)
                                                                       gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         <  1:209622/89‑3 (MQ=255)
                                                                       gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         <  1:775987/89‑3 (MQ=255)
                                                                       gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGGCATAATCCCCAGGATTATgtc                                                                                                                                         >  2:209622/1‑87 (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)
                                                                                       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                                                      <  1:748184/149‑1 (MQ=255)
                                                                                              gATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                      >  4:524043/1‑149 (MQ=255)
                                                                                               aTAAAGGCTTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatga                                                     >  3:738450/1‑149 (MQ=255)
                                                                                               aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  4:160653/106‑1 (MQ=255)
                                                                                               aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                <  1:343900/106‑1 (MQ=255)
                                                                                               aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                >  2:343900/1‑106 (MQ=255)
                                                                                               aTAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAAc                                                                                                >  3:160653/1‑106 (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)
                                                                                                       aTCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTgatgattgatcgg                                             >  1:283848/1‑142 (MQ=255)
                                                                                                          aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                  <  1:574663/125‑1 (MQ=255)
                                                                                                          aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                  >  2:574663/1‑125 (MQ=255)
                                                                                                           aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                     >  4:772677/1‑73 (MQ=255)
                                                                                                           aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcct                                         <  4:314536/149‑1 (MQ=255)
                                                                                                             gtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       <  1:261260/149‑1 (MQ=255)
                                                                                                                 cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   >  3:719059/1‑149 (MQ=255)
                                                                                                                 cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                   >  1:763078/1‑149 (MQ=255)
                                                                                                                  ttGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATa                                  >  1:604859/1‑149 (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                                       >  2:76557/1‑135 (MQ=255)
                                                                                                                           aaCTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcctgc                                       <  1:76557/135‑1 (MQ=255)
                                                                                                                                tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  2:223497/1‑149 (MQ=255)
                                                                                                                                tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  1:293072/1‑149 (MQ=255)
                                                                                                                                tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  3:431654/1‑149 (MQ=255)
                                                                                                                                tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                    >  1:14916/1‑149 (MQ=255)
                                                                                                                                     tGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTc               <  2:174800/149‑1 (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: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)
                                                                                                                                                  cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACg                      >  4:643222/1‑129 (MQ=255)
                                                                                                                                                  cTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACg                      <  3:643222/129‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCT  >  NC_000913/3166007‑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: