Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,963,902 T→G 7.4% Y216D (TAT→GAT)  PP_5229 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,963,9020TG 7.4% 267.7 / 11.1 94Y216D (TAT→GAT) PP_5229hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base G (5/2);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CGTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTCGC  >  NC_002947/5963756‑5964050
                                                                                                                                                  |                                                                                                                                                    
cGTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCtata                                                                                                                                                   <  2:829357/150‑1 (MQ=255)
       gTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGa                                                                                                                                            <  1:1028161/150‑1 (MQ=255)
        ttCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAg                                                                                                                                           <  1:1122346/150‑1 (MQ=255)
         tCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCACATGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAgg                                                                                                                                          <  1:761784/150‑1 (MQ=255)
             tGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCt                                                                                                                                      >  2:486534/1‑150 (MQ=255)
                agaaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTgg                                                                                                                                   <  2:769081/150‑1 (MQ=255)
                  aaCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGcc                                                                                                                                 >  1:989710/1‑150 (MQ=255)
                   aCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTgg                                                                                                                                   <  1:684971/147‑1 (MQ=255)
                   aCAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTgg                                                                                                                                   >  2:684971/1‑147 (MQ=255)
                    cAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCcgc                                                                                                                               >  2:1950464/1‑150 (MQ=255)
                       cgcCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGtt                                                                                                                            <  2:220286/150‑1 (MQ=255)
                        gcCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtg                                                                                                                           <  1:252657/150‑1 (MQ=255)
                         ccGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                          <  2:534691/150‑1 (MQ=255)
                         ccGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                          <  1:765720/150‑1 (MQ=255)
                         ccGGGAAGTAGCCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgt                                                                                                                          <  1:384578/150‑1 (MQ=255)
                          cGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTtgtg                                                                                                                         <  2:989710/150‑1 (MQ=255)
                          cGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATAGCGCCGAGGTGCTTGGCTGCGTtgtg                                                                                                                         >  1:2307681/1‑150 (MQ=255)
                                gTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCAcgc                                                                                                                   >  1:305151/1‑150 (MQ=255)
                                 tAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAgg                                                                                                                                          <  1:779388/126‑1 (MQ=255)
                                 tAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCGATATAGCCGAgg                                                                                                                                          >  2:779388/1‑126 (MQ=255)
                                     cTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGACGTGCTTGGCCGCGTTGTGCCACGCGTCAc                                                                                                              <  1:1261597/150‑1 (MQ=255)
                                       ccAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAccc                                                                                                            >  1:1662962/1‑150 (MQ=255)
                                        cAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCa                                                                                                           >  1:2057701/1‑150 (MQ=255)
                                        cAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCa                                                                                                           >  1:1248685/1‑150 (MQ=255)
                                        cAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCa                                                                                                           <  1:2047551/150‑1 (MQ=255)
                                                tGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTgc                                                                                                   <  2:2057701/150‑1 (MQ=255)
                                                 ggCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTgcg                                                                                                  >  2:150111/1‑150 (MQ=255)
                                                  gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGt                                                                                                                 >  1:124454/1‑134 (MQ=255)
                                                  gCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGt                                                                                                                 <  2:124454/134‑1 (MQ=255)
                                                   cAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGcc                                                                                                >  1:2287962/1‑150 (MQ=255)
                                                       tctTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCtt                                                                                                                                     >  2:431106/1‑109 (MQ=255)
                                                       tctTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCtt                                                                                                                                     <  1:431106/109‑1 (MQ=255)
                                                       tctTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGt                                                                                            >  2:2189578/1‑150 (MQ=255)
                                                       tctGGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCTTCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGt                                                                                            >  2:1090650/1‑150 (MQ=255)
                                                       tattGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGAGCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGt                                                                                            >  2:834730/3‑150 (MQ=255)
                                                        ctTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTACATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTa                                                                                           >  2:1257459/1‑150 (MQ=255)
                                                         ttGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTAc                                                                                          >  2:117331/1‑150 (MQ=255)
                                                              tACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtga                                                                                     <  1:2005885/150‑1 (MQ=255)
                                                                ccAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTgc                                                                                                   <  1:1347535/134‑1 (MQ=255)
                                                                ccAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTgc                                                                                                   >  2:1347535/1‑134 (MQ=255)
                                                                 cAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCtgatga                                                                                  >  1:921102/1‑150 (MQ=255)
                                                                   ggggCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGaaa                                                                                >  1:1830349/1‑150 (MQ=255)
                                                                   ggggCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGaaa                                                                                <  2:693064/150‑1 (MQ=255)
                                                                   ggggCTCTATGGGGTTCCCGCCATCGGGAGACGGAACCCGCCACACTCCAACCCTATCATCGGTACCCGTTCCCTCCTCTACATCGCCGACGTGCTTCCGCTCGTCGTGCCGCCCGGCACCCACACCCTGCCCCCGGTACGCTGATGaaa                                                                                <  2:123840/150‑1 (MQ=255)
                                                                        tGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCa                                                                           <  2:1248685/150‑1 (MQ=255)
                                                                           aTGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCt                                                                        >  1:2384288/1‑150 (MQ=255)
                                                                             ggggTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTgg                                                                      >  2:74774/1‑150 (MQ=255)
                                                                              gggTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGa                                                                     >  1:449051/1‑150 (MQ=255)
                                                                              gggTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGa                                                                     >  1:9193/1‑150 (MQ=255)
                                                                              gggTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGa                                                                     >  1:1409465/1‑150 (MQ=255)
                                                                                gTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAGACGCCAGCTGGAgg                                                                   <  2:1846443/150‑1 (MQ=255)
                                                                                    tGGCCATCGGCAGCCGAGACCCGCAACACTACAAGAGCATCGTCGGGACCGTTTTCCTCGGCTAGATCGCCGAGGTGGTTGACCGCAATGTGCGAAACGCCATACAGACGCTGCGCCCGAGACGTTAATTAAACCCCAGCCTGATGCttg                                                               >  2:1012483/1‑149 (MQ=255)
                                                                                     ggCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTAt                                                              >  2:2409723/1‑150 (MQ=255)
                                                                                      gCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTAtt                                                             <  2:1653070/150‑1 (MQ=255)
                                                                                       ccATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAc                                                                                                              <  1:1037987/100‑1 (MQ=255)
                                                                                       ccATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAc                                                                                                              >  2:1037987/1‑100 (MQ=255)
                                                                                        cATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGc                                                                                                                                  >  1:633506/1‑79 (MQ=255)
                                                                                        cATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGc                                                                                                                                  <  2:633506/79‑1 (MQ=255)
                                                                                              cAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCacac                                                     <  2:461548/150‑1 (MQ=255)
                                                                                                   gcgACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGTTGAAACGCCAGCTGGAGGCTTATTGCGCACACCTgcg                                                >  1:2308395/1‑150 (MQ=255)
                                                                                                     gACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGATTATTGCGCAGACATGTGCa                                              >  1:340273/1‑150 (MQ=255)
                                                                                                       cccGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAAc                                            <  1:74774/150‑1 (MQ=255)
                                                                                                       cccGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAAc                                            >  2:1458665/1‑150 (MQ=255)
                                                                                                        ccGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACg                                           <  1:2409723/150‑1 (MQ=255)
                                                                                                          gCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAg                                         >  1:743597/1‑150 (MQ=255)
                                                                                                           cAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                        <  2:2384288/150‑1 (MQ=255)
                                                                                                            aaCACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCa                                                                           >  2:410010/1‑114 (MQ=255)
                                                                                                            aaCACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCa                                                                           <  1:410010/114‑1 (MQ=255)
                                                                                                             acacTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgcgc                                      <  1:117331/150‑1 (MQ=255)
                                                                                                              cacTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGcc                                     >  2:49294/1‑150 (MQ=255)
                                                                                                                cTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAg                                   <  1:834730/150‑1 (MQ=255)
                                                                                                                 tACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTTGAGGCGTATTGCGCACACCTGCGCAGCGAGCGCCAgg                                  >  2:421093/1‑150 (MQ=255)
                                                                                                                  aCAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgt                                 <  1:2189578/150‑1 (MQ=255)
                                                                                                                  aCAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTACATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGgt                                 <  1:1257459/150‑1 (MQ=255)
                                                                                                                      gagCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTAc                                                                                          >  2:2007918/1‑89 (MQ=255)
                                                                                                                      gagCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTAc                                                                                          <  1:2007918/89‑1 (MQ=255)
                                                                                                                      gagCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCg                             <  1:1159711/150‑1 (MQ=255)
                                                                                                                         cagcGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGACGACTGGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAg                          <  1:150111/150‑1 (MQ=255)
                                                                                                                                 gTACCCTTGTGCCCGGCGGTATCGCCGAGGGGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTc                  >  2:1241852/1‑150 (MQ=255)
                                                                                                                                    cccTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCa               <  1:1090650/150‑1 (MQ=255)
                                                                                                                                       ttttCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAggg            <  2:2287962/150‑1 (MQ=255)
                                                                                                                                       ttttCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAggg            >  2:2251228/1‑150 (MQ=255)
                                                                                                                                       ttttCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAggg            >  1:1759434/1‑150 (MQ=255)
                                                                                                                                        tttCCTCGGCTATATCGCCGAGGTGCTTGGCCACGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGc           <  1:486534/150‑1 (MQ=255)
                                                                                                                                           ccTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTAc        >  1:1631338/1‑150 (MQ=255)
                                                                                                                                            cTCGGCTATATCGCCGAGGTGCTTGGCACCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGACCCGACAGCTGGAGGCTACTAGCCCTGGACTGCGCAACGAGTTCCAGGTGTCGGAGCACACCTCGCAGGGCTAcc       <  1:421093/150‑1 (MQ=255)
                                                                                                                                            cTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAAc                                            <  2:562194/113‑1 (MQ=255)
                                                                                                                                            cTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAAc                                            >  1:562194/1‑113 (MQ=255)
                                                                                                                                            cTCGGCGATATAGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTAcc       >  1:1926941/1‑150 (MQ=255)
                                                                                                                                             tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                        <  1:1397963/116‑1 (MQ=255)
                                                                                                                                             tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAgc                                        >  2:1397963/1‑116 (MQ=255)
                                                                                                                                             tCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCg      >  2:1934249/1‑150 (MQ=255)
                                                                                                                                                 cTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAccc                                                                                                            >  1:2149122/1‑44 (MQ=255)
                                                                                                                                                 cTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCAccc                                                                                                            <  2:2149122/44‑1 (MQ=255)
                                                                                                                                                 cTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTcgc  <  2:1409465/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGTTCCTGTTCGGTGAAGAACAGCGCCGGGAAGTAGGCTCCAGCGCCGTGGCAGCTCTTGAGTACCAGGGGCTGCATGGGGTTCTGGCCATCGGCAGCCGCGACCCGCAACACTACAAGAGCAGCGTCGGTACCCTTTTCCTCGGCTATATCGCCGAGGTGCTTGGCCGCGTTGTGCCACGCGTCACCCAGACCCTGCGCCCGGTACGCTGATGAAACGCCAGCTGGAGGCTTATTGCGCACACCTGCGCAACGAGCGCCAGGTGTCGGAGCACACCTCGCAGGGCTACCGTCGC  >  NC_002947/5963756‑5964050

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