Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 680,674 C→T 8.8% P64L (CCG→CTG)  PP_0583 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947680,6740CT 8.8% 186.8 / 6.9 68P64L (CCG→CTG) PP_0583AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base T (2/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ACCTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGGCGCGCG  >  NC_002947/680527‑680823
                                                                                                                                                   |                                                                                                                                                     
acCTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGc                                                                                                                                                     >  2:26844/1‑150 (MQ=255)
 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                                                                                                                                         <  1:457825/150‑1 (MQ=255)
            cccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGgtgt                                                                                                                                         <  2:115337/150‑1 (MQ=255)
                     cccGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCccg                                                                                                                                >  1:506573/1‑150 (MQ=255)
                            cGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGc                                                                                                                         >  2:1234690/1‑150 (MQ=255)
                             gAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGCGGTggggcaa                                                                                                                                           >  2:306127/1‑128 (MQ=255)
                             gAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGCGGTggggcaa                                                                                                                                           <  1:306127/131‑4 (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                                                                                                                            >  2:1020746/1‑130 (MQ=255)
                                             gATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGa                                                                                                                            <  1:1020746/130‑1 (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:945361/1‑150 (MQ=255)
                                                                                                                       acacGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCgg                              >  1:435709/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                                                                              <  1:216206/95‑1 (MQ=255)
                                                                                                                              aCGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGc                                                                              >  2:216206/1‑95 (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:319526/1‑150 (MQ=255)
                                                                                                                                         tttGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGAcgccg            >  1:456075/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)
                                                                                                                                                   |                                                                                                                                                     
ACCTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGCAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGGCGCGCG  >  NC_002947/680527‑680823

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