Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,289,861:1 +A 5.6% coding (1128/2349 nt) PP_1126 → hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,289,8611.A 5.6% 362.3 / 14.3 89coding (1128/2349 nt)PP_1126hydrolase
Reads supporting (aligned to +/- strand):  ref base . (47/37);  new base A (2/3);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGC  >  NC_002947/1289714‑1290005
                                                                                                                                                    |                                                                                                                                                
gTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAATGCATGTTCAACGCCTTGCCCGCTGAGGCCGATGCACTGGCCGCTTTCTATGGCGCAAGCAACAACGACGATGGCAATGCCtc                                                                                                                                                <  1:574645/150‑1 (MQ=255)
 ttCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGACCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGccatcg                                                                                                                                               >  2:1391907/1‑146 (MQ=255)
   cTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCtcgtc                                                                                                                                             <  1:1155229/150‑1 (MQ=255)
    tATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAGCGACGATTGCAATGCCtcgtcg                                                                                                                                            >  2:103937/1‑150 (MQ=255)
      ttCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCa                                                                                                                                          >  2:714028/1‑150 (MQ=255)
       tCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCaa                                                                                                                                         <  1:1700229/150‑1 (MQ=255)
       tCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCAGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGACAGCAACAACGACGATGGCAATGCCTCGTCGCaa                                                                                                                                         <  1:2232187/150‑1 (MQ=255)
       tCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGACCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcgttg                                                                                                                                         <  1:1391907/150‑4 (MQ=255)
         ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       >  2:1747064/1‑150 (MQ=255)
         ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       <  1:2367397/150‑1 (MQ=255)
         ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACATATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       >  1:2108697/1‑150 (MQ=255)
           tGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGgc                                                                                                                                     <  1:1453598/150‑1 (MQ=255)
                gTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAAt                                                                                                                                >  1:1927289/1‑150 (MQ=255)
                   aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTTCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTc                                                                                                                             >  2:583474/1‑150 (MQ=255)
                   aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcgttgtt                                                                                                                                       >  2:1315591/1‑135 (MQ=255)
                   aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcgttgtt                                                                                                                                       <  1:1315591/140‑6 (MQ=255)
                         gTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGTGCGCAATTTCAAGCAg                                                                                                                       <  1:1340222/150‑1 (MQ=255)
                         gTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAg                                                                                                                       <  1:169171/150‑1 (MQ=255)
                         gTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAg                                                                                                                       <  1:1857926/150‑1 (MQ=255)
                              ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcg                                                                                                                  >  1:1076588/1‑150 (MQ=255)
                              ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcg                                                                                                                  >  1:1544607/1‑150 (MQ=255)
                              ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcg                                                                                                                                            <  1:1681186/124‑1 (MQ=255)
                              ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcg                                                                                                                                            >  2:1681186/1‑124 (MQ=255)
                                       ggggcgGACCTATGTCGCCGGGCGCGAGGGCATGTGCAAGGCCTGGCTCGCGGAGGGCGATGCACTGGCCGCTTTCTAGGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGa                                                                                                         <  1:1083800/149‑1 (MQ=255)
                                              aCCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGATGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAc                                                                                                  <  2:735038/150‑1 (MQ=255)
                                              aCCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAc                                                                                                  >  2:1574346/1‑150 (MQ=255)
                                               ccTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACCTCTAcg                                                                                                 >  1:2199064/1‑150 (MQ=255)
                                                cTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAcgc                                                                                                >  1:2082461/1‑150 (MQ=255)
                                                        ccGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAcc                                                                                        <  1:427489/150‑1 (MQ=255)
                                                             cgcgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAg                                                                                   <  1:999203/150‑1 (MQ=255)
                                                              gcgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGAc                                                                                                        <  2:750444/128‑1 (MQ=255)
                                                              gcgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGAc                                                                                                        >  1:750444/1‑128 (MQ=255)
                                                                gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       <  2:1032104/95‑1 (MQ=255)
                                                                gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       >  1:1032104/1‑95 (MQ=255)
                                                                gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGt                                                                                >  2:438254/1‑150 (MQ=255)
                                                                  aGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTAc                                                                              <  1:2059550/150‑1 (MQ=255)
                                                                      aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTAc                                                                              >  1:1958832/1‑146 (MQ=255)
                                                                      aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTAc                                                                              <  2:1958832/146‑1 (MQ=255)
                                                                      aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                          >  2:2184033/1‑150 (MQ=255)
                                                                      aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                          >  2:2302964/1‑150 (MQ=255)
                                                                      aTGTTCAAGGCCTGGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAGGCCGCCAACCTGCAGCGTACCTAt                                                                          >  2:729930/1‑150 (MQ=255)
                                                                       tGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATg                                                                         <  1:338128/150‑1 (MQ=255)
                                                                        gTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATgg                                                                        >  2:1722303/1‑150 (MQ=255)
                                                                         ttCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGc                                                                       >  2:1684245/1‑150 (MQ=255)
                                                                         ttCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGc                                                                       >  1:307550/1‑150 (MQ=255)
                                                                            aaGGCCTTGCCCGCTGAGGGCGCTGCGCGGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaaa                                                                    >  1:931642/1‑150 (MQ=255)
                                                                            aaGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaaa                                                                    >  2:1963677/1‑150 (MQ=255)
                                                                             aGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAc                                                                   <  1:321052/150‑1 (MQ=255)
                                                                             aGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGAGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAc                                                                   >  2:1474516/1‑150 (MQ=255)
                                                                               gCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAc                                                                                                  >  1:1245258/1‑117 (MQ=255)
                                                                               gCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAc                                                                                                  <  2:1245258/117‑1 (MQ=255)
                                                                                  ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCCATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGc                                                              >  2:1783199/1‑150 (MQ=255)
                                                                                  ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGc                                                              >  1:48987/1‑150 (MQ=255)
                                                                                  ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGc                                                              <  1:729930/150‑1 (MQ=255)
                                                                                  ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGc                                                              >  2:684628/1‑150 (MQ=255)
                                                                                    gCCCGCTGAGGGCGATGTACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCg                                                            >  2:1229613/1‑150 (MQ=255)
                                                                                      ccGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGc                                                          >  1:1724084/1‑150 (MQ=255)
                                                                                       cGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGcc                                                         >  2:2347595/1‑150 (MQ=255)
                                                                                        gCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGccc                                                        <  1:2347595/150‑1 (MQ=255)
                                                                                           gAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGAGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCa                                                     <  1:1474516/150‑1 (MQ=255)
                                                                                              ggNGATGCACTGGCCGCTTTCTATGGCNAAAGCAACAACGNCGATGGNANTGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTg                                                  >  1:1435476/1‑150 (MQ=255)
                                                                                              ggCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTg                                                  >  1:1639233/1‑150 (MQ=255)
                                                                                                cGATGAGCTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                                <  1:714028/150‑1 (MQ=255)
                                                                                                   tGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAggg                                             >  1:359724/1‑150 (MQ=255)
                                                                                                    gCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgc                                            >  2:2401003/1‑150 (MQ=255)
                                                                                                     cACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAc                                                                                         <  2:1668252/104‑1 (MQ=255)
                                                                                                     cACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAc                                                                                         >  1:1668252/1‑104 (MQ=255)
                                                                                                      aCTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgcca                                          <  1:1574346/150‑1 (MQ=255)
                                                                                                       cTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccag                                         >  2:1198430/1‑150 (MQ=255)
                                                                                                           ccGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTTCCTATGGCAAACCCTGCCGGCCCGCGATGGAGGGccagccag                                     >  2:738858/1‑150 (MQ=255)
                                                                                                           ccGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGccagccag                                     >  2:1992344/1‑150 (MQ=255)
                                                                                                            cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagcca                                      <  2:311130/148‑1 (MQ=255)
                                                                                                            cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagcca                                      >  1:311130/1‑148 (MQ=255)
                                                                                                            cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                    <  1:583474/150‑1 (MQ=255)
                                                                                                             ggtttCTATGTAGAAATCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCagccagcc                                   <  2:1267186/148‑1 (MQ=255)
                                                                                                                ttCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGAt                                <  2:946927/150‑1 (MQ=255)
                                                                                                                 tCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGAtt                               <  2:2108697/150‑1 (MQ=255)
                                                                                                                 tCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGAtt                               <  2:1927289/150‑1 (MQ=255)
                                                                                                                   tATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgc                             >  1:339648/1‑150 (MQ=255)
                                                                                                                             gCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCgg                   >  1:2141398/1‑150 (MQ=255)
                                                                                                                               aacaacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaaa                                                                    <  1:1460778/99‑1 (MQ=255)
                                                                                                                               aacaacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaaa                                                                    >  2:1460778/1‑99 (MQ=255)
                                                                                                                               aacaacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTg                                                  >  2:603063/1‑117 (MQ=255)
                                                                                                                               aacaacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTg                                                  <  1:603063/117‑1 (MQ=255)
                                                                                                                                   acgacgATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGa                                 >  2:1285278/1‑130 (MQ=255)
                                                                                                                                   acgacgATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGa                                 <  1:1285278/130‑1 (MQ=255)
                                                                                                                                   acgacgATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTgg             <  1:1229613/150‑1 (MQ=255)
                                                                                                                                    cgacgaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGa            <  2:2141398/150‑1 (MQ=255)
                                                                                                                                        gaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc        >  2:7161/1‑150 (MQ=255)
                                                                                                                                        gaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc        >  2:438041/1‑150 (MQ=255)
                                                                                                                                          tGGCAATGCCTCGTGGCAACGGGGCAGTCCCAAGCCGGTGGGCTGGGTGGATGTCTACGCCGCGAATCTGCAGCGGACCTATGGCAAACCCGGCCGGCCCGCAGTGGAGGGCCAGACCGCCGATTGCCCGCTGCCGGTGCTGGATGCCGa      <  2:759974/150‑1 (MQ=255)
                                                                                                                                             cAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCg   >  2:1057312/1‑150 (MQ=255)
                                                                                                                                              aaTGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGc  <  2:307550/150‑1 (MQ=255)
                                                                                                                                                  ccTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc        >  1:2138578/1‑140 (MQ=255)
                                                                                                                                                  ccTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc        <  2:2138578/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGC  >  NC_002947/1289714‑1290005

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: