Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,414,110 A→G 6.4% K261E (AAA→GAA)  PP_5643 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,414,1100AG 6.4% 180.7 / 2.3 63K261E (AAA→GAA) PP_5643hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base G (2/2);  total (33/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 = 1.00e+00

ACCTAATAAATTTCTCGGGGTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATACCTTCACCACTATCTGA  >  NC_002947/4413964‑4414258
                                                                                                                                                  |                                                                                                                                                    
aCCTAATAAATTTCTCGGGGTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCaa                                                                                                                                                     >  2:94846/1‑148 (MQ=255)
    aataaATTTCTCGGGGTCTGGTGCTACTTGGGAAGATAAAGGGGTCATCACAGAGCATCAGCGGCTAGTGCAGCTTAAAACGAGTAGAGTCCTAGTATCGTCCGCAGGGTCTGCTCTTCGACAGGTCAGAGCGAAGGTAAGCAAAAcc                                                                                                                                                 >  4:289406/1‑148 (MQ=255)
        aaTTTCTCGGGGTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCcaca                                                                                                                                             <  5:315366/148‑1 (MQ=255)
         aTTTCTCGGGGTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACt                                                                                                                                           <  1:170817/149‑1 (MQ=255)
              tcGGGGTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTa                                                                                                                                      >  1:316945/1‑149 (MQ=255)
                ggggTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGTGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTAc                                                                                                                                     <  5:305731/148‑1 (MQ=255)
                         tGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATgc                                                                                                                           <  1:94846/149‑1 (MQ=255)
                             aCTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGt                                                                                                                       >  3:116982/1‑149 (MQ=255)
                              cTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTAcc                                                                                                                                    <  5:295626/135‑1 (MQ=255)
                              cTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTAcc                                                                                                                                    >  6:295626/1‑135 (MQ=255)
                              cTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGt                                                                                                                       <  2:316945/148‑1 (MQ=255)
                               ttGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGt                                                                                                                       >  7:162445/1‑147 (MQ=255)
                                tGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGt                                                                                                                       >  3:186927/1‑146 (MQ=255)
                                   gatgatAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGAGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCAAACTGAGTACCAGCACAAGCGCGTTAGAtt                                                                                                                 >  2:333556/1‑149 (MQ=255)
                                    atgatAAA‑GGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTa                                                                                                               >  4:261937/1‑149 (MQ=255)
                                          aaGGGGTCATCACAGAGCATCAGCGGCTTGTGCATCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACTTGCGCGTTAGATTTAGGACt                                                                                                          <  1:333556/149‑1 (MQ=255)
                                            ggggTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCTAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGa                                                                                                        <  1:82675/149‑1 (MQ=255)
                                             gggTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCTGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAg                                                                                                       <  3:80813/149‑1 (MQ=255)
                                             gggTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAg                                                                                                       >  6:324388/1‑149 (MQ=255)
                                             gggTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAg                                                                                                       <  3:261937/149‑1 (MQ=255)
                                              ggTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGAGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGt                                                                                                      >  2:115255/1‑149 (MQ=255)
                                                tcatcaCAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCaaaa                                                                                                                                                   <  7:102959/102‑1 (MQ=255)
                                                tcatcaCAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCaaaa                                                                                                                                                   >  8:102959/1‑102 (MQ=255)
                                                  atcaCAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCAgg                                                                                                  >  8:140825/1‑149 (MQ=255)
                                                    cacaGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGc                                                                                                 >  8:227379/1‑148 (MQ=255)
                                                      caGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGa                                                                                               >  8:6908/1‑148 (MQ=255)
                                                       agagCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGaaa                                                                                             <  8:162445/149‑1 (MQ=255)
                                                           cATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGAt                                                                                                                  >  2:243587/1‑124 (MQ=255)
                                                           cATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGAt                                                                                                                  <  1:243587/124‑1 (MQ=255)
                                                             tCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTa                                                                                          >  7:121076/1‑146 (MQ=255)
                                                                gCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTTCTCTTCGAAAGGTCCCAGCGAAGGTAAGCAAAACCCACACTGACTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGAcg                                                                                    >  3:343094/1‑149 (MQ=255)
                                                                 cGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGAcgc                                                                                   <  7:6908/149‑1 (MQ=255)
                                                                   gCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGa                                                                                 <  7:140825/149‑1 (MQ=255)
                                                                        tgCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTAtt                                                                            >  5:223023/1‑149 (MQ=255)
                                                                          cAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGa                                                                          >  6:238286/1‑149 (MQ=255)
                                                                              ttAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGacac                                                                       >  5:341889/1‑148 (MQ=255)
                                                                                   ccGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCAc                                                                 >  3:123360/1‑149 (MQ=255)
                                                                                          gCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGaaaa                                                          >  8:43633/1‑149 (MQ=255)
                                                                                           cGGCCTAGTAGGGTCAGGAGGGTCTGCTCTTCGACAGGTCCGAGCGAATGTAATCAAAACCCACACTGATTACCAGCACATGCGCGTTATATTTAGGTCTTAGTCCGGCGAAATTAAATACGCGACTATTGACACTATCACCGTAAAAc                                                         <  4:343094/149‑1 (MQ=255)
                                                                                           cGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCGAAACCCACACTGAGTACCa                                                                                                                                   <  7:16430/75‑1 (MQ=255)
                                                                                           cGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCGAAACCCACACTGAGTACCa                                                                                                                                   <  7:22042/75‑1 (MQ=255)
                                                                                           cGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCGAAACCCACACTGAGTACCa                                                                                                                                   >  8:16430/1‑75 (MQ=255)
                                                                                           cGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCGAAACCCACACTGAGTACCa                                                                                                                                   >  8:22042/1‑75 (MQ=255)
                                                                                             gCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTc                                                       >  2:222620/1‑149 (MQ=255)
                                                                                             gCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGGCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTc                                                       <  6:302997/149‑1 (MQ=255)
                                                                                                tagtagGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTa                                                                              <  1:318160/123‑1 (MQ=255)
                                                                                                tagtagGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTa                                                                              >  2:318160/1‑123 (MQ=255)
                                                                                                 agtagGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCg                                                      >  5:49044/1‑146 (MQ=255)
                                                                                                     gggTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCt                                               <  5:324388/149‑1 (MQ=255)
                                                                                                     gggTCAGAAGGTTCTGCTGTTCGCGAGGTCCGAGCGAATGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGGTAGATTTAGGACTGAGTCCGGCTAAAGTAAAGACGCGACTATTTACACTAGCACCTGAAAACTCGTATAGCt                                               <  6:252026/149‑1 (MQ=255)
                                                                                                        tCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGc                                                >  4:218220/1‑145 (MQ=255)
                                                                                                        tCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGc                                                <  3:218220/145‑1 (MQ=255)
                                                                                                        tCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGa                                            >  6:172217/1‑149 (MQ=255)
                                                                                                          agcagGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGAt                                           <  4:186927/148‑1 (MQ=255)
                                                                                                            cagGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGAt                                           >  7:31310/1‑146 (MQ=255)
                                                                                                            cagGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGAt                                           <  8:31310/146‑1 (MQ=255)
                                                                                                                     ctctTCGACAGGTCCGATCGAAGGTAAGCAAAACCCACACTGATTACCATCACATGCGCGTTAGATTTAGGACTGACTCAGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCaaa                               <  3:289406/149‑1 (MQ=255)
                                                                                                                     ctctTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCaa                                <  6:341889/148‑1 (MQ=255)
                                                                                                                           gACATGTCCGAGCGAAGGTAAGCAAAACCCACACTTAGTACAAGCAGATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTTATAATTACATCAAACGaca                          <  4:123360/148‑1 (MQ=255)
                                                                                                                            aCAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTa                                                    >  5:149971/1‑121 (MQ=255)
                                                                                                                            aCAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTa                                                    <  6:149971/121‑1 (MQ=255)
                                                                                                                                 tCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATa                   >  4:340069/1‑149 (MQ=255)
                                                                                                                                   cGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATAcc                 <  1:222620/149‑1 (MQ=255)
                                                                                                                                      gagAAGGTAAGCAAAAACCCCACTGAGTACCAGCACATGCGCGTTAGATTTAGTACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATACCt                <  8:121076/145‑1 (MQ=255)
                                                                                                                                        gAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATACCTTcac            >  4:328427/1‑149 (MQ=255)
                                                                                                                                          aGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATACCTTcacca          >  4:83223/1‑149 (MQ=255)
                                                                                                                                                  aaaaCCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATACCTTCACCACTATCTGa  <  3:83223/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCTAATAAATTTCTCGGGGTCTGGTGCTACTTGGGATGATAAAGGGGTCATCACAGAGCATCAGCGGCTTGTGCAGCTTAAACCGAGTAGCGGCCTAGTAGGGTCAGCAGGGTCTGCTCTTCGACAGGTCCGAGCGAAGGTAAGCAAAACCCACACTGAGTACCAGCACATGCGCGTTAGATTTAGGACTGAGTCCGGCGAAAGTAAAGACGCGACTATTGACACTAGCACCGGAAAACTCGTAGAGCTTGATAAGTACATCAAACGACACACAATACCTTCACCACTATCTGA  >  NC_002947/4413964‑4414258

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