Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,200,645 G→A 6.1% G131G (GGC→GGT)  PP_3690 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,200,6450GA 6.1% 196.2 / 2.0 65G131G (GGC→GGT) PP_3690hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (2/2);  total (35/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 = 9.06e-01

TGAAACATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCAAAATCGAAGAGC  >  NC_002947/4200501‑4200789
                                                                                                                                                |                                                                                                                                                 
tGAAACATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGc                                                                                                                                                  >  6:76091/1‑146 (MQ=255)
  aaaCATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCa                                                                                                                                                >  3:108291/1‑146 (MQ=255)
   aaCATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAgaga                                                                                                                                            >  4:303840/1‑149 (MQ=255)
    aCATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCACCAGAGAc                                                                                                                                           >  4:290134/1‑149 (MQ=255)
    aCATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCACCAGAGAc                                                                                                                                           >  4:294847/1‑149 (MQ=255)
     cATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACg                                                                                                                                          >  1:174471/1‑149 (MQ=255)
     cATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACg                                                                                                                                          >  4:298726/1‑149 (MQ=255)
        cGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATTATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGt                                                                                                                                       >  8:286915/1‑149 (MQ=255)
              cTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCga                                                                                                                                 <  3:303840/149‑1 (MQ=255)
                  cTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACt                                                                                                                             <  5:159830/149‑1 (MQ=255)
                     gTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCa                                                                                                                                                >  7:133307/1‑127 (MQ=255)
                     gTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCa                                                                                                                                                <  8:133307/127‑1 (MQ=255)
                     gTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACTGGTTTCCGAGACTg                                                                                                                            <  4:108291/147‑1 (MQ=255)
                                gtaggaCAGTCGAGTGGTTATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGAAATGATGGAGCTAATTCTGAAGATCGCCAGAGACGGGTTTCCGAGACTGGAGAGCattgatt                                                                                                               <  7:282600/147‑1 (MQ=255)
                                ggaggaCAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCattgatt                                                                                                               <  3:298726/149‑1 (MQ=255)
                                   ggaCAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCattgat                                                                                                                >  4:245734/1‑145 (MQ=255)
                                   ggaCAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCattgat                                                                                                                <  3:245734/145‑1 (MQ=255)
                                   ggaCAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCACCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAg                                                                                                            <  3:294847/149‑1 (MQ=255)
                                   ggaCAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCACCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAg                                                                                                            <  3:290134/149‑1 (MQ=255)
                                    gaCAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGc                                                                                                           >  8:58283/1‑149 (MQ=255)
                                        gTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGGCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTc                                                                                                       <  5:179746/149‑1 (MQ=255)
                                                gTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGc                                                                                               <  7:58283/149‑1 (MQ=255)
                                                 tcatcaACAAACTAAAGACCAGATACCTACAGTAAAAGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTATCAACCTCGCAAGAGACCGTTTTCCGAGAATGGAGAGCATTGATTCAGACTTAAACGACGCa                                                                                              >  5:100949/1‑149 (MQ=255)
                                                     gAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCattgat                                                                                                                <  6:285013/127‑1 (MQ=255)
                                                     gAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCattgat                                                                                                                >  5:285013/1‑127 (MQ=255)
                                                     gAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAAc                                                                                          >  1:93523/1‑149 (MQ=255)
                                                       aCAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACTACGCAGAACTc                                                                                        >  4:290578/1‑149 (MQ=255)
                                                           aCTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGttt                                                                                    <  5:76091/149‑1 (MQ=255)
                                                             tCAAGACCAGATTCGATCAGTAAACGAATGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCATCAACGACGCAGAACTCGTTTCg                                                                                  >  2:184612/1‑149 (MQ=255)
                                                                 gACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGgaga                                                                                                                        >  6:303024/1‑107 (MQ=255)
                                                                 gACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGgaga                                                                                                                        <  5:303024/107‑1 (MQ=255)
                                                                   ccAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACt                                                                            <  2:174471/149‑1 (MQ=255)
                                                                     aGATTCGTTCAGTAAACGATTGGATGGGAAGTGAATGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTAc                                                                          >  1:201389/1‑149 (MQ=255)
                                                                           gTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATTTAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCg                                                                    <  3:290578/149‑1 (MQ=255)
                                                                               aGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGttt                                                                                    <  2:50701/129‑1 (MQ=255)
                                                                               aGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGttt                                                                                    >  1:50701/1‑129 (MQ=255)
                                                                                gTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTa                                                                           >  2:287479/1‑137 (MQ=255)
                                                                                gTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTa                                                                           <  1:287479/137‑1 (MQ=255)
                                                                                                gAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACAGCCTccgcc                                                >  7:165891/1‑148 (MQ=255)
                                                                                                     gagtgTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCa                                           >  3:243368/1‑148 (MQ=255)
                                                                                                      tgtgtCGTAGATATAGTCTTGTTGGAGCGAATTCTGAAGCTGGCCACAGTCGGGTTTCGTAGTCTTGAGATCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTATCTACCGCCCGTAGTATGCACCGCCTCCTGCGTCCaaa                                         <  6:100949/148‑1 (MQ=255)
                                                                                                      agtgTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGc                                                                       <  5:78491/119‑1 (MQ=255)
                                                                                                      agtgTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGc                                                                       >  6:78491/1‑119 (MQ=255)
                                                                                                            gTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGc                                   >  8:192239/1‑149 (MQ=255)
                                                                                                                      aCATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATc                        <  4:243368/149‑1 (MQ=255)
                                                                                                                        atgatgGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGc                                       >  7:199840/1‑133 (MQ=255)
                                                                                                                        atgatgGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGc                                       <  8:199840/133‑1 (MQ=255)
                                                                                                                        atgatgGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTa                      >  1:22674/1‑149 (MQ=255)
                                                                                                                              gAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCt                >  3:78396/1‑149 (MQ=255)
                                                                                                                               aGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCt                >  3:112915/1‑148 (MQ=255)
                                                                                                                                gCGAATTCTGAAGCTCGCCATAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCa              <  7:286915/149‑1 (MQ=255)
                                                                                                                                gCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTg                                    >  1:92082/1‑128 (MQ=255)
                                                                                                                                gCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTg                                    <  2:92082/128‑1 (MQ=255)
                                                                                                                                  gAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCaa             >  6:138788/1‑148 (MQ=255)
                                                                                                                                   aaTTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCaaaa           >  2:192996/1‑149 (MQ=255)
                                                                                                                                     ttCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCa                                                                                              >  6:305036/1‑65 (MQ=255)
                                                                                                                                     ttCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCa                                                                                              <  5:305036/65‑1 (MQ=255)
                                                                                                                                     ttCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATc                        <  7:10502/134‑1 (MQ=255)
                                                                                                                                     ttCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATc                        >  8:10502/1‑134 (MQ=255)
                                                                                                                                        tGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑aaa                            <  3:253299/127‑1 (MQ=255)
                                                                                                                                        tGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑aaa                            >  4:253299/1‑127 (MQ=255)
                                                                                                                                        tGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCAAAATCg        >  4:73244/1‑147 (MQ=255)
                                                                                                                                          aaGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCAAAATCGaa      <  4:112915/147‑1 (MQ=255)
                                                                                                                                          aaGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCAAAATCGaa      <  3:73244/147‑1 (MQ=255)
                                                                                                                                           aGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCAAAATCGAagag   <  4:78396/149‑1 (MQ=255)
                                                                                                                                                gCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGTACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCTTGCCCGAAAAAATCTACCCCCTCAAAATCGAAGAGc  <  1:184612/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGAAACATCGCATCCTGACTCGTCCTTGTGAAGGAGGACAGTCGAGTGGTCATGAACAAACTCAAGACCAGATTCGTTCAGTAAACGATTGGATGGGAAGTGAGTGTCGTAGATATAGACATGATGGAGCGAATTCTGAAGCTCGCCAGAGACGGGTTTCCGAGACTGGAGAGCATTGATTGAGCCTTCAACGACGCAGAACTCGTTTCGGTAACTACCGCCCGTAGTATCCACCGCCTCCGCCGTCCAAAGCGTGCCCG‑AAAAATCTACCCCCTCAAAATCGAAGAGC  >  NC_002947/4200501‑4200789

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