Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,701,978 C→T 6.1% D136D (GAC→GAT PP_3267 → Clp protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,701,9780CT 6.1% 244.5 / 4.1 81D136D (GAC→GATPP_3267Clp protease
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (3/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

TGGGCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGA  >  NC_002947/3701836‑3702121
                                                                                                                                              |                                                                                                                                                
tGGGCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTgcgc                                                                                                                                            <  7:19342/149‑1 (MQ=255)
         tCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTgcgc                                                                                                                                            <  4:214360/140‑1 (MQ=255)
         tCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTgcgc                                                                                                                                            >  3:214360/1‑140 (MQ=255)
           ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGTATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                 <  8:244399/149‑1 (MQ=255)
           ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                   >  6:134361/1‑147 (MQ=255)
           ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                   >  6:318040/1‑147 (MQ=255)
           ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                 <  5:237757/149‑1 (MQ=255)
            tctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                   >  3:244118/1‑146 (MQ=255)
              ttCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                   >  3:254079/1‑144 (MQ=255)
                ccTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                  >  4:318293/1‑143 (MQ=255)
                ccTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                  <  3:318293/143‑1 (MQ=255)
                  ttGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTg                                                                                                                          >  7:12122/1‑149 (MQ=255)
                    gCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                         <  4:254079/148‑1 (MQ=255)
                    gCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGc                                                                                                                        >  1:317884/1‑149 (MQ=255)
                     ccGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCg                                                                                                                       <  8:12122/149‑1 (MQ=255)
                     ccGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCg                                                                                                                       <  5:134361/149‑1 (MQ=255)
                            gCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCa                                                                                                                <  3:122626/149‑1 (MQ=255)
                                        ttGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCgg                                                                                                    >  8:233047/1‑149 (MQ=255)
                                              tttCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                         >  1:26232/1‑122 (MQ=255)
                                              tttCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                         <  2:26232/122‑1 (MQ=255)
                                                  aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTAt                                                                                                                                      >  4:279476/1‑105 (MQ=255)
                                                  aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGtttt                                                                                          >  2:263364/1‑149 (MQ=255)
                                                    cAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                         >  3:1804/1‑148 (MQ=255)
                                                     aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                       <  1:263364/149‑1 (MQ=255)
                                                     aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                       <  5:318040/149‑1 (MQ=255)
                                                            ttttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGt                                                                                                                                        <  2:98626/93‑1 (MQ=255)
                                                            ttttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGt                                                                                                                                        >  1:98626/1‑93 (MQ=255)
                                                              ttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACcgc                                                                               <  4:244118/148‑1 (MQ=255)
                                                                 cACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGtt                                                                            <  2:317884/148‑1 (MQ=255)
                                                                 cACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTa                                                                           >  8:92007/1‑149 (MQ=255)
                                                                    ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                         >  5:13881/1‑132 (MQ=255)
                                                                    ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                         <  6:13881/132‑1 (MQ=255)
                                                                    ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                          <  3:274891/147‑1 (MQ=255)
                                                                    ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                          >  4:274891/1‑147 (MQ=255)
                                                                            ttGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGcc                                                                 >  7:95080/1‑148 (MQ=255)
                                                                                   ttGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATcc                                                            >  3:246297/1‑146 (MQ=255)
                                                                                     gCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                         <  4:246297/146‑1 (MQ=255)
                                                                                     gCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAgg                                                       <  8:95080/148‑1 (MQ=255)
                                                                                         ggTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGAGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACt                                                  >  3:54971/1‑149 (MQ=255)
                                                                                          gTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTg                                                 >  8:80732/1‑149 (MQ=255)
                                                                                           ttCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                         <  6:77094/109‑1 (MQ=255)
                                                                                             cGGTGGTCCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGgc                                               <  4:54971/148‑1 (MQ=255)
                                                                                               gTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGcgcc                                            >  3:177407/1‑149 (MQ=255)
                                                                                               gTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGcgcc                                            >  4:88931/1‑149 (MQ=255)
                                                                                                  gACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCcgcg                                         <  3:88931/149‑1 (MQ=255)
                                                                                                      aTTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAg                                     >  8:57256/1‑149 (MQ=255)
                                                                                                          gcgGATGTCCGAGTATGCGATGAGACTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGgcgc                                 >  8:75118/1‑149 (MQ=255)
                                                                                                           cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                <  4:177407/149‑1 (MQ=255)
                                                                                                           cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                <  7:233047/149‑1 (MQ=255)
                                                                                                             gATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATCTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTgg                              <  4:1804/149‑1 (MQ=255)
                                                                                                                  ccGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                         >  7:245398/1‑149 (MQ=255)
                                                                                                                   cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                        <  7:57256/149‑1 (MQ=255)
                                                                                                                   cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                        >  2:270413/1‑149 (MQ=255)
                                                                                                                    gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                  >  3:242293/1‑107 (MQ=255)
                                                                                                                    gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                  <  4:242293/107‑1 (MQ=255)
                                                                                                                    gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                       <  7:75118/149‑1 (MQ=255)
                                                                                                                          gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                     <  4:192790/82‑1 (MQ=255)
                                                                                                                          gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                     >  3:192790/1‑82 (MQ=255)
                                                                                                                          gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACggg                      >  6:299709/1‑144 (MQ=255)
                                                                                                                          gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACggg                      <  5:299709/144‑1 (MQ=255)
                                                                                                                                 gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                         <  1:19058/134‑1 (MQ=255)
                                                                                                                                 gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                         >  2:19058/1‑134 (MQ=255)
                                                                                                                                 gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGa           >  2:16068/1‑148 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                       <  4:66970/135‑1 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                       >  3:66970/1‑135 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                 <  6:270855/141‑1 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                 <  4:314963/141‑1 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                 >  5:270855/1‑141 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                 >  3:314963/1‑141 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa         >  6:21869/1‑149 (MQ=255)
                                                                                                                                  ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa         >  4:282453/1‑149 (MQ=255)
                                                                                                                                    tGGATTACGACCTGCGCCTGTATGCGCAGATAGTGGCGGAGCGCACGAAAGGCGCGGCCGAGCTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGCATGATCCATGAAATc       >  2:142269/1‑149 (MQ=255)
                                                                                                                                      gATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcga     <  5:21869/149‑1 (MQ=255)
                                                                                                                                       aTTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac    <  8:245398/149‑1 (MQ=255)
                                                                                                                                        ccaCGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac    >  7:324329/3‑148 (MQ=255)
                                                                                                                                           cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                >  5:353215/1‑117 (MQ=255)
                                                                                                                                           cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                <  6:353215/117‑1 (MQ=255)
                                                                                                                                           cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                        >  7:276019/1‑125 (MQ=255)
                                                                                                                                           cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                        <  8:276019/125‑1 (MQ=255)
                                                                                                                                           cGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCTCAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga  <  1:270413/148‑1 (MQ=255)
                                                                                                                                           cGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac    >  2:16461/1‑145 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TGGGCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGA  >  NC_002947/3701836‑3702121

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: