Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 612,755 G→C 8.6% P477R (CCG→CGG)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947612,7550GC 8.6% 254.9 / 11.6 81P477R (CCG→CGG) dxs1‑deoxyxylulose‑5‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (32/42);  new base C (4/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCG  >  NC_002947/612608‑612901
                                                                                                                                                    |                                                                                                                                                  
aTGGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTTTTTTTTTTTTCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccggg                                                                                                                                                  >  2:811382/1‑150 (MQ=255)
aTGGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccggg                                                                                                                                                  >  2:811574/1‑150 (MQ=255)
aTGGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccggg                                                                                                                                                  >  2:370350/1‑150 (MQ=255)
     cTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCttt                                                                                                                                             <  2:1236158/150‑1 (MQ=255)
     cTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                             <  1:926118/150‑1 (MQ=255)
      tCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTg                                                                                                                                            >  1:527730/1‑150 (MQ=255)
          ccAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAACGGCCCGGCCGCTGTAcg                                                                                                                                        <  1:1304060/150‑2 (MQ=255)
           cAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAg                                                                                                                                       <  1:753741/150‑1 (MQ=255)
           cAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAg                                                                                                                                       <  2:1337342/150‑1 (MQ=255)
                gTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGc                                                                                                                                  <  1:370350/150‑1 (MQ=255)
                    gCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGt                                                                                                                              >  1:696616/1‑150 (MQ=255)
                     ccGAACACCCGAGTGGCCATTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGGGGGTTTCCGCTGATCGGCGCATTCGGGCACGTCCCCCGCATCTACTTTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTg                                                                                                                             <  1:1391128/150‑1 (MQ=255)
                             cagcagGGCCACTTTGTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCAt                                                                                                                     >  1:343963/1‑150 (MQ=255)
                               gcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCt                                                                                                                   >  2:116158/1‑150 (MQ=255)
                                 ggggCCCCTTTTTCGCCCTGGCGGCGCACCACTCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTg                                                                                                                 <  1:262651/149‑1 (MQ=255)
                                 agGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTg                                                                                                                 <  2:1289759/150‑1 (MQ=255)
                                     ccACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCa                                                                                                             <  1:215711/150‑1 (MQ=255)
                                               gCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTtctc                                                                                                   <  1:116158/150‑1 (MQ=255)
                                                cccTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcg                                                                                                  >  1:1382893/1‑150 (MQ=255)
                                                 ccTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcgt                                                                                                 <  2:723036/150‑1 (MQ=255)
                                                      cggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTg                                                                                            <  1:284367/150‑1 (MQ=255)
                                                       ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAACGGCCCGGCCGCTGTAc                                                                                                                                         <  2:1353436/104‑1 (MQ=255)
                                                       ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAACGGCCCGGCCGCTGTAc                                                                                                                                         >  1:1353436/1‑104 (MQ=255)
                                                         cggACCACGCCCTTGCCGATCTCCAGTGGTTGGAGGTCTCCGCTGATCGGCGCAGTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                         <  1:23739/150‑1 (MQ=255)
                                                         cggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGAGCAGTTCGTTCTCGTCGCTgggg                                                                                         <  1:226369/150‑1 (MQ=255)
                                                         cggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTCATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                         <  1:587063/150‑1 (MQ=255)
                                                         cggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCAGCTGAGCGGCGCATTCGGGCCGGTACCACGCGAGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                         <  1:1163582/150‑1 (MQ=255)
                                                              cACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCg                                                                                                                                >  1:605/1‑106 (MQ=255)
                                                              cACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCg                                                                                                                                <  2:605/106‑1 (MQ=255)
                                                              cACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtcat                                                                                    <  2:712063/150‑1 (MQ=255)
                                                              cACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtcat                                                                                    >  1:64682/1‑150 (MQ=255)
                                                                   ccTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGGGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               >  1:661649/1‑150 (MQ=255)
                                                                   ccTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               <  2:1382893/150‑1 (MQ=255)
                                                                     ttGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGc                                                                             <  1:331632/150‑1 (MQ=255)
                                                                      tGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCa                                                                            <  2:343963/150‑1 (MQ=255)
                                                                         cGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               <  2:539792/144‑1 (MQ=255)
                                                                         cGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               >  1:539792/1‑144 (MQ=255)
                                                                          gATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGcc                                                                        >  1:197536/1‑150 (MQ=255)
                                                                            tctcCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCg                                                                      <  2:1200214/150‑1 (MQ=255)
                                                                             ctcCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCgg                                                                     <  1:211646/150‑1 (MQ=255)
                                                                               ccAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGa                                                                   >  2:1323462/1‑150 (MQ=255)
                                                                               ccAGTGGTTGCAGGTCCCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGCCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGa                                                                   <  1:1023618/150‑1 (MQ=255)
                                                                                cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGc                                                                                                                       >  1:806075/1‑97 (MQ=255)
                                                                                cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGc                                                                                                                       <  2:806075/97‑1 (MQ=255)
                                                                                     gTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAgc                                                             >  2:81730/1‑150 (MQ=255)
                                                                                        gCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGc                                                                                                                       >  2:877430/1‑89 (MQ=255)
                                                                                        gCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGc                                                                                                                       <  1:877430/89‑1 (MQ=255)
                                                                                        gCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCa                                                          >  1:353702/1‑150 (MQ=255)
                                                                                          aGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtc                                                                                      >  1:170558/1‑120 (MQ=255)
                                                                                          aGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtc                                                                                      <  2:170558/120‑1 (MQ=255)
                                                                                            gTCTCCGCTGATCGGCGCATTCGGGCCGGTACCGCGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               <  2:311044/125‑1 (MQ=255)
                                                                                            gTCTCCGCTGATCGGCGCATTCGGGCCGGTACCGCGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               >  1:311044/1‑125 (MQ=255)
                                                                                            gTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTa                                                      >  2:259001/1‑150 (MQ=255)
                                                                                               tcCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGa                                                   <  2:353702/150‑1 (MQ=255)
                                                                                                  gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcgtc                                                                                                >  1:481776/1‑102 (MQ=255)
                                                                                                  gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcgtc                                                                                                <  2:481776/102‑1 (MQ=255)
                                                                                                    tGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGt                                              <  2:197536/150‑1 (MQ=255)
                                                                                                         ggCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAg                                         >  1:488102/1‑150 (MQ=255)
                                                                                                           cgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               <  1:1063748/110‑1 (MQ=255)
                                                                                                           cgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               >  2:1063748/1‑110 (MQ=255)
                                                                                                             cATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGc                                     <  2:527730/150‑1 (MQ=255)
                                                                                                               ttCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtca                                                                                     >  2:674405/1‑100 (MQ=255)
                                                                                                               ttCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtca                                                                                     <  1:674405/100‑1 (MQ=255)
                                                                                                                  gggCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTtctc                                                                                                   <  2:893908/83‑1 (MQ=255)
                                                                                                                  gggCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTtctc                                                                                                   >  1:893908/1‑83 (MQ=255)
                                                                                                                  gggCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCg                                >  2:385234/1‑150 (MQ=255)
                                                                                                                    gCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTg                              <  1:1323462/150‑1 (MQ=255)
                                                                                                                     ccGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTgg                             <  1:81730/150‑1 (MQ=255)
                                                                                                                     acGGTACCCGCGTGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATACCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTgg                             >  1:340687/2‑150 (MQ=255)
                                                                                                                          aCCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCgg                        <  1:259001/150‑1 (MQ=255)
                                                                                                                             aCGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGcc                     >  1:423151/1‑150 (MQ=255)
                                                                                                                               gcgtggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGctg                    >  2:1100829/5‑148 (MQ=255)
                                                                                                                               gcgtggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCg                    >  1:631066/5‑150 (MQ=255)
                                                                                                                                 gggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               >  1:575563/1‑88 (MQ=255)
                                                                                                                                 gggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               <  2:575563/88‑1 (MQ=255)
                                                                                                                                   ccacctgTACAACGGCCCGGCCGCTGTACGGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCt                >  2:534625/7‑150 (MQ=255)
                                                                                                                                   ccacctgTACAACGGCCCGGCCGCTGTACGGTGGCCTGTCGTTAACGTCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTc               >  2:984769/7‑150 (MQ=255)
                                                                                                                                    gTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTc               <  1:85709/150‑1 (MQ=255)
                                                                                                                                     tACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCg              >  1:1282415/1‑150 (MQ=255)
                                                                                                                                       ccGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGcc            >  1:997350/1‑150 (MQ=255)
                                                                                                                                        tgTACAACGGCCCGGCCGCTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGCTACGACCTGTCGTACCTGc                                     <  1:560058/123‑1 (MQ=255)
                                                                                                                                        tgTACAACGGCCCGGCCGCTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGCTACGACCTGTCGTACCTGc                                     >  2:560058/2‑124 (MQ=255)
                                                                                                                                          tACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGAc         >  1:926871/1‑150 (MQ=255)
                                                                                                                                            cAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCa       <  1:315972/150‑1 (MQ=255)
                                                                                                                                                 gccgggccgTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCg  <  1:385234/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATGGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCG  >  NC_002947/612608‑612901

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