Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 680,677 A→G 9.0% Q65R (CAG→CGG)  PP_0583 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947680,6770AG 9.0% 180.3 / 7.6 67Q65R (CAG→CGG) PP_0583AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (2/4);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

CCTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGGCGCGCG  >  NC_002947/680528‑680823
                                                                                                                                                     |                                                                                                                                                  
ccTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCa                                                                                                                                                    <  2:1045495/150‑1 (MQ=255)
 cTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAg                                                                                                                                                   >  2:1257340/1‑150 (MQ=255)
  tGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGt                                                                                                                                                  >  1:456369/1‑150 (MQ=255)
     ccAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTggg                                                                                                                                               <  2:456369/150‑1 (MQ=255)
     ccAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTggg                                                                                                                                               <  1:992635/150‑1 (MQ=255)
       aGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTggggg                                                                                                                                             >  1:1201221/1‑150 (MQ=255)
       aGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTggggg                                                                                                                                             >  2:68759/1‑150 (MQ=255)
        gACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGCGGTggggca                                                                                                                                            <  2:975782/150‑3 (MQ=255)
          ccccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCACCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGCGGTggggcaaa                                                                                                                                          <  1:574661/150‑5 (MQ=255)
           cccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGgtgt                                                                                                                                         >  2:890830/1‑150 (MQ=255)
           cccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGgtgt                                                                                                                                         <  2:115337/150‑1 (MQ=255)
           cccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGgtgt                                                                                                                                         <  1:457825/150‑1 (MQ=255)
                    cccGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCccg                                                                                                                                >  1:506573/1‑150 (MQ=255)
                           cGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGc                                                                                                                         >  2:1234690/1‑150 (MQ=255)
                            gAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGCGGTggggcaa                                                                                                                                           <  1:306127/131‑4 (MQ=255)
                            gAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGCGGTggggcaa                                                                                                                                           >  2:306127/1‑128 (MQ=255)
                              aaGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGa                                                                                                                      <  2:1201221/150‑1 (MQ=255)
                                  cTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCCCCGCCGTGGGGGGTGTGGTCCCCGGCGGATGCCGAGGAt                                                                                                                  >  2:586577/1‑150 (MQ=255)
                                       cgcCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGt                                                                                                             <  1:68759/150‑1 (MQ=255)
                                          cGGATCGTGGACCCCTCGGCACCAGCACGACTGGGTGCAATTCCCCTACGCCATCAGCGGCGTACTCGGCGGGTACACACGCGACGGCAGTTACGTGGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGt                                                                                                          <  2:635718/150‑1 (MQ=255)
                                            gATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGa                                                                                                                            <  1:1020746/130‑1 (MQ=255)
                                            gATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGa                                                                                                                            >  2:1020746/1‑130 (MQ=255)
                                                gTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTc                                                                                                    <  1:955348/150‑1 (MQ=255)
                                                  ggACCCCTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGa                                                                                                  <  1:1252898/150‑1 (MQ=255)
                                                   gACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGAt                                                                                                 <  2:1159220/150‑1 (MQ=255)
                                                        cTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAg                                                                                            <  2:317851/150‑1 (MQ=255)
                                                            ccaccaGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCga                                                                                        <  1:890830/150‑1 (MQ=255)
                                                              accaGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCgaga                                                                                      >  1:128938/1‑150 (MQ=255)
                                                                    aCGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCa                                                                                <  1:783954/150‑1 (MQ=255)
                                                                    aCGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCa                                                                                <  1:1257340/150‑1 (MQ=255)
                                                                      gACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGAGTGAAACCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGc                                                                              >  1:655836/1‑150 (MQ=255)
                                                                      caCTGGGTACAATCCTCCCACGCCATCACCGCCGTACTCGGCTTGTACACCCGCCACGCCAGTTACTCTGCCCCACCGCAGCGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGc                                                                              <  2:1092243/149‑1 (MQ=255)
                                                                           ggTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTa                                                                         >  1:1236407/1‑150 (MQ=255)
                                                                            gTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTAc                                                                        <  1:1234690/150‑1 (MQ=255)
                                                                                aaTTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTgc                                                                    >  1:50437/1‑150 (MQ=255)
                                                                                 aTTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTgcg                                                                   >  1:365935/1‑150 (MQ=255)
                                                                                        tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGc                                                                                                                    >  1:700885/1‑94 (MQ=255)
                                                                                        tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGc                                                                                                                    <  2:700885/94‑1 (MQ=255)
                                                                                        tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTgcgcgccgag                                                            >  2:232093/1‑144 (MQ=255)
                                                                                        tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGAt                                                            >  1:913565/1‑150 (MQ=255)
                                                                                        tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGGTGATCCCCGCCGCTGCCGAGCATGAAGTTGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGAt                                                            >  1:607323/1‑150 (MQ=255)
                                                                                                     gcgTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATggg                                               >  2:701240/1‑150 (MQ=255)
                                                                                                           tCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCgg                                         >  1:547397/1‑150 (MQ=255)
                                                                                                                gtgtACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACaa                                    <  2:60195/150‑1 (MQ=255)
                                                                                                                  gtACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATg                                  <  2:506573/150‑1 (MQ=255)
                                                                                                                    acacacGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGcc                                <  2:547397/150‑1 (MQ=255)
                                                                                                                      acacGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCgg                              >  1:435709/1‑150 (MQ=255)
                                                                                                                      acacGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCgg                              >  1:945361/1‑150 (MQ=255)
                                                                                                                      acacGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGCAACAATGCCgg                              <  2:50437/150‑1 (MQ=255)
                                                                                                                        acGCGACGTCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGt                            <  2:1218882/150‑1 (MQ=255)
                                                                                                                           cgACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCt                         <  1:652341/150‑1 (MQ=255)
                                                                                                                             aCGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGc                                                                              >  2:216206/1‑95 (MQ=255)
                                                                                                                             aCGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGc                                                                              <  1:216206/95‑1 (MQ=255)
                                                                                                                               ggCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCTGAACAATGCCGGGTGCTGGaa                     <  1:26844/150‑1 (MQ=255)
                                                                                                                               ggCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGaa                     >  1:388740/1‑150 (MQ=255)
                                                                                                                                 cAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGt                   <  1:586577/150‑1 (MQ=255)
                                                                                                                                 cAGTTACTTTGCCCCACCGCAGTGGGGGATGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGt                   >  2:722315/1‑150 (MQ=255)
                                                                                                                                  aGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTg                  >  2:832295/1‑150 (MQ=255)
                                                                                                                                        tttGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGAcgccg            >  1:456075/1‑150 (MQ=255)
                                                                                                                                        tttGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGAcgccg            >  1:319526/1‑150 (MQ=255)
                                                                                                                                           cccccACTGCGGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGa                                                                                                  <  1:943182/60‑1 (MQ=255)
                                                                                                                                           cccccACTGCGGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGa                                                                                                  >  2:943182/2‑61 (MQ=255)
                                                                                                                                              ccACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGgcg      >  1:490435/1‑150 (MQ=255)
                                                                                                                                               cACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCg                               >  2:56875/1‑124 (MQ=255)
                                                                                                                                               cACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCg                               <  1:56875/124‑1 (MQ=255)
                                                                                                                                               cACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGgcgc     <  1:832295/150‑1 (MQ=255)
                                                                                                                                               cACCGCAGTGGGGGATGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGgcgc     <  1:722315/150‑1 (MQ=255)
                                                                                                                                                  cGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGGcgcgcg  >  1:332815/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGGCGCGCG  >  NC_002947/680528‑680823

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