Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,857,531 C→A 5.1% P22H (CCT→CAT)  PP_3407 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,857,5310CA 5.1% 252.9 / 2.4 78P22H (CCT→CAT) PP_3407hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base A (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTG  >  NC_002947/3857386‑3857679
                                                                                                                                                 |                                                                                                                                                     
ccAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTAcc                                                                                                                                                       >  1:90328/1‑146 (MQ=255)
ccAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTAcc                                                                                                                                                       <  2:90328/146‑1 (MQ=255)
 cAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACctgctg                                                                                                                                                  >  1:931289/1‑150 (MQ=255)
    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                                                                                                                                       >  2:468408‑M1/121‑149 (MQ=255)
            tccggccctgctgaagggcctgggacgcaccggtttcccgtttcaccccctgaggggcacctgaacgcgaggcataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                       <  1:468408‑M1/29‑1 (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                                                                                                            <  2:1350298/148‑1 (MQ=255)
                                        tACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCaa                                                                                                            >  1:1350298/1‑148 (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                                                        >  1:464082/1‑150 (MQ=255)
                                                                                          cGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTg                                                        >  2:1803393/1‑150 (MQ=255)
                                                                                                cTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACgg                                                  <  1:1032020/150‑1 (MQ=255)
                                                                                                cTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACgg                                                  <  1:580250/150‑1 (MQ=255)
                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                 <  1:1527958/150‑1 (MQ=255)
                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                 >  2:577731/1‑150 (MQ=255)
                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                 <  2:913758/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                                         >  2:172632‑M1/11‑133 (MQ=255)
                                                                                                                          tacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                         <  1:172632‑M1/123‑1 (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)
                                                                                                                                                 |                                                                                                                                                     
CCAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTG  >  NC_002947/3857386‑3857679

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