Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,999,708 G→C 5.3% Q131H (CAG→CAC)  PP_3528 → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,999,7080GC 5.3% 246.4 / 2.6 76Q131H (CAG→CAC) PP_3528ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTGGAGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACATTACCCCCGCATACGTCAACAATGCCGC  >  NC_002947/3999562‑3999854
                                                                                                                                                  |                                                                                                                                                  
cTGGAGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGa                                                                                                                                                  <  1:42951/149‑1 (MQ=255)
   gAGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGacca                                                                                                                                               >  7:10747/1‑149 (MQ=255)
   gAGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGacca                                                                                                                                               >  2:287022/1‑149 (MQ=255)
   gAGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGacc                                                                                                                                                <  7:162691/148‑1 (MQ=255)
    aGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGaccac                                                                                                                                              >  2:340615/1‑149 (MQ=255)
     ggCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGaccacc                                                                                                                                             <  4:238218/149‑1 (MQ=255)
       cATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTg                                                                                                                                           <  5:47061/149‑1 (MQ=255)
        aTTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGTCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTgg                                                                                                                                          >  5:264832/1‑149 (MQ=255)
        aTTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTgg                                                                                                                                          >  2:189726/1‑149 (MQ=255)
        aTTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTgg                                                                                                                                          <  5:180471/149‑1 (MQ=255)
           gAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGc                                                                                                                                         >  5:292983/1‑147 (MQ=255)
           gAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGc                                                                                                                                         >  6:286465/1‑147 (MQ=255)
            aGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCa                                                                                                                                       <  5:354301/148‑1 (MQ=255)
              cACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCCTGCTCGCGCACGac                                                                                                                                                 <  7:270024/136‑1 (MQ=255)
              cACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCCTGCTCGCGCACGac                                                                                                                                                 >  8:270024/1‑136 (MQ=255)
                 cGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGc                                                                                                                                 <  5:387049/149‑1 (MQ=255)
                 cGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGc                                                                                                                                 <  4:40700/149‑1 (MQ=255)
                       cATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGAcc                                                                                                                           >  2:220939/1‑149 (MQ=255)
                        aTCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCt                                                                                                                          <  5:131868/149‑1 (MQ=255)
                            aCCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGa                                                                                                                        >  4:272097/1‑147 (MQ=255)
                                 ggCCACGCAGCCGATGACCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTCTCGTCTACCTTTCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGc                                                                                                                   >  6:161332/1‑147 (MQ=255)
                                   ccACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGc                                                                                                               >  7:94462/1‑149 (MQ=255)
                                        cAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGtgtg                                                                                                          <  1:340615/149‑1 (MQ=255)
                                          gCCGATGCCCCCTAGGTGTTCGCCCAGGCCGCAGGCAAGCCAGTGGTCTACCTGGATGCCGAACAACCCTACCCACGACGCATCGGCCTGGTCGTGCGCGAGCAGCACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTCGc                                                                                                        >  6:283132/1‑149 (MQ=255)
                                            cGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCa                                                                                                      <  6:292983/149‑1 (MQ=255)
                                                  ccccTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGcc                                                                                                >  4:130259/1‑149 (MQ=255)
                                                   cccTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCcg                                                                                               <  3:130259/149‑1 (MQ=255)
                                                    acTCGGTGTTCGCCCATGCCAAAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCAACCCACACCGCATCCGACTCGTCGTGCGCCAGCAGCACCACCTGGCCACTGAGCAGGACCTGAAAGGCAAGCGTGTGGCCACTGTCcgc                                                                                              >  4:95976/2‑149 (MQ=255)
                                                           gtTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAgg                                                                                                                    >  5:133381/1‑120 (MQ=255)
                                                           gtTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAgg                                                                                                                    <  6:133381/120‑1 (MQ=255)
                                                              cGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGTCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAAt                                                                                     <  6:264832/148‑1 (MQ=255)
                                                               gCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGc                                                                                   <  1:189726/149‑1 (MQ=255)
                                                                      ccGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTAc                                                                            <  5:286465/149‑1 (MQ=255)
                                                                       cGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTAcc                                                                           >  6:177515/1‑149 (MQ=255)
                                                                       cGCAGGCAAGCCTGTGGCCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGcc                                                                                  >  2:265920/1‑142 (MQ=255)
                                                                       cGCAGGCAAGCCTGTGGCCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCACGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGcc                                                                                  <  1:265920/142‑1 (MQ=255)
                                                                             cAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTgg                                                                     <  3:272097/149‑1 (MQ=255)
                                                                             cAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTgg                                                                     <  1:220939/149‑1 (MQ=255)
                                                                              aaGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGa                                                                                                                        <  7:226266/97‑1 (MQ=255)
                                                                              aaGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGa                                                                                                                        >  8:226266/1‑97 (MQ=255)
                                                                              aaGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGc                                                                    <  8:94462/149‑1 (MQ=255)
                                                                               aGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCCTGCTCGCGCACGaccagg                                                                                                                                             <  5:25571/75‑3 (MQ=37)
                                                                               aGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCCTGCTCGCGCACGaccagg                                                                                                                                             >  6:25571/1‑73 (MQ=37)
                                                                                      ggTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGcc                                                                                  <  8:1927/127‑1 (MQ=255)
                                                                                      ggTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGcc                                                                                  >  7:1927/1‑127 (MQ=255)
                                                                                      ggTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGcc                                                                                  >  7:7487/1‑127 (MQ=255)
                                                                                      ggTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGcc                                                                                  <  8:7487/127‑1 (MQ=255)
                                                                                        tCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTg                                                          >  8:307582/1‑149 (MQ=255)
                                                                                        tCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTg                                                          >  4:252863/1‑149 (MQ=255)
                                                                                            ccTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTAcc                                                                           >  4:345593/1‑128 (MQ=255)
                                                                                            ccTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTAcc                                                                           <  3:345593/128‑1 (MQ=255)
                                                                                            ccTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTgcag                                                       >  8:125869/1‑148 (MQ=255)
                                                                                            ccTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTgca                                                        <  3:121378/147‑1 (MQ=255)
                                                                                            ccTGGCGGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTgca                                                        <  8:36783/147‑1 (MQ=255)
                                                                                               ggCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGa                                                                                                                                                  >  3:75502/1‑54 (MQ=255)
                                                                                               ggCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGa                                                                                                                                                  <  4:75502/54‑1 (MQ=255)
                                                                                                gCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGc                                                  <  8:10747/149‑1 (MQ=255)
                                                                                                  tGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAg                                                                               >  1:149349/1‑118 (MQ=255)
                                                                                                  agCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAg                                                                               <  2:149349/117‑1 (MQ=255)
                                                                                                    ccGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTgg                                               <  7:307582/148‑1 (MQ=255)
                                                                                                      gAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGAc                                                                                                                            >  5:341971/1‑69 (MQ=255)
                                                                                                      gAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGAc                                                                                                                            <  6:341971/69‑1 (MQ=255)
                                                                                                      gAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGcc                                                                   >  3:266070/1‑126 (MQ=255)
                                                                                                      gAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGcc                                                                   <  4:266070/126‑1 (MQ=255)
                                                                                                       aacaacCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCt                                            >  5:285956/1‑148 (MQ=255)
                                                                                                        acaacCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCa                                          <  3:95976/149‑1 (MQ=255)
                                                                                                        acaacCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCa                                          <  7:125869/149‑1 (MQ=255)
                                                                                                               tACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTAc                                     >  1:359994/1‑147 (MQ=255)
                                                                                                                 cccACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAgg                                 <  6:51045/149‑1 (MQ=255)
                                                                                                                           aTCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGt                                        >  3:350720/1‑132 (MQ=255)
                                                                                                                           aTCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGt                                        <  4:350720/132‑1 (MQ=255)
                                                                                                                               gCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACATTACCCCCGCa                   >  7:316767/1‑149 (MQ=255)
                                                                                                                                ccTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACATTACCCCCGCa                   >  2:195679/1‑148 (MQ=255)
                                                                                                                                    gTCGTGGGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACATTACCCCCGCATACGt              <  3:42851/149‑1 (MQ=255)
                                                                                                                                           gcgAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCAGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACAATACCCCCGCATACGTCAACAAt       >  1:51077/1‑149 (MQ=255)
                                                                                                                                                 aGGACGAGCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCAGGGGCCGCGGCTGGATTGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACATTACCCCCGCTTACGTCAACAATgcctc  <  2:51077/148‑3 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGGAGGCATTGAGCACCGGCGCCATCGACCTCGGCCACGCAGCCGATGCCCCCTCGGTGTTCGCCCAGGCCGCAGGCAAGCCTGTGGTCTACCTGGCTGCCGAACAACCCTACCCACGCGGCATCGGCCTGGTCGTGCGCGAGCAGGACCACCTGGCCAGTGTGCAGGACCTGAAAGGCAAGCGTGTGGCCACGGGCCGCGGCTGGAATGCCCAGTACCTGCTGGCCGTAGCCCTGCAGCAGGCTGGGCTCAGTTACAAGGACATTACCCCCGCATACGTCAACAATGCCGC  >  NC_002947/3999562‑3999854

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