Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,160,645 G→C 7.7% R168G (CGT→GGT)  ygiQ ← Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,160,6450GC 7.7% 409.2 / 20.8 117R168G (CGT→GGT) ygiQRadical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (48/60);  new base C (3/6);  total (51/66)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 117 total aligned reads displayed.

ACTAATTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCATGACGTAAAC  >  NC_000913/3160502‑3160788
                                                                                                                                               |                                                                                                                                               
 cTAATTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAg                                                                                                                                           >  3:801660/1‑149 (MQ=255)
  tAATTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGa                                                                                                                                          <  3:309784/149‑1 (MQ=255)
  tAATTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGa                                                                                                                                          >  1:521322/1‑149 (MQ=255)
   aaTTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCatta                                                                                                                                         <  1:366718/149‑4 (MQ=255)
    aTTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACt                                                                                                                                        >  2:252127/1‑149 (MQ=255)
    aTTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCattat                                                                                                                                        <  1:851930/149‑5 (MQ=255)
     ttGGCTCGCCCATCGCCAGACGATGCGCCACTTGCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTa                                                                                                                                       <  4:950342/149‑1 (MQ=255)
     ttGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTa                                                                                                                                       >  3:1076989/1‑149 (MQ=255)
     ttGACTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCattatg                                                                                                                                       <  3:1051428/149‑6 (MQ=255)
        gCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGcc                                                                                                                                    >  1:45466/1‑149 (MQ=255)
          tCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTc                                                                                                                                  >  4:966995/1‑149 (MQ=255)
            gCCCATCGCCAGACGATGCGCCCCTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCaa                                                                                                                                <  1:593362/149‑1 (MQ=255)
            gCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCaa                                                                                                                                >  1:445673/1‑149 (MQ=255)
              ccATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATa                                                                                                                              >  2:161283/1‑149 (MQ=255)
                aTCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATAcc                                                                                                                            >  3:718441/1‑149 (MQ=255)
                aTCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCattatg                                                                                                                                       >  1:1087228/1‑133 (MQ=255)
                aTCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCattatg                                                                                                                                       <  2:1087228/138‑6 (MQ=255)
                 tCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATAccg                                                                                                                           >  3:934164/1‑149 (MQ=255)
                      aGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAg                                                                                                                      <  1:252127/149‑1 (MQ=255)
                      aGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAg                                                                                                                      <  1:839598/149‑1 (MQ=255)
                       gACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTc                                                                                                                                  >  2:164848/1‑136 (MQ=255)
                        aCGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGa                                                                                                                    >  1:83963/1‑149 (MQ=255)
                         cGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGAt                                                                                                                   <  4:391140/149‑1 (MQ=255)
                         cGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGAt                                                                                                                   <  3:311883/149‑1 (MQ=255)
                         cGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGAt                                                                                                                   <  1:214343/149‑1 (MQ=255)
                          gATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATc                                                                                                                  >  2:835606/1‑149 (MQ=255)
                            tGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCAc                                                                                                                <  1:1087407/149‑1 (MQ=255)
                            tGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCAc                                                                                                                <  1:189084/149‑1 (MQ=255)
                             gcgcCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCAcc                                                                                                               <  2:150020/149‑1 (MQ=255)
                             gcgcCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCAcc                                                                                                               <  4:718441/149‑1 (MQ=255)
                                ccACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGt                                                                                                            >  1:769825/1‑149 (MQ=255)
                                ccACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGt                                                                                                            >  2:1035310/1‑149 (MQ=255)
                                        accaGCGTACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGa                                                                                                                                          >  2:597184/1‑111 (MQ=255)
                                        accaGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGa                                                                                                                                          <  1:597184/111‑1 (MQ=255)
                                        accaGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCttt                                                                                                    <  2:492941/149‑1 (MQ=255)
                                        accaGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCat                                                                                                                                           >  3:686791/1‑109 (MQ=255)
                                        accaGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCCGTACCGCGCat                                                                                                                                           <  4:686791/110‑2 (MQ=255)
                                          caGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTcc                                                                                                  >  4:1018696/1‑149 (MQ=255)
                                               gACGCTCACCGATACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGAGCACCGGTACATCTTTCCACGcc                                                                                             >  3:612728/1‑149 (MQ=255)
                                                aCGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCct                                                                                            <  2:974421/149‑1 (MQ=255)
                                                  gCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCctct                                                                                          >  4:571516/1‑149 (MQ=255)
                                                    tCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCttt                                                                                        >  1:1024911/1‑149 (MQ=255)
                                                      aCCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTAc                                                                                      >  2:578093/1‑149 (MQ=255)
                                                        cGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACaa                                                                                    <  3:216379/149‑1 (MQ=255)
                                                         gTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAAc                                                                                   <  3:571516/149‑1 (MQ=255)
                                                            aCCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCt                                                                                <  3:362164/149‑1 (MQ=255)
                                                            aCCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCt                                                                                <  4:801660/149‑1 (MQ=255)
                                                            aCCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCt                                                                                <  4:80298/149‑1 (MQ=255)
                                                              cAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAg                                                                                                                      >  1:459586/1‑109 (MQ=255)
                                                              cAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAg                                                                                                                      <  2:459586/109‑1 (MQ=255)
                                                              cAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTgg                                                                              <  4:35690/149‑1 (MQ=255)
                                                              cAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTgg                                                                              <  1:726643/149‑1 (MQ=255)
                                                               aaaCATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACaa                                                                                    <  2:26596/142‑1 (MQ=255)
                                                               aaaCATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACaa                                                                                    >  1:26596/1‑142 (MQ=255)
                                                                 aCATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGta                                                                           >  4:821993/1‑149 (MQ=255)
                                                                   aTCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGtata                                                                         <  3:266512/149‑1 (MQ=255)
                                                                        cATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCa                                                                    <  3:103240/149‑1 (MQ=255)
                                                                        cATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCa                                                                    <  4:923709/149‑1 (MQ=255)
                                                                          tGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAgt                                                                  <  1:835606/149‑1 (MQ=255)
                                                                           gTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAgtg                                                                 >  2:35207/1‑149 (MQ=255)
                                                                             cGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGtgtg                                                               <  2:355744/149‑1 (MQ=255)
                                                                              ggCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTgg                                                              >  4:61892/1‑149 (MQ=255)
                                                                                cTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGgcg                                                            >  1:421385/1‑149 (MQ=255)
                                                                                cTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGgcg                                                            >  4:743030/1‑149 (MQ=255)
                                                                                 tttCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGgcgc                                                           >  4:887094/1‑149 (MQ=255)
                                                                                      aaTCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATcc                                                      <  2:953806/149‑1 (MQ=255)
                                                                                          caccaGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGgc                                                  <  1:35207/149‑1 (MQ=255)
                                                                                           accaGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTACACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGgcg                                                 >  2:909979/1‑149 (MQ=255)
                                                                                              agctcagAACGGCGCACGGTATCGGACAAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCtt                                              <  4:934164/143‑1 (MQ=255)
                                                                                              aGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTAGAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCtt                                              <  4:626260/149‑1 (MQ=255)
                                                                                                 aCGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTAcc                                           >  3:759837/1‑149 (MQ=255)
                                                                                                  cGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCCCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTAccc                                          >  4:75753/1‑149 (MQ=255)
                                                                                                   ggAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCcg                                         >  2:6457/1‑149 (MQ=255)
                                                                                                    gAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCcgc                                        >  2:468204/1‑149 (MQ=255)
                                                                                                      aCGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGa                                      <  1:750461/149‑1 (MQ=255)
                                                                                                      aCGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGa                                      <  4:1076989/149‑1 (MQ=255)
                                                                                                          cgcACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGtt                                  >  1:697147/1‑149 (MQ=255)
                                                                                                           gcACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTa                                 <  1:1032330/149‑1 (MQ=255)
                                                                                                           gcACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTa                                 <  3:743030/149‑1 (MQ=255)
                                                                                                                gTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGtt                                  <  2:734495/143‑1 (MQ=255)
                                                                                                                gTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGtt                                  >  1:734495/1‑143 (MQ=255)
                                                                                                                  aTCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCg                          >  4:953166/1‑149 (MQ=255)
                                                                                                                   tCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCgg                                                    <  1:1019530/122‑1 (MQ=255)
                                                                                                                   tCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCgg                                                    >  2:1019530/1‑122 (MQ=255)
                                                                                                                     ggACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCgtgt                       <  3:358121/149‑1 (MQ=255)
                                                                                                                     ggACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCgtgt                       >  1:905693/1‑149 (MQ=255)
                                                                                                                       aCCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCGCCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAATCAGTGTGTCGCGATCCGGGCGCTTACCCGCGGCGTTATCCGGCGTGTAg                     <  4:667871/149‑1 (MQ=255)
                                                                                                                       aCCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAg                     <  4:511536/149‑1 (MQ=255)
                                                                                                                       aCCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAg                     <  2:486350/149‑1 (MQ=255)
                                                                                                                         cAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGc                   >  3:574420/1‑149 (MQ=255)
                                                                                                                           gTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCAt                 >  4:474271/1‑149 (MQ=255)
                                                                                                                            tAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATc                <  4:248624/149‑1 (MQ=255)
                                                                                                                               tCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTc             <  3:295842/149‑1 (MQ=255)
                                                                                                                                cATAATGCGCGGTACGGCGCAGACTAGCCTGAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCa            <  3:645924/149‑1 (MQ=255)
                                                                                                                                 aTAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGTGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCAt           <  2:298121/149‑1 (MQ=255)
                                                                                                                                 aTAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCAt           <  3:821993/149‑1 (MQ=255)
                                                                                                                                 aTAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCAt           <  2:769825/149‑1 (MQ=255)
                                                                                                                                   tctGCGCCGTACCGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCATGa         >  1:476494/3‑149 (MQ=255)
                                                                                                                                    aTGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCATGAc        >  1:1091478/1‑149 (MQ=255)
                                                                                                                                    aTGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCATGAc        <  2:355125/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ACTAATTGGCTCGCCCATCGCCAGACGATGCGCCACTTCCACCAGCGGACGCTCACCGTTACCAAACATCAGCATGTCGGCTTTCGAATCCACCAGCACGGAACGGCGCACGGTATCGGACCAGTAATCATAATGCGCGGTACGGCGCAGACTAGCCTCAATACCGCCGAGGATCACCGGTACATCTTTCCACGCCTCTTTACAACGCTGGGTATAAACCAGTGTGGCGCGATCCGGGCGCTTACCCGCGACGTTATCCGGCGTGTAGGCATCGTCATGACGTAAAC  >  NC_000913/3160502‑3160788

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: