Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,797 T→G 11.4% intergenic (+29/+19) yphA → / ← yphB DoxX family inner membrane protein/mutarotase superfamily protein, YphB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,673,7970TG11.4% 170.5 / 7.8 73intergenic (+29/+19)yphA/yphBDoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base G (6/3);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCCAGACAGAAGAA  >  NC_000913/2673650‑2673944
                                                                                                                                                   |                                                                                                                                                   
tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                    <  1:625612/149‑1 (MQ=255)
 gACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                   >  1:1006669/1‑149 (MQ=255)
          tGCGGTTGGGCCAAATCTGATTAATTTCTGGAAGAATGTCTGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGTGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                                                                                          <  2:414539/149‑2 (MQ=38)
             ggTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTTCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                                                                       <  1:912094/149‑1 (MQ=255)
              gTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTTCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                                                                      >  2:770350/1‑149 (MQ=255)
                   gCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAtt                                                                                                                                 >  1:352567/1‑149 (MQ=255)
                     cAACTATGATTAATTTCTGGAAGAATGTCAGTATCGCCGGCGCGTTCGTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCTTTAGGATAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                                                                                               <  1:761363/149‑1 (MQ=255)
                      aaaTATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGTATTTAATAAATTACa                                                                                                                              >  1:891409/1‑149 (MQ=255)
                        atatGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                                                            <  1:589393/149‑1 (MQ=255)
                            gATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                                                                                        >  2:287687/1‑149 (MQ=255)
                             aTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActc                                                                                                                       >  2:851006/1‑149 (MQ=255)
                              ttAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActct                                                                                                                      >  2:13810/1‑149 (MQ=255)
                                 aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                                   >  1:1079493/1‑149 (MQ=255)
                                 aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                                   >  2:460306/1‑149 (MQ=255)
                                 aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTAGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAAAGCCACTCTAcc                                                                                                                   >  1:1021694/1‑149 (MQ=255)
                                        gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                                                                                            >  1:1106597/1‑149 (MQ=255)
                                          agaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                                                                                          <  2:336425/149‑1 (MQ=255)
                                              tGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCCGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATATAATAAATTACAGCCACTCTACCCGCAACGACAtca                                                                                                      >  2:1129155/1‑148 (MQ=255)
                                                tCAGGATCGCTGGCGCGTTCTTGCTATTTGCTATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATctct                                                                                                    <  2:28660/149‑1 (MQ=255)
                                                  aGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                                                                                  >  1:115007/1‑149 (MQ=255)
                                                   gTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGaa                                                                                                 <  2:1067419/149‑1 (MQ=255)
                                                        gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATTCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAggtgt                                                                                            >  2:544536/1‑149 (MQ=255)
                                                             cgcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTCCCCGCACCGCCATCTCTGAAGTTGTTGAtt                                                                                       >  2:773634/1‑149 (MQ=255)
                                                              gcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                                                                                      >  2:857683/1‑149 (MQ=255)
                                                              gcgTGCTTGCTATTGGCAATTACCGGGCCGGGGGCAAGTTCTCTCGATCGGCGTTAGGAAAGTCGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                                                                                      <  1:541369/149‑1 (MQ=255)
                                                                    ttGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACTTCGCTGAAGTTGTTGATTCCCCTgg                                                                                >  1:329210/1‑149 (MQ=255)
                                                                    ttGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGGGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTgg                                                                                <  1:113790/149‑1 (MQ=255)
                                                                        tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGca                                                                            >  2:390401/1‑148 (MQ=255)
                                                                         aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa                                                                           <  1:785083/149‑1 (MQ=255)
                                                                         aTTGGCAATTAACGGGCCGGGTGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa                                                                           <  2:9939/149‑1 (MQ=255)
                                                                          ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAATATGACGGATGCCGCGTGAACGCCTTATCCGGCATATTATAAATTACAGACACTCTACACGAAACGCAATCTCTGAAGTTGTTTATTCCCATGGTGCTAg                                                                          >  1:105727/1‑149 (MQ=255)
                                                                          ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCCACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                                                                          >  1:749880/1‑149 (MQ=255)
                                                                          ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                                                                          >  1:807983/1‑149 (MQ=255)
                                                                             gCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGCACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcgc                                                                       >  2:76685/1‑149 (MQ=255)
                                                                               aaTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACTCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAGGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTTGCGCaa                                                                     >  2:525856/1‑149 (MQ=255)
                                                                               aaTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCaa                                                                     >  2:560424/1‑149 (MQ=255)
                                                                                aTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAAt                                                                    <  1:671893/149‑1 (MQ=255)
                                                                                 ttACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATg                                                                   <  2:1079493/149‑1 (MQ=255)
                                                                                  tACCGGGCAGGGGGCAATTGCTCTCGATCGGCGGTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAACTTACAGCCACTCTACCCTCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga                                                                  <  2:912666/149‑1 (MQ=255)
                                                                                     cgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt                                                               >  1:407321/1‑149 (MQ=255)
                                                                                           gggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAcccc                                                         >  2:431727/1‑149 (MQ=255)
                                                                                            ggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAccccc                                                        >  2:321765/1‑149 (MQ=255)
                                                                                               gCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTc                                                     <  1:460306/149‑1 (MQ=255)
                                                                                                 aaTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCGTATCCGGCATTTAATAAATTAAAGCAACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTACCCTGGTGCTAGCGCAATGAGATCAcccaattagg                                                   >  2:547090/1‑142 (MQ=255)
                                                                                                       tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgat                                             <  1:857683/149‑1 (MQ=255)
                                                                                                          cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGTATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgat                                          <  2:891409/149‑1 (MQ=255)
                                                                                                           gATCGGCGTTATGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGAAACGACATCTCTGACGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatg                                         <  1:770350/149‑1 (MQ=255)
                                                                                                           gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatg                                         >  2:1010602/1‑149 (MQ=255)
                                                                                                            aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCGCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatga                                        <  1:797052/149‑1 (MQ=255)
                                                                                                               ggCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatca                                     <  2:105285/149‑1 (MQ=255)
                                                                                                                 cGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCGTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGCCATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatcatc                                   >  2:149845/1‑149 (MQ=255)
                                                                                                                   ttAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCg                                 <  1:653922/149‑1 (MQ=255)
                                                                                                                     aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgc                               <  2:214573/149‑1 (MQ=255)
                                                                                                                      ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcg                              >  2:233469/1‑149 (MQ=255)
                                                                                                                       attaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCGCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                             >  1:1083705/6‑149 (MQ=38)
                                                                                                                       attaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                             >  1:670839/6‑149 (MQ=38)
                                                                                                                         aaGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCAt                           <  1:1172035/149‑1 (MQ=255)
                                                                                                                          aaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATAATCCGGCGCATg                          >  1:1092687/3‑149 (MQ=38)
                                                                                                                           gATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATgg                         <  2:456760/149‑1 (MQ=255)
                                                                                                                           gATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATgg                         <  1:996022/149‑1 (MQ=255)
                                                                                                                               ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCa                     >  2:153167/1‑149 (MQ=255)
                                                                                                                                cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCACCGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTcct                    >  2:723659/1‑147 (MQ=255)
                                                                                                                                cGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCag                    >  1:388600/1‑148 (MQ=255)
                                                                                                                                    tGCGGCGGGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCCACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCGAGCGCAATGAGCTCACCCCCTTCCGGGCGATGAGGATCATCCGGCGCATGGCTCATCGGt                <  2:1153362/149‑1 (MQ=255)
                                                                                                                                     gcggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGtt               >  1:176692/1‑149 (MQ=255)
                                                                                                                                      cggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTc              <  2:352567/149‑1 (MQ=255)
                                                                                                                                       ggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGGTcc             >  1:759329/1‑149 (MQ=255)
                                                                                                                                          gttcaCGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCcaga          <  2:670839/145‑1 (MQ=255)
                                                                                                                                          gttcaCGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCcaga          <  2:1083705/145‑1 (MQ=255)
                                                                                                                                          gTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCcaga          <  1:250791/149‑1 (MQ=255)
                                                                                                                                           tGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACCTCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTGCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCcagac         >  2:445820/1‑149 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCAGGCGTTCcagacag       >  1:127262/1‑149 (MQ=255)
                                                                                                                                               cGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCCAGACagaa     <  2:206240/149‑1 (MQ=255)
                                                                                                                                                  cTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCCAGACAgaagaa  >  1:248013/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTCCAGACAGAAGAA  >  NC_000913/2673650‑2673944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: