Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,346,306 C→T 6.5% R65H (CGC→CAC)  PP_2942 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,346,3060CT 6.5% 179.1 / 2.4 61R65H (CGC→CAC) PP_2942two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

TGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGACTCCTCGTCGTC  >  NC_002947/3346166‑3346445
                                                                                                                                            |                                                                                                                                            
tGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCg                                                                                                                                      <  1:83061/149‑1 (MQ=255)
      ggTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCtcat                                                                                                                                 >  2:189234/1‑148 (MQ=255)
        tGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCAt                                                                                                                                        <  4:136777/139‑1 (MQ=255)
        tGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCAt                                                                                                                                        >  3:136777/1‑139 (MQ=255)
        tGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCAc                                                                                                                              >  5:254632/1‑149 (MQ=255)
          gATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCa                                                                                                                            <  6:27806/149‑1 (MQ=255)
           aTCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCAc                                                                                                                              >  2:347105/1‑146 (MQ=255)
               gACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATc                                                                                                                       >  4:37586/1‑148 (MQ=255)
               gACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATc                                                                                                                       >  5:186547/1‑148 (MQ=255)
                 ccGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAAc                                                                                                                    >  2:276549/1‑149 (MQ=255)
                  cGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAAcc                                                                                                                   >  1:57898/1‑149 (MQ=255)
                    tCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCgg                                                                                                                 >  6:202331/1‑149 (MQ=255)
                       ttGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATc                                                                                                                       >  6:302077/1‑140 (MQ=255)
                       ttGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATc                                                                                                                       <  5:302077/140‑1 (MQ=255)
                          aCGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGTGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTc                                                                                                           >  4:232977/1‑149 (MQ=255)
                          aCGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGTGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTc                                                                                                           >  4:237697/1‑149 (MQ=255)
                           cGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCg                                                                                                          <  5:306374/149‑1 (MQ=255)
                                    tgatTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTgc                                                                                                 >  6:370242/1‑149 (MQ=255)
                                          cgcTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGt                                                                                           <  3:113238/149‑1 (MQ=255)
                                          cgcTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGt                                                                                           <  6:254632/149‑1 (MQ=255)
                                                  tCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTc                                                                                   <  5:202331/149‑1 (MQ=255)
                                                  tCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTc                                                                                   <  3:316326/149‑1 (MQ=255)
                                                      cATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTcgccg                                                                               <  6:186547/149‑1 (MQ=255)
                                                             ggTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCg                                                                        >  2:324347/1‑149 (MQ=255)
                                                                 aGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGTGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCa                                                                     <  3:232977/148‑1 (MQ=255)
                                                                 aGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGTGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCa                                                                     <  3:237697/148‑1 (MQ=255)
                                                                   cAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTcg                                                                                  <  1:214392/133‑1 (MQ=255)
                                                                   cAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTcg                                                                                  >  2:214392/1‑133 (MQ=255)
                                                                    aaGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCtt                                                                 >  1:22533/1‑149 (MQ=255)
                                                                     aGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTc                                                                <  3:37586/149‑1 (MQ=255)
                                                                          cTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGt                                                            <  2:22533/148‑1 (MQ=255)
                                                                           tGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTgg                                                                           >  2:321782/1‑132 (MQ=255)
                                                                           tGACGGTCGATGGGTGATGCTGATTGATCTGTGCCACCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTgg                                                                           <  1:321782/132‑1 (MQ=255)
                                                                            gACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCa                                                          >  8:49597/1‑148 (MQ=255)
                                                                             aCGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCa                                                          <  1:189234/147‑1 (MQ=255)
                                                                                   gATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTg                                                  >  4:259216/1‑149 (MQ=255)
                                                                                    aTGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTgg                                                 <  5:370242/149‑1 (MQ=255)
                                                                                     tGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGc                                                <  6:287124/149‑1 (MQ=255)
                                                                                          gATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACg                                            >  4:263499/1‑148 (MQ=255)
                                                                                                gattgatCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGc                                     >  8:4749/1‑149 (MQ=255)
                                                                                                 attgatCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGcc                                    >  7:130242/1‑149 (MQ=255)
                                                                                                    gatCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGAc                                 >  1:77268/1‑149 (MQ=255)
                                                                                                    gatCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGGCATAGGGTTGGCCACGCAATAGCCGa                                  <  1:276549/148‑1 (MQ=255)
                                                                                                     atCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCTCGCATCGCTCATCACCA‑GAGCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGAc                                 <  3:224024/148‑1 (MQ=255)
                                                                                                          tgCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGct                              <  3:259216/146‑2 (MQ=255)
                                                                                                          ggCCCGCAAGGGGGCCCCGCCCAAACCAGGCATGCGCGCATCGGTCATCATCACAATCACTCCCCACATCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCa                              <  1:324347/146‑1 (MQ=255)
                                                                                                            ccAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACt                          <  2:77268/148‑1 (MQ=255)
                                                                                                             cAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑CATCAACCGGGCCCTCCGCCATCAGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTCGCCACGCAATAGCCGACGCAGACTg                         >  3:5284/1‑148 (MQ=255)
                                                                                                              aGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCCGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGc                        >  3:125800/1‑148 (MQ=255)
                                                                                                              aGCAAGGTGGCCCCGTCCATACAAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGc                        >  1:155203/1‑148 (MQ=255)
                                                                                                                cAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTaa                     <  2:57898/149‑1 (MQ=255)
                                                                                                                cAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTa                      >  4:309908/1‑148 (MQ=255)
                                                                                                                 aaGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCCGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAg                    <  4:125800/149‑1 (MQ=255)
                                                                                                                  aGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAgg                   <  7:4749/149‑1 (MQ=255)
                                                                                                                      ggCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCTCCT‑GATCAACCGGGCCCTCGGCCATCTGCGAAAGTCCCTGTTCGCCGCTGGTCGCCAGCTTCACGTGATAGGGTTGTGCACGCAATAGCCGACGCTGCATGCTAAGGCTCg               >  8:167316/1‑149 (MQ=255)
                                                                                                                      ggCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCg               <  7:49597/149‑1 (MQ=255)
                                                                                                                        ccccGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGAc             <  1:347105/149‑1 (MQ=255)
                                                                                                                          ccGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGActc           >  6:50677/1‑149 (MQ=255)
                                                                                                                                aTACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGACTCCtcgtc     >  4:176028/1‑149 (MQ=255)
                                                                                                                                   ccAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGACTCCTcgtcgt   <  3:263499/148‑1 (MQ=255)
                                                                                                                                   ccAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGACTCCTCgtcgtc  >  3:233217/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TGTAACGGTGGATCTGACCGTCATTGACGGCCTTGATGATTGCGCTTGGGTCGGCATAGCCGGTCAGCAAGATCCTGACGGTCGATGGGTGATGCTGATTGATCTGTGCCAGCAAGGTGGCCCCGTCCATACCAGGCATGCGCGCATCGCTCATCACCA‑GATCAACCGGGCCCTCGGCCATCTGCGCCAGTGCCTGTTCGCCGCTGGTCGCCAGCTTCACGTCATAGGGTTGGCCACGCAATAGCCGACGCAGACTGCTAAGGATCGACTCCTCGTCGTC  >  NC_002947/3346166‑3346445

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: ‑