Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,582,559 A→G 6.6% intergenic (+136/+31) PP_4910 → / ← PP_4911 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,582,5590AG 6.6% 225.7 / 4.2 76intergenic (+136/+31)PP_4910/PP_4911hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base G (3/2);  total (40/36)
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.35e-01

CGCGCCTGACAGGTAGACGCGACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGG  >  NC_002947/5582415‑5582705
                                                                                                                                                 |                                                                                                                                                  
cgcgCCTGACAGGTAGACGCGACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGc                                                                                                                                                <  7:122828/149‑1 (MQ=255)
  cgcCTGACAGGTAGACGCGACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCgg                                                                                                                                              <  2:241332/149‑1 (MQ=255)
        aCAGGTAGACGCGACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCATTCGCTGCGGGCtttt                                                                                                                                        >  5:90143/1‑148 (MQ=255)
            gTAGACGCGACATCTATAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCgg                                                                                                                                              <  8:1476/139‑1 (MQ=255)
            gTAGACGCGACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGt                                                                                                                                    <  3:19195/149‑1 (MQ=255)
            gTAGACGCGACATCTAAAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCgg                                                                                                                                              >  7:1476/1‑139 (MQ=255)
                  gcgACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTgg                                                                                                                              >  5:192889/1‑149 (MQ=255)
                  gcgACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTgg                                                                                                                              >  4:85113/1‑149 (MQ=255)
                    gACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCATTCGCTGCGGGCTTTTTTGTATCTgccg                                                                                                                            <  6:90143/149‑4 (MQ=38)
                     aCATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCg                                                                                                                           <  4:82978/149‑1 (MQ=255)
                            cAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGt                                                                                                                          >  7:152778/1‑143 (MQ=255)
                            cAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGt                                                                                                                          <  8:152778/143‑1 (MQ=255)
                            cAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGa                                                                                                                      <  2:73852/147‑1 (MQ=255)
                            cAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGAt                                                                                                                     <  1:63422/148‑1 (MQ=255)
                            cAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATc                                                                                                                    <  6:141981/149‑1 (MQ=255)
                                  aTGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCa                                                                                                                   <  2:156861/144‑1 (MQ=255)
                                  aTGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCa                                                                                                                   >  1:156861/1‑144 (MQ=255)
                                       ttGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAAGTAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGAt                                                                                                           >  1:83101/1‑148 (MQ=255)
                                         gAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAAc                                                                                                       >  6:155446/1‑149 (MQ=255)
                                            ccAGGCTGGTGTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTc                                                                                                    <  6:192889/149‑1 (MQ=255)
                                             cAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTcc                                                                                                   <  3:11191/149‑1 (MQ=255)
                                             cAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTc                                                                                                    >  6:17874/1‑148 (MQ=255)
                                              aGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCa                                                                                                  >  1:159489/1‑149 (MQ=255)
                                                    gTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTAt                                                                                                                                  >  6:199596/1‑111 (MQ=255)
                                                    gTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTAt                                                                                                                                  <  5:199596/111‑1 (MQ=255)
                                                        ccTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGc                                                                                         >  5:236655/1‑148 (MQ=255)
                                                          ttGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGaa                                                                                      >  5:320824/1‑149 (MQ=255)
                                                          ttGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGAGaa                                                                                      >  8:213367/1‑149 (MQ=255)
                                                             cGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCg                                                                                   <  5:17874/149‑1 (MQ=255)
                                                              gATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGGGGGCttt                                                                                                                                         <  5:229671/94‑1 (MQ=255)
                                                              gATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGAGGGCttt                                                                                                                                         >  6:229671/1‑94 (MQ=255)
                                                                  aCCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAg                                                                              >  5:315295/1‑149 (MQ=255)
                                                                   cctgcctgGAGAAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCACGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAgg                                                                             >  6:116040/1‑149 (MQ=255)
                                                                    ctgcctgGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGc                                                                            >  6:309011/1‑149 (MQ=255)
                                                                       cctgGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCAt                                                                         >  2:58768/1‑149 (MQ=255)
                                                                                 gATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTg                                                               <  6:315295/149‑1 (MQ=255)
                                                                                 gATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTg                                                               <  7:197293/149‑1 (MQ=255)
                                                                                 gATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCAGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTg                                                               >  8:251662/1‑149 (MQ=255)
                                                                                  aTCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTgg                                                              <  1:58768/149‑1 (MQ=255)
                                                                                   tCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTgg                                                              <  2:83101/148‑1 (MQ=255)
                                                                                     cAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGtt                                                           >  8:99782/1‑149 (MQ=255)
                                                                                      aGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGACATCGGGATTTTGGGttt                                                          >  5:96698/1‑149 (MQ=255)
                                                                                          aaTAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTcc                                                      >  5:146191/1‑149 (MQ=255)
                                                                                                cgcAAACCCGCTCCAACGAAGACCGGCAGATACAAAAAAGCCAGCACCGAAAGCGGGCTTTGTCGTATCTGCCCGATTGATCAGAACCCATAAACTACCATATCGGAGCGAATCGTCGACGCCATCGCGATTATGGGTTTCTCCGcccc                                                 >  8:226981/1‑149 (MQ=255)
                                                                                                 gcAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGAtttt                                                                <  5:182162/133‑1 (MQ=255)
                                                                                                 gcAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGAtttt                                                                >  6:182162/1‑133 (MQ=255)
                                                                                                   aaaCCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGGATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCctt                                              >  6:70240/1‑149 (MQ=255)
                                                                                                    aaCCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCcttc                                             <  5:116040/149‑1 (MQ=255)
                                                                                                      cccGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCcttctt                                           <  7:251662/149‑1 (MQ=255)
                                                                                                      cccGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCcttctt                                           >  1:135525/1‑149 (MQ=255)
                                                                                                      cccGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCcttct                                            >  4:294748/1‑148 (MQ=255)
                                                                                                        cGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGGATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCa                                         <  5:70240/149‑1 (MQ=255)
                                                                                                              tACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAg                                   >  6:264681/1‑149 (MQ=255)
                                                                                                                cGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCgg                                                    >  8:44754/1‑130 (MQ=255)
                                                                                                                cGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCgg                                                    <  7:44754/130‑1 (MQ=255)
                                                                                                                 gAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCg                                <  2:159489/149‑1 (MQ=255)
                                                                                                                 gAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCg                                >  2:213357/1‑149 (MQ=255)
                                                                                                                    gACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTg                             <  7:99782/149‑1 (MQ=255)
                                                                                                                      ccGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCa                           <  6:236655/149‑1 (MQ=255)
                                                                                                                       ggccAGATACGACAAAGCCCGCATTCGCTGCGGGCTTTTTTGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAg                          >  4:270008/4‑149 (MQ=37)
                                                                                                                         gCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCg                        >  3:123083/1‑149 (MQ=255)
                                                                                                                            gATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGc                     <  5:309011/149‑1 (MQ=255)
                                                                                                                             aTACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCt                    >  5:292794/1‑149 (MQ=255)
                                                                                                                               aCAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCTTCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCttt                  >  5:251663/1‑149 (MQ=255)
                                                                                                                                 aaaaaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCa                           >  1:118084/1‑138 (MQ=255)
                                                                                                                                 aaaaaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCa                           <  2:118084/138‑1 (MQ=255)
                                                                                                                                 aaaaaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCttttt                >  2:270081/1‑149 (MQ=255)
                                                                                                                                 aaaaaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCttttt                >  3:224988/1‑149 (MQ=255)
                                                                                                                                   aaaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGTGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCg              <  7:213367/149‑1 (MQ=255)
                                                                                                                                   aaaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCg              <  6:146191/149‑1 (MQ=255)
                                                                                                                                    aaaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCgg             <  5:264681/149‑1 (MQ=255)
                                                                                                                                     aaGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCGGt            <  5:155446/149‑1 (MQ=255)
                                                                                                                                       gCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCGGTGc          <  3:85113/149‑1 (MQ=255)
                                                                                                                                           ggtgCGAATGCGGGCTTTGTCGTATGTGGCCTTTTGATCTGAAGGGATATACTTCCATTTCGGTGCGAATAGGAGAGGCCATGTGGATTTTGGGTTTCTGAGGAGCCTTCTTATCCTGAGTCGGTGAAGGGGGCTATTTAGGTGcgtat      >  3:162172/4‑145 (MQ=255)
                                                                                                                                           gCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCAATATCGGTGCTAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCAGGCCCCTTCTTCAACTGAGTCGGTGCATCGGGCTTTTTCGGTTCctct      >  3:162169/1‑149 (MQ=255)
                                                                                                                                               cGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCgg  <  2:135525/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCGCCTGACAGGTAGACGCGACATCTACAGCGCATGTCTTGAGCCAGGCTGGTTTCCTTGCGATCACCTGCCTGGCGCAAGATCCAGCCAATAGCGCAAACCCGCTCCTACGAAGACCGGCAGATACAAAAAAGCCCGCAGCGAATGCGGGCTTTGTCGTATCTGGCCGTTTGATCAGAACCGATAAACTTCCATATCGGTGCGAATCGGCGAGGCCATCGGGATTTTGGGTTTCTCCGGCCCCTTCTTCACCTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGG  >  NC_002947/5582415‑5582705

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