Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,664,502 T→C 6.0% L122P (CTG→CCG) ‡ acnA‑II → aconitate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,664,5020TC 6.0% 198.0 / 2.3 66L122P (CTG→CCG) ‡acnA‑IIaconitate hydratase
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base C (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

CCCGGCCCGCGTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGC  >  NC_002947/2664356‑2664647
                                                                                                                                                  |                                                                                                                                                 
cccGGCCCGCGTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTg                                                                                                                                                  <  2:296824/148‑1 (MQ=255)
   ggCCCGCGTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCg                                                                                                                                              <  2:4681/149‑1 (MQ=255)
    gCCCGCGTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGt                                                                                                                                             <  6:134741/149‑1 (MQ=255)
      ccGCGTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTc                                                                                                                                            >  2:229000/1‑148 (MQ=255)
       cgcgTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTc                                                                                                                                            >  8:190002/1‑147 (MQ=255)
              tgtgCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCggt                                                                                                                                   >  5:51240/1‑149 (MQ=255)
                tgCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCggtgg                                                                                                                                 >  7:269404/1‑149 (MQ=255)
                tgCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCggtg                                                                                                                                  <  5:28303/148‑1 (MQ=255)
                tgCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGACGGCCAGCGAGTggtcga                                                                                                                                  <  5:195663/148‑6 (MQ=255)
                  ccATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCGCAGGTCAACCCGGAGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGtt                                                                                                                               >  3:261580/1‑149 (MQ=255)
                  ccATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCggtggt                                                                                                                                <  3:232557/148‑1 (MQ=255)
                  ccATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGtt                                                                                                                               <  7:277750/149‑1 (MQ=255)
                     tGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTc                                                                                                                             >  8:276340/1‑148 (MQ=255)
                           ccTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCg                                                                                                                       >  8:54153/1‑148 (MQ=255)
                             ttGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAg                                                                                                                    >  1:165828/1‑149 (MQ=255)
                                      cTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCTACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGaa                                                                                                           <  7:247596/149‑1 (MQ=255)
                                      cTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGaa                                                                                                           <  6:69360/149‑1 (MQ=255)
                                                 gACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACcgcg                                                                                                 <  6:51240/148‑1 (MQ=255)
                                                   ccTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCa                                                                                              >  8:35543/1‑149 (MQ=255)
                                                       gCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATc                                                                                            <  8:269404/147‑1 (MQ=255)
                                                            ccTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCg                                                                                           >  7:177524/1‑143 (MQ=255)
                                                            ccTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATc                                                                                     <  7:54153/149‑1 (MQ=255)
                                                            ccTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATc                                                                                     <  2:113901/149‑1 (MQ=255)
                                                                    aCGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAAc                                                                             >  7:262274/1‑149 (MQ=255)
                                                                      gCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGa                                                                           <  7:276340/149‑1 (MQ=255)
                                                                       ccATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACg                                                                            >  5:170345/1‑147 (MQ=255)
                                                                       ccATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACg                                                                            >  2:169372/1‑147 (MQ=255)
                                                                       ccATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGaa                                                                          <  7:138366/149‑1 (MQ=255)
                                                                        cATTGCCGACAAAGGCGGCGACCCAGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCATTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGaaaa                                                                                                         <  2:228453/117‑1 (MQ=255)
                                                                        cATTGCCGACAAAGGCGGCGACCCAGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCAATCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGaaaa                                                                                                         >  1:228453/1‑117 (MQ=255)
                                                                         aTTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGa                                                                        <  7:190002/149‑1 (MQ=255)
                                                                             ccGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGaa                                                                                                      >  4:191659/1‑115 (MQ=255)
                                                                             ccGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGaa                                                                                                      <  3:191659/115‑1 (MQ=255)
                                                                             ccGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGACGGCCAGCGAGTggtc                                                                                                                                    >  2:239352/1‑82 (MQ=38)
                                                                             ccGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGACGGCCAGCGAGTggtc                                                                                                                                    <  1:239352/85‑4 (MQ=38)
                                                                               gACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCtt                                                                  >  6:203378/1‑149 (MQ=255)
                                                                                         gcgACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTtcct                                                         >  8:79066/1‑146 (MQ=255)
                                                                                           gACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCaa                                                      >  5:202904/1‑149 (MQ=255)
                                                                                            aCCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAAc                                                     >  6:236145/1‑149 (MQ=255)
                                                                                                     aGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGAcc                                               >  6:191870/1‑146 (MQ=255)
                                                                                                       gTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaa                                          <  8:177524/149‑1 (MQ=255)
                                                                                                        tCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaag                                         <  5:203378/149‑1 (MQ=255)
                                                                                                            cccGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGt                                     >  3:244355/1‑149 (MQ=255)
                                                                                                             ccGGTGGTGCCGGTTCAAGTGATCGTCGACCACTGGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCAACAACTGGACCAAGAAGGCGtt                                    <  4:261580/149‑1 (MQ=255)
                                                                                                             ccGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGtt                                    >  1:4065/1‑149 (MQ=255)
                                                                                                               ggtggtGCCGGTGAAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTc                                   >  4:96267/1‑148 (MQ=255)
                                                                                                                 tggtgCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAg                                <  4:244355/149‑1 (MQ=255)
                                                                                                                 tggtgCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACCTCATCAACTGGACCAAGAAGGCGTTCAAg                                <  2:165828/149‑1 (MQ=255)
                                                                                                                      ccGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAAc                             <  8:167871/147‑1 (MQ=255)
                                                                                                                      ccGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAAc                             >  7:167871/1‑147 (MQ=255)
                                                                                                                         gTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTc                                                                 <  7:211541/108‑1 (MQ=255)
                                                                                                                         gTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTc                                                                 >  8:211541/1‑108 (MQ=255)
                                                                                                                             aaCTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGt                     >  7:194782/1‑148 (MQ=255)
                                                                                                                             aaCTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTg                    <  1:169372/149‑1 (MQ=255)
                                                                                                                              aCTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGCCCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATAGCCGCAACGAAGCGCGCTTCCACTTCATCAACTGGACGAAGAAGGCGTTCAAGCACGTCGACGt                     >  8:65084/1‑147 (MQ=255)
                                                                                                                              aCTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAAAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGa                               <  3:277609/137‑1 (MQ=255)
                                                                                                                              aCTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAAAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGa                               >  4:277609/1‑137 (MQ=255)
                                                                                                                                tGATCGTCGACCACTCGCTGGCCGTCGAGTGCTGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaa                                          <  1:283839/124‑1 (MQ=255)
                                                                                                                                tGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaa                                          >  2:283839/1‑124 (MQ=255)
                                                                                                                                  aTCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCa               <  3:96267/149‑1 (MQ=255)
                                                                                                                                      caccgcACTCGACGGCCAGCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCg           >  5:230655/6‑149 (MQ=255)
                                                                                                                                         aCCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGAACCAGCCggg         >  1:259390/1‑148 (MQ=255)
                                                                                                                                            aCTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAAc     <  6:202904/149‑1 (MQ=255)
                                                                                                                                             cTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACg    <  7:35543/149‑1 (MQ=255)
                                                                                                                                              tCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACgg   >  1:174530/1‑149 (MQ=255)
                                                                                                                                               cGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGc  <  8:262274/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCCGGCCCGCGTGGTGTGCCATGACATCCTTGGCCAGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGACCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAGGCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGC  >  NC_002947/2664356‑2664647

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