Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,167,904 A→T 5.3% P328P (CCA→CCT PP_5407 → TnsD‑like transposition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,167,9040AT 5.3% 228.0 / 2.6 75P328P (CCA→CCTPP_5407TnsD‑like transposition protein
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base T (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

CCAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGT  >  NC_002947/6167760‑6168050
                                                                                                                                                |                                                                                                                                                   
ccAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGa                                                                                                                                                 <  2:1371/149‑1 (MQ=255)
  aaaTGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAgg                                                                                                                                               >  1:339923/1‑149 (MQ=255)
  aaaTGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAgg                                                                                                                                               >  3:95522/1‑149 (MQ=255)
    aTGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCt                                                                                                                                             <  4:336766/149‑1 (MQ=255)
       aCCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTaaa                                                                                                                                         <  6:371410/149‑1 (MQ=255)
        cccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAAc                                                                                                                                        <  8:178035/149‑1 (MQ=255)
         ccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACt                                                                                                                                       >  1:164517/1‑149 (MQ=255)
         ccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACt                                                                                                                                       >  4:141419/1‑149 (MQ=255)
         ccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCAGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACt                                                                                                                                       >  2:316214/1‑149 (MQ=255)
           gCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCGCTGGGAGGCTagtcact                                                                                                                                      <  2:313919/149‑7 (MQ=255)
            cAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaa                                                                                                                                    <  2:268251/149‑1 (MQ=255)
             aaGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaac                                                                                                                                   <  5:1702/149‑1 (MQ=255)
                ccccGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaactaa                                                                                                                                <  2:197439/149‑1 (MQ=255)
                  ccGTTGTGACTGCCATCCGTTGAAGCTTCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAaca                                                                                                                              <  8:247650/149‑1 (MQ=255)
                            tGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGTCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTaa                                                                                                                    <  3:235033/149‑1 (MQ=255)
                              ccATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCGGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaact                                                                                                                                  >  7:171402/1‑133 (MQ=255)
                              ccATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCGGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaact                                                                                                                                  <  8:171402/133‑1 (MQ=255)
                              ccATCCGTTGAAGCATCTGAGACTGATTTGGAGGCTATTCGGGAGCTTAGAATCATTCGTAAAGGCCTATCAACATTCCGCGAACAAGCAGCAACTTCCAGTGACTAGCCTCACAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGc                                                                                                                  >  5:259771/1‑149 (MQ=255)
                               cATCCGTTTAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTATAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGc                                                                                                                  <  7:164267/148‑1 (MQ=255)
                                aTCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGcc                                                                                                                 >  7:300524/1‑148 (MQ=255)
                                        aaGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGc                                                                                                        >  3:21778/1‑149 (MQ=255)
                                          gCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCAATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGACTATTTCAGcc                                                                                                       >  3:85445/1‑148 (MQ=255)
                                               tGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaagg                                                                                                  >  2:311015/1‑148 (MQ=255)
                                                 aGTCTGATTTGGTTGCTATTCTGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCATCGAGGCTAAACTTAACTAACACTTTCTTTAATCCTATTTCAGCCCaagga                                                                                                 <  8:114125/147‑1 (MQ=255)
                                                   tCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAaggaagga                                                                                             <  8:300524/149‑1 (MQ=255)
                                                              ggCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGcc                                                                                  >  4:65080/1‑149 (MQ=255)
                                                              ggCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGAAGCAACTTACAGTGACTAGCCTCCCAGCGAGGCTAAACTTTACTAAAACTTTATTTAAGCTTATTTCAGACAAAGGAAGGAGAACTGGTgac                                                                                  >  7:39080/1‑147 (MQ=255)
                                                               gCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGccc                                                                                 >  5:244488/1‑149 (MQ=255)
                                                               gCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGccc                                                                                 >  1:47208/1‑149 (MQ=255)
                                                               gCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGccc                                                                                 >  1:52564/1‑149 (MQ=255)
                                                                 tATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGcccc                                                                                >  4:217153/1‑148 (MQ=255)
                                                                     cGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            <  4:95522/148‑1 (MQ=255)
                                                                      gggAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAAcc                                                                          <  5:269592/149‑1 (MQ=255)
                                                                       ggAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAAccc                                                                         <  2:164517/149‑1 (MQ=255)
                                                                               gAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACt                                                                 <  3:141419/149‑1 (MQ=255)
                                                                                aaTCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            <  3:365934/137‑1 (MQ=255)
                                                                                aaTCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTa                                                                <  7:247330/149‑1 (MQ=255)
                                                                                aaTCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGATAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            >  4:365934/1‑137 (MQ=255)
                                                                                 ataaTTAGTTAAGGCATATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAt                                                               <  6:259771/146‑1 (MQ=255)
                                                                                       cGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAAt                                                          <  4:85445/148‑1 (MQ=255)
                                                                                            aGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTg                                                    <  3:65080/149‑1 (MQ=255)
                                                                                             ggCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaactaac                                                                                                                               >  7:163555/1‑73 (MQ=255)
                                                                                             ggCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaactaac                                                                                                                               <  8:163555/73‑1 (MQ=255)
                                                                                               ccTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAAc                                                                                         >  4:357599/1‑109 (MQ=255)
                                                                                                  aTCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAGTTAATGACTTGAAg                                                 >  2:176191/1‑146 (MQ=255)
                                                                                                     aaCATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGaa                                            >  8:72807/1‑148 (MQ=255)
                                                                                                        aTGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtc                                        <  2:339923/149‑1 (MQ=255)
                                                                                                            cgcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTATAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctc                                     >  1:81941/1‑149 (MQ=255)
                                                                                                            cgcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAt                                                               <  3:147612/122‑1 (MQ=255)
                                                                                                            cgcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAt                                                               >  4:147612/1‑122 (MQ=255)
                                                                                                            cgcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCATATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctcc                                    <  7:229674/149‑1 (MQ=255)
                                                                                                            cccctaCCAGCATCAACATACAGCGCCTATCCCCCCACCGAGGCTAAAAATAAATAACACTTTCATTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctcc                                    <  4:202701/144‑1 (MQ=255)
                                                                                                             gcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCcg                                   >  3:180272/1‑149 (MQ=255)
                                                                                                                   cagcagCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACt                             >  7:141761/1‑149 (MQ=255)
                                                                                                                   cagcagCAACTTCCAGTGACTAGCCTCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACt                             >  2:88923/1‑148 (MQ=255)
                                                                                                                      cagcaACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGATACCCAAAAAACTATTTAATGACTTGAAGGAGCAAATCGTCCGCGAACtgct                          >  4:233725/1‑147 (MQ=255)
                                                                                                                      cagcaACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACtgat                          <  7:72807/149‑1 (MQ=255)
                                                                                                                      cagcaACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACtgat                          <  4:21778/149‑1 (MQ=255)
                                                                                                                         caACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGAc                       >  1:359369/1‑149 (MQ=255)
                                                                                                                            cTTCCAGTGACTAGTCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTgg                    <  1:311015/149‑1 (MQ=255)
                                                                                                                                ttaagtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgcgc                 >  7:59966/8‑148 (MQ=255)
                                                                                                                                aattGACTAGCCTACCAGCTAGTCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAATGAGAACTTGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAATGAGAAAATCCTCCGCGAACTTATGACTGgcgc                 <  8:39080/145‑1 (MQ=255)
                                                                                                                                    gtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGc             >  5:51875/4‑148 (MQ=255)
                                                                                                                                      cTAGCCTCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaaa         <  1:88923/149‑1 (MQ=255)
                                                                                                                                       tAGCCTCCCAGCGAGGCTAAACTTAAGTAAGACTTTCTTTAAGCCTATTTCAGCCCATGGAAGGAGAACTGGTGCCCCGATACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaat          >  5:354827/1‑147 (MQ=255)
                                                                                                                                       tAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaa          >  2:285065/1‑148 (MQ=255)
                                                                                                                                         gCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaaaaa       <  6:51875/149‑1 (MQ=255)
                                                                                                                                          ccTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAg      <  3:217153/149‑1 (MQ=255)
                                                                                                                                            tCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAACGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaaaaa       <  4:71373/146‑1 (MQ=255)
                                                                                                                                            tCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAAGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaaaaa       >  3:71373/1‑146 (MQ=255)
                                                                                                                                             cccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACg   >  6:170407/1‑149 (MQ=255)
                                                                                                                                              ccAGCGAGGCTAAACTTAACTACCACTTTCTTTAAGCCTATTTCAGACCAAGGATGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATGCTCCGCTAACTGATGACTGTCGCCAGCAAAAAAGCCGt  >  6:89710/1‑149 (MQ=255)
                                                                                                                                              ccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTAAATTACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGAc    <  1:316214/147‑1 (MQ=255)
                                                                                                                                               cAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCa               <  6:88678/135‑1 (MQ=255)
                                                                                                                                               cAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCa               >  5:88678/1‑135 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGT  >  NC_002947/6167760‑6168050

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