Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,425,141 G→C 10.8% intergenic (‑199/‑601) PP_5647 ← / → PP_3917 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,425,1410GC10.8% 204.6 / 14.1 75intergenic (‑199/‑601)PP_5647/PP_3917hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base C (2/6);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

ATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCT  >  NC_002947/4424994‑4425284
                                                                                                                                                     |                                                                                                                                               
atCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATTTAAGCGGGACAGCGCCAGTAGAGAACCTAGCCCGATCACTGAGCCCGTCCAGCAGACTGAAGTGCAAGTAGCAAAAAGCCCAGcgcgt                                                                                                                                                >  4:20121/1‑148 (MQ=255)
atCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCgg                                                                                                                                                >  1:47662/1‑149 (MQ=255)
 tCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCACAGCGCggg                                                                                                                                               >  2:114553/1‑149 (MQ=255)
         gCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGgggctgggct                                                                                                                                       >  8:118871/1‑149 (MQ=255)
         gCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCCCCgcgctgggc                                                                                                                                        <  7:125400/148‑1 (MQ=255)
           ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTc                                                                                                                                     >  2:295152/1‑149 (MQ=255)
            tAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                    >  6:142079/1‑149 (MQ=255)
              actGCGACCGGTGGAAACTGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCACGTCCAGCAGACAGACGTTCAAGTAGCAAAAAGCCCAGCCCCGCGCTGGGCttttt                                                                                                                                   <  1:358518/146‑4 (MQ=255)
               aTGGGACCGGTTGATAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGATCCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCCCCGCGCTGGGCtttttgc                                                                                                                                 <  6:244937/149‑6 (MQ=255)
               aTGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTc                                                                                                                                 <  1:648/149‑1 (MQ=255)
               aTGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCCCCGCGCTGGGCtttttgc                                                                                                                                 <  6:239473/149‑6 (MQ=255)
                 gCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcag                                                                                                                               >  2:22304/1‑149 (MQ=255)
                  cGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcag                                                                                                                               <  8:217649/148‑1 (MQ=255)
                   gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGTCCCAGCGCGGGGCTGGGCTTCTGTcagc                                                                                                                              <  1:22304/149‑1 (MQ=255)
                   gACCGGTGGAAAATCATCGGCCCTAATCATTGCTGAAAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGACCAGCAGACAGAAGTTAAAGTAGCAAAAATCCCAGCGCGCGGCTGGGCTTATGTcag                                                                                                                               >  2:248438/1‑147 (MQ=255)
                    aCCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                            <  1:248438/149‑1 (MQ=255)
                     ccGGTGGAAAATGATCTGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCaga                                                                                                                           <  4:241711/149‑1 (MQ=255)
                       ggTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCagag                                                                                                                          <  5:130450/148‑1 (MQ=255)
                       ggTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCagag                                                                                                                          >  6:130450/1‑148 (MQ=255)
                             aaaTGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGAc                                                                                                                    >  1:229729/1‑148 (MQ=255)
                             aaaTGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACt                                                                                                                   <  5:194179/149‑1 (MQ=255)
                                      gCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGc                                                                                                          >  5:187554/1‑149 (MQ=255)
                                       cccTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCg                                                                                                         <  6:295703/149‑1 (MQ=255)
                                       cccTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCg                                                                                                         <  7:118871/149‑1 (MQ=255)
                                       cccTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGCTGGGCTTCt                                                                                                                                    <  5:69431/121‑1 (MQ=255)
                                       cccTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGCTGGGCTTCt                                                                                                                                    >  6:69431/1‑121 (MQ=255)
                                         cTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGtt                                                                                                       >  4:118843/1‑149 (MQ=255)
                                               cTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAAtt                                                                                                  <  6:18278/148‑1 (MQ=255)
                                                ccgCTGATATCTATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTAc                                                                                                >  3:36726/3‑149 (MQ=255)
                                                  gCTGATAGCTAGCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGAGCACTGAGCCCGTCCAGCAGACAGAAGTGCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCt                                                                                                                        <  5:72451/123‑1 (MQ=255)
                                                  gCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCCGCAGACAGAAGTGCAAGGAGCAAAAAGCCCCGCGCGGGGCGGGGCTTCTGTCAGCAGAGcg                                                                                                                        >  6:72451/1‑122 (MQ=255)
                                                  gCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGc                                                                                              >  4:18766/1‑149 (MQ=255)
                                                      atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                          >  4:345546/1‑149 (MQ=255)
                                                               aGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTc                                                                                 <  1:114553/149‑1 (MQ=255)
                                                               aGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTc                                                                                 <  1:295152/149‑1 (MQ=255)
                                                                     aCAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATgg                                                                           <  1:263203/149‑1 (MQ=255)
                                                                       aGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTg                                                                                           >  2:252662/1‑131 (MQ=255)
                                                                       aGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTg                                                                                           <  1:252662/131‑1 (MQ=255)
                                                                         cgcCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGgcgc                                                                       <  3:18766/149‑1 (MQ=255)
                                                                             aGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTtgt                                                                   <  3:118843/149‑1 (MQ=255)
                                                                             aGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGAGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTtgt                                                                   >  7:317603/1‑149 (MQ=255)
                                                                                 gagaACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGa                                                                 >  8:251791/1‑147 (MQ=255)
                                                                                  agaACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAg                                                              >  6:13162/1‑149 (MQ=255)
                                                                                  agaACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCAAATGGGCGCTTGTGATAg                                                              >  2:242213/1‑149 (MQ=255)
                                                                                    aaCCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCt                                                            <  2:229729/149‑1 (MQ=255)
                                                                                      ccGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCta                                                           >  3:258935/1‑148 (MQ=255)
                                                                                       cGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTATCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCt                                                            >  1:68556/1‑147 (MQ=255)
                                                                                       cGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCt                                                            <  2:68556/146‑1 (MQ=255)
                                                                                        gAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAg                                                        <  7:251791/149‑1 (MQ=255)
                                                                                                    tGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTg                                            >  7:114014/1‑149 (MQ=255)
                                                                                                           gTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTCCCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGcac                                     >  8:306345/1‑147 (MQ=255)
                                                                                                            tCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCt                                    <  4:258935/149‑1 (MQ=255)
                                                                                                              cagcagACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTc                                   >  4:263629/1‑148 (MQ=255)
                                                                                                                  agacagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATgg                              <  6:187554/149‑1 (MQ=255)
                                                                                                                  agacagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATgg                              >  7:215701/1‑149 (MQ=255)
                                                                                                                        aaGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCgg                        <  5:142079/149‑1 (MQ=255)
                                                                                                                          gTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTc                                                                                 <  6:313570/90‑1 (MQ=255)
                                                                                                                          gTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTc                                                                                 >  5:313570/1‑90 (MQ=255)
                                                                                                                          gTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCg                         <  4:144448/146‑1 (MQ=255)
                                                                                                                          gTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCg                         >  3:144448/1‑146 (MQ=255)
                                                                                                                           ttCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTg                                            >  4:158355/1‑126 (MQ=255)
                                                                                                                           ttCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTg                                            <  3:158355/126‑1 (MQ=255)
                                                                                                                           ttCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTgg                     <  1:242213/149‑1 (MQ=255)
                                                                                                                                    gCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGt              <  8:114014/148‑1 (MQ=255)
                                                                                                                                      aaaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTc           >  6:362468/1‑149 (MQ=255)
                                                                                                                                      aaaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTc           >  8:233378/1‑149 (MQ=255)
                                                                                                                                      aaaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTc           >  6:226667/1‑149 (MQ=255)
                                                                                                                                       agaaGCCCAGCCCCGCGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGc                                                >  6:315969/3‑111 (MQ=255)
                                                                                                                                       agaaGCCCAGCCCCGCGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGc                                                <  5:315969/109‑1 (MQ=255)
                                                                                                                                       aaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCt          >  8:32451/1‑149 (MQ=255)
                                                                                                                                         aaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgc        <  8:215701/149‑1 (MQ=255)
                                                                                                                                         aaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgc        <  5:226667/149‑1 (MQ=255)
                                                                                                                                         aaGCCCAGCCCCGCGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCtt                >  2:180265/1‑141 (MQ=255)
                                                                                                                                         aaGCCCAGCCCCGCGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCtt                <  1:180265/141‑1 (MQ=255)
                                                                                                                                              ccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                         <  7:292688/112‑1 (MQ=255)
                                                                                                                                              ccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                         >  8:292688/1‑112 (MQ=255)
                                                                                                                                               cAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTTTTTCTGCGCGcc   >  4:70810/1‑149 (MQ=255)
                                                                                                                                                aGCGCGCGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGATATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCt  >  6:92526/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCT  >  NC_002947/4424994‑4425284

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