Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,459 T→G 22.0% L159R (CTG→CGG)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,4590TG22.0% 295.4 / 24.6 183L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (60/76);  new base G (32/12);  total (95/88)
Fisher's exact test for biased strand distribution p-value = 1.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

CGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4163310‑4163608
                                                                                                                                                     |                                                                                                                                                     
cGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCt                                                                                                                                                       <  2:1458051/150‑1 (MQ=255)
   gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcg                                                                                                                                                    <  2:3444678/150‑1 (MQ=255)
    tCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcgc                                                                                                                                                   <  2:2709065/150‑1 (MQ=255)
      ggTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCggcgcggcgccg                                                                                                                                                 >  2:4171334/1‑148 (MQ=255)
        tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                               >  2:3629161/1‑150 (MQ=255)
        tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                               <  2:1618698/150‑1 (MQ=255)
        tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcc                                                                                                                                               >  2:4774286/1‑149 (MQ=255)
           tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGCGCAGCGccc                                                                                                                                            >  1:267916/1‑150 (MQ=255)
           tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                            <  1:523342/150‑1 (MQ=255)
           tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                            <  2:6497709/150‑1 (MQ=255)
           tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGGTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGACTCGGCGCGGCGCCGCGccc                                                                                                                                            >  2:6840431/1‑150 (MQ=255)
             gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGTGCGCAGCGccccc                                                                                                                                          <  2:3505479/150‑1 (MQ=255)
             gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccc                                                                                                                                          >  2:850635/1‑150 (MQ=255)
              tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                         >  2:5255990/1‑150 (MQ=255)
              tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                         <  2:4944778/150‑1 (MQ=255)
                 tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGGGGAGCGcccccccaa                                                                                                                                      >  1:2318050/1‑147 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccccaaaa                                                                                                                                    >  1:3814296/1‑145 (MQ=255)
                   aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAaaaa                                                                                                                                    >  1:3635985/1‑150 (MQ=255)
                    ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATTAAAGTGTTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCGGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                   <  2:7376860/150‑1 (MQ=255)
                    ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCTGCGCCCCCCATaaaa                                                                                                                                   <  2:5731604/150‑1 (MQ=255)
                      aGTACCCTGCTGATCGGCCTGTGGTGTGGAATGAAAGTGTTCAAGCTCGACCGCTATACCGCCTTGCTGACGGCCGCCGGCAGCGCATTCTGCGGCGCCGCCGCCGTGCTGGATTTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAgc                                                                                                                                 <  2:6995330/150‑1 (MQ=255)
                          cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCAATCTGCGGCGCCGCCGCGGTGCTGGCTTTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCa                                                                                                                             <  2:1404965/150‑1 (MQ=255)
                           ccTGCTGATCGGCCTGTGGTGTGGAATGAAAGTGTTCAAGCTCGCCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAt                                                                                                                            <  2:5581521/150‑1 (MQ=255)
                            ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCATGCTCGACCGCGATACGGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATg                                                                                                                           <  2:3737344/150‑1 (MQ=255)
                            ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCCGCGCAGACCCCCCCAAAAAAAAGCCaag                                                                                                                           >  1:2924944/1‑148 (MQ=255)
                             tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCCCCGCGCCCCCCCAAAAAGCGcccggg                                                                                                                          >  2:5372753/1‑146 (MQ=255)
                                tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAg                                                                                                                       >  2:2430906/1‑150 (MQ=255)
                                tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAGGGccg                                                                                                                       >  2:6179152/1‑148 (MQ=255)
                                 gATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCGCCCAGGCAGt                                                                                                                      >  1:4298301/1‑150 (MQ=255)
                                  aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGGCGCCCCCCCCCCCAAAAAAGCGCCAAGGCAGTc                                                                                                                     >  1:456482/1‑150 (MQ=255)
                                      gCCTGGGGTGTGGCATGAATGTGTTCATGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                 <  2:2460703/150‑1 (MQ=255)
                                        cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGc                                                                                                               <  2:267916/150‑1 (MQ=255)
                                        cTGTGGTGTGGAATGAATGTTTTCACGCTCGACCGCTTTACCGCCTTGCTGACCGCCGCGGGCAGCGCAATCTGCGGCGCCGCCGCGGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGc                                                                                                               <  1:1806845/150‑1 (MQ=255)
                                           tggtgtggCATTAATGTTTACAAGCTCGACCGCGATACCGCCTGGCTGCCCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCGGTGCTGGAATTCTAATCGGCGCGGCCCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                            <  1:6707243/150‑1 (MQ=255)
                                           tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCCCCCGGGGCGTCGGCAGCGtc                                                                                                            >  1:4592167/1‑149 (MQ=255)
                                             gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGGAGCGTAGt                                                                                                          >  2:2177358/1‑150 (MQ=255)
                                              tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                         <  2:4925751/150‑1 (MQ=255)
                                                tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCCCCCCCCCAAAAAGCCCCCAGGCCGACGGAAGCGGAGTggg                                                                                                       >  1:5739799/1‑148 (MQ=255)
                                                 ggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCACAAAGAGCCAAGGCAGGCCGGAGCGTAGttgtg                                                                                                      >  1:318600/1‑146 (MQ=255)
                                                  gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCCAAAAAACCCCAGGGCCCTCGGGAGCCGAGGGCTGt                                                                                                     >  1:3714218/1‑150 (MQ=255)
                                                   cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                    >  2:5652358/1‑150 (MQ=255)
                                                    aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                   >  2:4236727/1‑150 (MQ=255)
                                                     tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                  >  2:5814297/1‑150 (MQ=255)
                                                      taatgttttCAAGCTCGACCGCTATACCGCTTTGCTGACCGCCGCCGGCAGCGCCATCGGGGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                 <  1:5510084/143‑1 (MQ=255)
                                                       aatttgtTCACGCCCGCCCGCGACCCGGCCTTGCGGCCCGCCGCGGGCCGCGCCAGGGGCGGCGCCGCCGCGGGGGGGGATTTCTCCGGGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                <  1:2190671/146‑1 (MQ=255)
                                                          gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGCGCCCCCCATAAAATCGCCCTGGCAGTCGGGAGCGTAGTGCTGGTAGGGAcc                                                                                             >  1:498929/1‑150 (MQ=255)
                                                          gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAATCGCCAAGGCAGTCGGCAGCGTAGTGGGGTTGGGTAcc                                                                                             >  1:983311/1‑150 (MQ=255)
                                                              tCACGCTCGACCGCGTTACCGCCTTGCTGACCGCCGCGGGCAGCGCCATGTGGGGCGCCGCCGCCGTGCTGGCATTCGACGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                         <  1:4237080/150‑1 (MQ=255)
                                                              tCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCCAAAAGCGCCAAGGCACTCGGAAGCCTAGTGGTGGTGGGGAccgagt                                                                                         >  1:4381733/1‑146 (MQ=255)
                                                               cacgCTCGACCGCTATCCCGCGTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                        <  2:2045234/147‑1 (MQ=255)
                                                               cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCACCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                        <  2:5110140/150‑1 (MQ=255)
                                                                aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAATCGTAGTGCTGTTCGGTACCTTGTCg                                                                                       >  2:5948265/1‑150 (MQ=255)
                                                                aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCCGCGCCCCCCCTAAAAGCGCCCAGGCAGGCGGAAGCGTAGTGCTGTTGGGTACCTGGTCg                                                                                       >  1:1139982/1‑150 (MQ=255)
                                                                    tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGTTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                   <  2:3709672/150‑1 (MQ=255)
                                                                    tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                   <  2:5359548/150‑1 (MQ=255)
                                                                    gcgcccGCGATCCGGCGTTGCTGCCCGCCGCGGGCAGCGCCTTGGGGGGCGCCGCCGCGGTGCGGGCATTCGAAGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                   <  1:6404586/146‑1 (MQ=255)
                                                                      gACCGCTATCCCGCTTTGCTGCCCGCGGCGGGCAGCGCCATCTGGGGCGCCGCCGCGGGGGGGGAATTCGACGCGGGGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                 <  1:85661/150‑1 (MQ=255)
                                                                               aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                                        <  2:2545310/150‑1 (MQ=255)
                                                                                 cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                      <  2:663643/150‑1 (MQ=255)
                                                                                  gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGCCCCCCCCAAAAACGCCCAGGCAGTCGGGAGCGTACTGGTGGTGGGGACCATGTGGATGTTGCTTTTACCAcag                                                                     >  1:3150584/1‑148 (MQ=255)
                                                                                   ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTgg                                                                    >  2:6121592/1‑150 (MQ=255)
                                                                                    cgtGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCGGGCATTCGACTCGGCGGTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc                                                                   <  2:4647703/148‑1 (MQ=255)
                                                                                      tGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCTTTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGa                                                                 <  2:4460377/150‑1 (MQ=255)
                                                                                         tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCa                                                              >  2:6864225/1‑150 (MQ=255)
                                                                                         tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCTAAAAGCGCCATGGCAGTCGGTAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATcc                                                              >  2:1289372/1‑149 (MQ=255)
                                                                                             cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                          >  2:7336254/1‑150 (MQ=255)
                                                                                                     gCAGCGCCATCTGCGGCGCCGCCGCCGGGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggc                                                  <  2:674418/150‑1 (MQ=255)
                                                                                                       aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCATTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATCAATGCCGgctggctg                                                >  2:1203185/1‑150 (MQ=255)
                                                                                                           ccATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATc                                            >  2:2993433/1‑150 (MQ=255)
                                                                                                             cccTGCGGCGCCGCCGCCGTGCTGGCATTCGACTGGGCGGTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                          <  1:778863/148‑1 (MQ=255)
                                                                                                             cacggcggcggcgCCGCAGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGCTGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                          >  2:5809631/5‑150 (MQ=255)
                                                                                                             aTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAAAAAAGCGCCATGGCCGTCGGAAGCGTAGTGCTGTTGGGGTCCATGTCGATGTTACTTTAACCACTGGCGATAAAATCCGGCTGGCTGCAACTg                                          >  1:1087008/1‑150 (MQ=255)
                                                                                                               cTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                        <  2:1276112/150‑1 (MQ=255)
                                                                                                                 gcggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCTTGTCGATGTTCCGTTACCCACGGGCGATCAATGCCGGCTGGCGGCATCAGGAtc                                      >  2:1506768/1‑149 (MQ=255)
                                                                                                                   ggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCATCAAAGCGCCCAGGCAGTCGGGCGCGGAGAGCAGTTGGGGGCCTTGTCGATGTTGCTATAACCACAGGCGCACAAAACCCGCAGGCTTCATCAGGATaaa                                    >  1:7294957/1‑148 (MQ=255)
                                                                                                                    gcgcCGCCGCGGTGCTGTCATTCTACTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCa                                   <  1:5874152/150‑1 (MQ=255)
                                                                                                                     cgccgccgcCGTGCTGGCTTTCGAATCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                                  <  1:7054345/150‑1 (MQ=255)
                                                                                                                       ccgccgcggTGCTGGAATTCGAATGGGGGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                <  1:4917008/150‑1 (MQ=255)
                                                                                                                       cccccgccgGGGTGGCATTCTACTGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                <  1:6849119/147‑1 (MQ=255)
                                                                                                                         gccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCAACTGGATACCATGGgc                              >  1:6077171/1‑150 (MQ=255)
                                                                                                                             ccgTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg                          <  2:2492150/150‑1 (MQ=255)
                                                                                                                                gggtGGCTTTCGACGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgc                       <  1:4577590/147‑1 (MQ=255)
                                                                                                                                      cATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgc                 <  2:5059128/150‑1 (MQ=255)
                                                                                                                                       aTTCGACTCGGCGCGGCGCAGCGCCCCCCCGAAAAGCGCCATGGCAGTCGGAAGCGTAGGGCGGGTCGGTACCGTGTCGATGTTCCGTTACCCACAGGCGATCAAAGCCGGCGGGCGGCATCTGGATACCATGGGCGCCGGgcggctgct                >  2:937915/1‑150 (MQ=255)
                                                                                                                                       aTTCGACCCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct                <  1:613360/150‑1 (MQ=255)
                                                                                                                                        ttggAATCGGCGCGGGGCTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTGACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg               <  2:456809/147‑1 (MQ=255)
                                                                                                                                         actaaTCGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg              <  1:3924510/147‑1 (MQ=255)
                                                                                                                                           tacggggCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc            <  1:3271739/145‑1 (MQ=255)
                                                                                                                                           gAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc            <  2:1447237/150‑1 (MQ=255)
                                                                                                                                           gAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGTTACATTAACCACATGCGATCAAATCCGGCTGGCAGCATAATGATACCATGGGCGCCCGGCTGCTGCTGggc            >  2:710990/1‑150 (MQ=255)
                                                                                                                                            aagcGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg           <  2:3914610/147‑1 (MQ=255)
                                                                                                                                            aaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg           <  2:2318050/150‑1 (MQ=255)
                                                                                                                                            aaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg           <  2:308038/150‑1 (MQ=255)
                                                                                                                                             aTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTTCTGTTCGGTACCTTGTCGATGTTACTTTAACCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg          >  1:4456246/1‑150 (MQ=255)
                                                                                                                                             aTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg          >  2:1527111/1‑150 (MQ=255)
                                                                                                                                               cggcgcggcgCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa        <  2:3677890/150‑1 (MQ=255)
                                                                                                                                                   gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  1:6559317/1‑150 (MQ=255)
                                                                                                                                                   gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  2:3404605/1‑150 (MQ=255)
                                                                                                                                                   gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt    >  2:3404468/1‑150 (MQ=255)
                                                                                                                                                     tGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTCCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc  >  2:5418815/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4163310‑4163608

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