Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,052,380 C→T 5.2% pseudogene (949/1053 nt) yeeL ← pseudogene, glycosyltransferase homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,052,3800CT 5.2% 259.6 / 3.2 77pseudogene (949/1053 nt)yeeLpseudogene, glycosyltransferase homology
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base T (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

ACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATCGAGGGTGGTAAGGAT  >  NC_000913/2052234‑2052528
                                                                                                                                                  |                                                                                                                                                    
accGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCgg                                                                                                                                                    <  1:543512/149‑1 (MQ=255)
accGTTAATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCgg                                                                                                                                                    <  3:491763/149‑1 (MQ=255)
 ccGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCggg                                                                                                                                                   <  3:537419/149‑1 (MQ=255)
    ttCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACCCTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCa                                                                                                                                                <  1:482823/149‑1 (MQ=255)
          aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATTCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGt                                                                                                                                          >  4:94772/1‑149 (MQ=255)
          aCAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGt                                                                                                                                          <  4:196839/149‑1 (MQ=255)
           cAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGtt                                                                                                                                         <  1:26174/149‑1 (MQ=255)
                      tACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGc                                                                                                                              <  1:433941/149‑1 (MQ=255)
                         tATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATa                                                                                                                           >  3:253595/1‑149 (MQ=255)
                              aTAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTg                                                                                                                      >  3:156880/1‑149 (MQ=255)
                              aTAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTg                                                                                                                      >  4:108143/1‑149 (MQ=255)
                                aaGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGca                                                                                                                    >  1:381851/1‑149 (MQ=255)
                                       gTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAg                                                                                                             <  4:284358/149‑1 (MQ=255)
                                        tctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGTGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAgg                                                                                                            >  4:137615/1‑149 (MQ=255)
                                        tctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGTGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAgg                                                                                                            >  2:118576/1‑149 (MQ=255)
                                        tctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAgg                                                                                                            >  3:285305/1‑149 (MQ=255)
                                        tctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAgg                                                                                                            >  2:251233/1‑149 (MQ=255)
                                         ctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGTGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAggg                                                                                                           <  3:137615/149‑1 (MQ=255)
                                         ctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGTGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAggg                                                                                                           <  1:118576/149‑1 (MQ=255)
                                         ctcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAggg                                                                                                           <  4:371466/149‑1 (MQ=255)
                                          tcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAgggg                                                                                                          >  3:434868/1‑149 (MQ=255)
                                          tcATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAgggg                                                                                                          <  4:264028/149‑1 (MQ=255)
                                              atTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGtatt                                                                                                      >  1:350503/1‑149 (MQ=255)
                                                 ctcGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATtata                                                                                                   <  3:353331/149‑1 (MQ=255)
                                                   cGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAgg                                                                                                 <  2:392990/149‑1 (MQ=255)
                                                        aaaCACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACtt                                                                                            <  4:434868/149‑1 (MQ=255)
                                                              tttttAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCata                                                                                      >  1:534688/1‑149 (MQ=255)
                                                               ttttAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCatat                                                                                     <  1:522428/149‑1 (MQ=255)
                                                                tttAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATg                                                                                    >  4:236604/1‑149 (MQ=255)
                                                                  tAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGtatt                                                                                                      <  2:225666/129‑1 (MQ=255)
                                                                  tAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGtatt                                                                                                      >  1:225666/1‑129 (MQ=255)
                                                                         tttGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGa                                                                           <  2:381851/149‑1 (MQ=255)
                                                                         tttGCGAATTCATCCAGACTTATATTGTCAATTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAg                                                                          >  2:767/1‑149 (MQ=255)
                                                                          ttGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAg                                                                          <  1:251233/149‑1 (MQ=255)
                                                                            gCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGtatta                                                                                                     >  3:209340/1‑120 (MQ=255)
                                                                            gCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGtatta                                                                                                     <  4:209340/120‑1 (MQ=255)
                                                                                   cattcAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGa                                                                 <  3:94772/145‑1 (MQ=255)
                                                                                   cATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGa                                                                 <  3:236604/149‑1 (MQ=255)
                                                                                   cATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGa                                                                 <  2:534688/149‑1 (MQ=255)
                                                                                   cATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGa                                                                 <  2:489016/149‑1 (MQ=255)
                                                                                        aGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTAt                                                            <  2:350503/149‑1 (MQ=255)
                                                                                                 ttGTCAATTTTTTGCCCAACGACGGATGAATCCGAGACATCGGCCATTGCGGGCCTTCGCGTTTTGTACTCAGCATAAACTGCAACAACGGGGGTATTATAGGCAGTTGCCAGATGTAATATTAAGGTGTCAGGAGTTATTGCGAtata                                                   >  4:520973/1‑149 (MQ=255)
                                                                                                  tGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATaa                                                  <  3:277270/149‑1 (MQ=255)
                                                                                                   gTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAg                                                 <  3:108143/149‑1 (MQ=255)
                                                                                                   gTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAg                                                 >  4:533287/1‑149 (MQ=255)
                                                                                                    tCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTaa                                                                                 <  2:122446/116‑1 (MQ=255)
                                                                                                    tCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTaa                                                                                 >  1:122446/1‑116 (MQ=255)
                                                                                                    tCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGc                                                >  3:245584/1‑149 (MQ=255)
                                                                                                    tCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGc                                                >  3:45605/1‑149 (MQ=255)
                                                                                                      aaTTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCAt                                              >  1:326893/1‑149 (MQ=255)
                                                                                                      aaTTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTAAGGAGTTATTGCGATATAAGCAt                                              >  4:68467/1‑149 (MQ=255)
                                                                                                          ttttGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCttt                                          <  4:45605/149‑1 (MQ=255)
                                                                                                            ttGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTa                                        <  2:326893/149‑1 (MQ=255)
                                                                                                               cccAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGa                                     <  1:767/149‑1 (MQ=255)
                                                                                                                ccAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGa                                     <  3:318128/148‑1 (MQ=255)
                                                                                                                ccAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGa                                     >  4:318128/1‑148 (MQ=255)
                                                                                                                ccAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGAt                                    >  2:211341/1‑149 (MQ=255)
                                                                                                                ccAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGAt                                    >  2:494111/1‑149 (MQ=255)
                                                                                                                 cAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGAtt                                   <  3:124569/149‑1 (MQ=255)
                                                                                                                 cAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGAtt                                   >  4:124569/1‑149 (MQ=255)
                                                                                                                     gacgacTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATtgctg                               >  1:26987/1‑149 (MQ=255)
                                                                                                                           tgactCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAatat                                                                              <  2:34871/96‑1 (MQ=255)
                                                                                                                           tgactCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAatat                                                                              >  1:34871/1‑96 (MQ=255)
                                                                                                                           tgactCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAAc                         <  2:50587/149‑1 (MQ=255)
                                                                                                                               tCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGttt                     >  1:476949/1‑149 (MQ=255)
                                                                                                                               tCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGttt                     >  2:54716/1‑149 (MQ=255)
                                                                                                                               tCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGttt                     <  2:252716/149‑1 (MQ=255)
                                                                                                                                 cGAGACGTCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTaa                   <  4:77319/149‑1 (MQ=255)
                                                                                                                                 cGAGACGTCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTaa                   >  3:77319/1‑149 (MQ=255)
                                                                                                                                   agaCATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAg                                                                   <  3:384247/99‑1 (MQ=255)
                                                                                                                                   agaCATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAg                                                                   >  4:384247/1‑99 (MQ=255)
                                                                                                                                   agaCATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATc                 <  2:41580/149‑1 (MQ=255)
                                                                                                                                     aCATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATCGa               <  3:320506/149‑1 (MQ=255)
                                                                                                                                      cATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGa                                                       <  3:175835/108‑1 (MQ=255)
                                                                                                                                      cATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGa                                                       >  4:175835/1‑108 (MQ=255)
                                                                                                                                       aTCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATCGAgg             >  3:424546/1‑149 (MQ=255)
                                                                                                                                             cATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGTGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATCGAGGGTGGTa       >  1:411293/1‑149 (MQ=255)
                                                                                                                                                  cGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATCGAGGGTGGTAAGGAt  <  2:473775/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTCTCATATTCTCGCCAAAACACTTTTTAATGCCTTTGCGAATTCATCCAGACTTATATTGTCAATTTTTTGCCCAACGACGACTGACTCCGAGACATCGGCCATTGCGGGCCATCGCGTTTTGTAATCAGCATAAATTGCAACAACAGGGGTATTATAGGCACTTGCCATATGTAATATTGAGGTGTCAGGAGTTATTGCGATATAAGCATCTTTTATGATTGCTGCTGAACGTTTAATCGAGGGTGGTAAGGAT  >  NC_000913/2052234‑2052528

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