Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,162,527 G→A 8.0% intergenic (+474/‑106) PP_2776 → / → PP_2777 homocysteine S‑methyltransferase family protein/acyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,162,5270GA 8.0% 172.4 / 3.8 63intergenic (+474/‑106)PP_2776/PP_2777homocysteine S‑methyltransferase family protein/acyl carrier protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/3);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01

CTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAG  >  NC_002947/3162388‑3162670
                                                                                                                                             |                                                                                                                                               
cTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCCGGATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGgag                                                                                                                                           >  1:98736/1‑148 (MQ=255)
cTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCCTAAGGGgagagt                                                                                                                                        <  1:163480/149‑3 (MQ=255)
cTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCCTAAGGGgagagt                                                                                                                                        <  5:227992/149‑3 (MQ=255)
 tATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCGAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACta                                                                                                                                       >  8:118743/1‑148 (MQ=255)
 tATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACtt                                                                                                                                       <  7:290136/149‑1 (MQ=255)
   tCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACttt                                                                                                                                      <  7:66417/148‑1 (MQ=255)
   tCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACttt                                                                                                                                      <  8:104850/148‑1 (MQ=255)
           tttCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGTCTAAAGTCTcccc                                                                                                                             <  4:34752/149‑2 (MQ=255)
           tttCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCg                                                                                                                             >  5:208820/1‑149 (MQ=255)
             tCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCgg                                                                                                                            <  3:124378/148‑1 (MQ=255)
              cAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGc                                                                                                                           <  8:329331/148‑1 (MQ=255)
                   aaTCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCg                                                                                                                      >  2:13263/1‑148 (MQ=255)
                    aTCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTa                                                                                                                                     <  4:38284/132‑1 (MQ=255)
                    aTCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTa                                                                                                                                     >  3:38284/1‑132 (MQ=255)
                    aTCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTa                                                                                                                                     <  2:33770/132‑1 (MQ=255)
                    aTCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTa                                                                                                                                     >  1:33770/1‑132 (MQ=255)
                        aTTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGa                                                                                                                 >  1:277234/1‑148 (MQ=255)
                             gTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTAATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGc                                                                                                           >  2:165063/1‑149 (MQ=255)
                                 ttttACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACtt                                                                                                                                       <  6:354760/117‑1 (MQ=255)
                                 ttttACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACtt                                                                                                                                       >  5:354760/1‑117 (MQ=255)
                                      cGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACa                                                                                                  <  6:208820/149‑1 (MQ=255)
                                          cccTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCtt                                                                                              >  4:36250/1‑149 (MQ=255)
                                                 cAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCt                                                                                                              <  1:249304/126‑1 (MQ=255)
                                                 cAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCt                                                                                                              >  2:249304/1‑126 (MQ=255)
                                                                     aTGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTattcat                                                                     >  5:247291/1‑149 (MQ=255)
                                                                      tGAATCAAAGGATTCAGAGCGCCCACGCATTTCTCACAAAAAGCACTGATGGCACCTTGACTCACCCATAAGGGGGAGACTTTAGAATCCCGGCACCCCATAGATCTGGCCTCAAAACACCTATATTTTGCACGATGACATTattcat                                                                     >  2:311614/1‑148 (MQ=255)
                                                                        aaTCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTattcattc                                                                   <  5:335285/148‑1 (MQ=255)
                                                                          tCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAAt                                                                >  4:333893/1‑149 (MQ=255)
                                                                          tCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAAt                                                                >  4:117783/1‑149 (MQ=255)
                                                                           cAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATg                                                               >  3:243467/1‑149 (MQ=255)
                                                                             aaGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCgtt                                                                                                                    <  8:336518/94‑3 (MQ=255)
                                                                             aaGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGAt                                                                                                                    >  7:336518/1‑94 (MQ=255)
                                                                                gTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTAATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCt                                                          <  1:165063/149‑1 (MQ=255)
                                                                                    cAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCtttt                                                       <  3:85339/148‑1 (MQ=255)
                                                                                     aGGGCGCGCTGGGATTTGTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGGCCGATGACTTTATTCATTCAATGCGGCTTTTGt                                                     <  1:13263/149‑1 (MQ=255)
                                                                                      gggCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTa                                                    <  3:36250/149‑1 (MQ=255)
                                                                                        gcgcgcTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAgg                                                  <  3:333893/149‑1 (MQ=255)
                                                                                          gcgcTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGa                                                <  2:98736/149‑1 (MQ=255)
                                                                                                ggATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAaca                                            <  4:243467/147‑1 (MQ=255)
                                                                                                ggATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACt                                          <  2:277234/149‑1 (MQ=255)
                                                                                                       tcACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGCCTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGa                                    >  7:329495/1‑148 (MQ=255)
                                                                                                        cacaAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGaa                                   >  7:43338/1‑148 (MQ=255)
                                                                                                         acaAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGaat                                  >  7:53205/1‑148 (MQ=255)
                                                                                                                gCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGtt                          <  4:116187/149‑1 (MQ=255)
                                                                                                                   cTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                   <  4:139338/121‑1 (MQ=255)
                                                                                                                   cTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                   >  3:139338/1‑121 (MQ=255)
                                                                                                                      aTGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAAt                     <  6:247291/148‑1 (MQ=255)
                                                                                                                       ttgCACGTTGACTATCCCGTTAGGGGGAGACTTTAGACTTCCGGCCCGCTATAGATCTGTCCTCATAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCTGCTTTTGTAGTGAAACACTCAATGTATATTCCGTTGGAAttt                   <  1:311614/147‑2 (MQ=255)
                                                                                                                               aaagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                   >  2:243776/5‑109 (MQ=255)
                                                                                                                               aaagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                   <  1:243776/105‑1 (MQ=255)
                                                                                                                                gACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACa               >  4:129013/1‑144 (MQ=255)
                                                                                                                                gACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACa               <  3:129013/144‑1 (MQ=255)
                                                                                                                                aagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGt           >  5:312906/4‑148 (MQ=255)
                                                                                                                                  ctctCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCGCGATAGATGTGGCCTCAAGACAGCGTTTTTTGGTCCGATGACGTTATTGATTGAATGCGGGTTTTGTAGGGAAAAGCTCAATGaataata                              >  8:189457/1‑126 (MQ=255)
                                                                                                                                  ctctACCCTTAGGGGGATACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGTCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATc                              <  7:189457/127‑1 (MQ=255)
                                                                                                                                   tctcCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGtcat        >  5:205728/1‑148 (MQ=255)
                                                                                                                                     tcCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGtca         >  3:3641/1‑145 (MQ=255)
                                                                                                                                     tcCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGtca         <  4:3641/145‑1 (MQ=255)
                                                                                                                                       cccTTATGGGGTGACTTTAGACTACCGGCCCCCGAAATATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATcagca   <  7:118743/149‑1 (MQ=255)
                                                                                                                                       cccTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATcag     <  8:43338/147‑1 (MQ=255)
                                                                                                                                        ccTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAAt                                      <  8:137931/113‑1 (MQ=255)
                                                                                                                                        ccTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAAt                                      >  7:137931/1‑113 (MQ=255)
                                                                                                                                         cTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATcagcag  >  8:233233/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                               
CTATCGTGCTGTTTCAAGCAATCGATTTTGTCATTTTACGACCCCTCTCCAAGCCCC‑‑ATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAG  >  NC_002947/3162388‑3162670

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