Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,857,538:1 +A 5.2% coding (72/231 nt) PP_3407 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,857,5381.A 5.2% 306.1 / 10.1 77coding (72/231 nt)PP_3407hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (36/37);  new base A (2/2);  total (38/39)
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.82e-01

AAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTGCTGCCC  >  NC_002947/3857390‑3857685
                                                                                                                                                     |                                                                                                                                                   
aaaCTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCCGGCCTGAGCTTGTCGTACCTGCTGTT‑c                                                                                                                                                    <  1:913791/150‑1 (MQ=255)
aaaCTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑c                                                                                                                                                    <  2:1508689/150‑1 (MQ=255)
  aCTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGa                                                                                                                                                  <  1:777292/150‑1 (MQ=255)
   cTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGAc                                                                                                                                                 <  1:993350/150‑1 (MQ=255)
   cTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCCACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGAc                                                                                                                                                 <  1:1365257/150‑1 (MQ=255)
    tGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACa                                                                                                                                                <  1:1703032/150‑1 (MQ=255)
    tGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCATGGGCCTGTCGTACCTGCTGTT‑CGACa                                                                                                                                                <  1:1202718/150‑1 (MQ=255)
       gggCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                             <  1:622327/150‑1 (MQ=255)
       gggCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                             <  1:1670888/150‑1 (MQ=255)
        tccggccctgctgaagggcctgggacgcaccggtttcccgtttcaccccctgaggggcacctgaacgcgaggcataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                             <  1:468408‑M1/29‑1 (MQ=255)
        tccggccctgctgaagggcctgggacgcaccggtttcccgtttcaccccctgaggggcacctgaacgcgaggcataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                             >  2:468408‑M1/121‑149 (MQ=255)
        ggCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCt                                                                                                                                            >  1:1625807/1‑150 (MQ=255)
          ccGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCg                                                                                                                                          >  2:12431/1‑150 (MQ=255)
           cGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGc                                                                                                                                         >  1:913758/1‑150 (MQ=255)
               gcgcAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGc                                                                                                                                     >  2:1430579/1‑150 (MQ=255)
                cgcAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCt                                                                                                                                    <  2:116118/150‑1 (MQ=255)
                cgcAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCGGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCt                                                                                                                                    <  1:1314603/150‑1 (MQ=255)
                  cAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGt                                                                                                                                  <  2:652405/150‑1 (MQ=255)
                    ggATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAGGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTAt                                                                                                                                >  1:1275624/1‑150 (MQ=255)
                        gAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACa                                                                                                                            <  1:1430579/150‑1 (MQ=255)
                             ggTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGCAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCGGTATGACATTGAc                                                                                                                       >  2:1032020/1‑150 (MQ=255)
                             ggTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGAc                                                                                                                       >  1:44860/1‑150 (MQ=255)
                              gTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACg                                                                                                                      <  2:609389/150‑1 (MQ=255)
                                 cgcTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCa                                                                                                                   >  1:1776890/1‑150 (MQ=255)
                                    tACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCaa                                                                                                                  >  1:1350298/1‑148 (MQ=255)
                                    tACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCaa                                                                                                                  <  2:1350298/148‑1 (MQ=255)
                                                   gACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCtctc                                                                                                 >  1:1207999/1‑150 (MQ=255)
                                                   gACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCtctc                                                                                                 >  1:1264705/1‑150 (MQ=255)
                                                   gACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTCCTGTT‑CGACAGCCTCGCACGCTGTATCACATTGACTGAAAGCCCGCCCCATGCtctc                                                                                                 >  2:344448/1‑150 (MQ=255)
                                                    aCAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑c                                                                                                                                                    <  1:1193644/98‑1 (MQ=255)
                                                    aCAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑c                                                                                                                                                    >  2:1193644/1‑98 (MQ=255)
                                                        gCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCACCCCATGCTCTCCTTAt                                                                                            <  2:77869/150‑1 (MQ=255)
                                                         cTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATg                                                                                           >  1:1007754/1‑150 (MQ=255)
                                                             cggcATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTc                                                                                       <  2:1007754/150‑1 (MQ=255)
                                                                    cGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCa                                                                                <  1:866741/150‑1 (MQ=255)
                                                                     gCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTGTGGCATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAg                                                                               <  2:739195/150‑1 (MQ=255)
                                                                     gCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAg                                                                               <  1:1417812/150‑1 (MQ=255)
                                                                      cccGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAgg                                                                              >  1:32676/1‑150 (MQ=255)
                                                                      cccGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAgg                                                                              >  1:1695577/1‑150 (MQ=255)
                                                                       ccGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGCAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGt                                                                             <  2:1176102/150‑1 (MQ=255)
                                                                            ggTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGcccc                                                                        >  2:580250/1‑150 (MQ=255)
                                                                             gTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCt                                                                       >  2:131596/1‑150 (MQ=255)
                                                                             gTGATCGCCCGGTTTCTGCACTTCTGCCTGCTAGAAGAAGACTTCATTCCTGAGCCTGTCGTACCTGCTGGC‑CGACAGCCTCGAACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCt                                                                       >  1:1077844/1‑150 (MQ=255)
                                                                                aTCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGACGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAg                                                                    <  2:1335398/150‑1 (MQ=255)
                                                                                 tagcaccgTTGCTGCGCTTATGCCTGATGGAAGAAGACTTCAGGTCTGCGCATTTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTTAgg                                                                   <  2:447814/143‑1 (MQ=255)
                                                                                 tCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAgg                                                                   >  1:61197/1‑150 (MQ=255)
                                                                                 tCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAgg                                                                   <  1:1797103/150‑1 (MQ=255)
                                                                                 ataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTc                                                                                                                                           <  1:1410699‑M1/31‑1 (MQ=255)
                                                                                 ataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTc                                                                                                                                           >  2:1410699‑M1/48‑78 (MQ=255)
                                                                                      cGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTg                                                              >  2:1803393/1‑150 (MQ=255)
                                                                                      cGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTg                                                              >  1:464082/1‑150 (MQ=255)
                                                                                            cTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACgg                                                        <  1:1032020/150‑1 (MQ=255)
                                                                                            cTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACgg                                                        <  1:580250/150‑1 (MQ=255)
                                                                                             tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                       <  2:913758/150‑1 (MQ=255)
                                                                                             tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                       >  2:577731/1‑150 (MQ=255)
                                                                                             tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                       <  1:1527958/150‑1 (MQ=255)
                                                                                             tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                       <  2:1776890/150‑1 (MQ=255)
                                                                                                  ttCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAAcc                                                  >  1:1420989/1‑150 (MQ=255)
                                                                                                     tgcctgcTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               <  2:61197/150‑1 (MQ=255)
                                                                                                     tgcctgcTGGAAGAAGACTTCAGGCCTGAGCCTGTAGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGGAACGGGAAACCGGt                                               >  1:238593/1‑150 (MQ=255)
                                                                                                        ctgctgGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGt                                                                             >  1:359693/1‑117 (MQ=255)
                                                                                                        ctgctgGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGt                                                                             <  2:359693/117‑1 (MQ=255)
                                                                                                         ggcgtcaatgtcatacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGGCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               >  2:746488‑M1/24‑146 (MQ=255)
                                                                                                         ggcgtcaatgtcatacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGGCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               <  1:746488‑M1/123‑1 (MQ=255)
                                                                                                            tgGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTccc                                        >  1:1449970/1‑150 (MQ=255)
                                                                                                             ggAAGAAGACTTCAGGGCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCa                                       <  2:95046/150‑1 (MQ=255)
                                                                                                             ggAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCAAACCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCa                                       >  1:1060707/1‑150 (MQ=255)
                                                                                                              gaagaaGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAg                                      >  2:646765/1‑150 (MQ=255)
                                                                                                               aagaagACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAgg                                     <  2:1449970/150‑1 (MQ=255)
                                                                                                                 gaagaCTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               >  2:64146/1‑138 (MQ=255)
                                                                                                                 gaagaCTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               <  1:64146/138‑1 (MQ=255)
                                                                                                                 gaagaCTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGcc                                   >  2:1101957/1‑150 (MQ=255)
                                                                                                                      tacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               <  1:172632‑M1/123‑1 (MQ=255)
                                                                                                                      tacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                               >  2:172632‑M1/11‑133 (MQ=255)
                                                                                                                      cTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGTTTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTCca                              <  1:1040276/150‑1 (MQ=255)
                                                                                                                           ggCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAgg                         <  2:44860/150‑1 (MQ=255)
                                                                                                                              cTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGcc                      >  1:900062/1‑150 (MQ=255)
                                                                                                                                   ccTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACa                 >  1:862197/1‑150 (MQ=255)
                                                                                                                                            ccTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGtg        <  1:12431/150‑1 (MQ=255)
                                                                                                                                                ctgTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTCCGTCCCAGGCCCTTCAGCAGGGCCGGCCAATGCCGGtgctgc    >  2:365000/1‑150 (MQ=255)
                                                                                                                                                  gTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTGCTGccc  >  2:1719342/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTGCTGCCC  >  NC_002947/3857390‑3857685

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