Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 2,251,239 G→A P28P (CCC→CCT LY180_10740 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065842,251,2390GA97.9% 98.9 / ‑4.9 47P28P (CCC→CCTLY180_10740hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/22);  minor base T (1/0);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTT  >  CP006584/2250951‑2251471
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
cGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCCCGGTCAGCATCGGTGCCGTCCTGCGCGGCAGCACCTTAATCAGTGCATCGACCTCGGAGTTTGCGACCCCAGCCAGACTGACACCGCCCAGCGTCCCGGCCTTAGGTTtcctc                                                                                                                                                                                                                                  >  2:554471/1‑298 (MQ=255)
         tCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGAcc                                                                                                                                                                                                                         <  1:325883/298‑1 (MQ=21)
           aTGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTg                                                                                                                                                                                                                       >  1:473132/1‑298 (MQ=21)
            tGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGt                                                                                                                                                                                                                      <  1:478590/298‑1 (MQ=21)
                   cggtttcggtGGCTGCCAGGGAAATATCACGGCCACCGCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATaa                                                                                                                                                                                                               <  2:280155/289‑1 (MQ=255)
                        cTGATGGCCGCCAATGAAAAATCACTGCCCCATCCGCCCATAGGTCATCCGACGAGAGAATTTTCGGGAACCGCACCGCACCCGGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCGGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGc                                                                                                                                                                                                          <  1:353910/298‑1 (MQ=11)
                         tGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATGAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCt                                                                                                                                                                                                         >  1:166708/1‑298 (MQ=21)
                           aTGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCCCTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTca                                                                                                                                                                                                       <  1:550587/298‑1 (MQ=21)
                           aTGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTca                                                                                                                                                                                                       >  1:140714/1‑298 (MQ=21)
                            tGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCTGCAGCACCTTAATCAGTGCATCGACTTCGGATTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTTGGTTTCCTCAGAGTGACCTGGTCGCTACCCCGCCcac                                                                                                                                                                                                      >  2:348909/1‑298 (MQ=11)
                                 gcCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGc                                                                                                                                                                                                 <  1:101726/298‑1 (MQ=21)
                                    aGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTtcat                                                                                                                                                                                                                                   <  2:602009/261‑1 (MQ=21)
                                    aGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTtcat                                                                                                                                                                                                                                   >  1:602009/1‑261 (MQ=21)
                                         aaTATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCCCCACAAGTTCAGGCAGTCCCAGTGCGGCGACCTCCTGCTCGGTCCGCCTCGGTACCGTCCTGCGCGGCAGCCCCTTAATCGGTGCATCGACTTCGGCGTTTCGCACCCCCGCCAAACTCACTCCGCGCACCGGCCCGGCCTTAGTTGTCCGCCCCGTCACCTGCTCCAACACCCGCTcaccaccttctacccg                                                                                                                                                                                         >  2:242168/1‑288 (MQ=11)
                                           taccACTGCCACACCCGCCTACAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACCAAAAACCGGCCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGAt                                                                                                                                                                                       <  2:262139/295‑1 (MQ=255)
                                               aCTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTc                                                                                                                                                                                   >  1:77780/1‑298 (MQ=21)
                                                acgCCACAGTCCGCCCGCAGGTGCTCGACGGATGGCTGTGTCGGGCAGGGCCCTGCCCCGGGTGCGGCGACAAAAACCCGGCCACCCTACCGGCCAGCGCCACACGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCAGCGGTGCCGGCATGCGCGGCAGCACCTTAATCAGTGCATCGGCTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTc                                                                                                                                                                                   <  2:39980/296‑1 (MQ=11)
                                                  gCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCGCCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTAACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAg                                                                                                                                                                                >  1:788516/1‑298 (MQ=21)
                                                  gCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAg                                                                                                                                                                                <  1:237687/298‑1 (MQ=21)
                                                               aTCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCCCAGCGGCCCGGCATCAGGTTTCATCCTCGTCACCTGTTCGATCACCTGCCTCCCACACCTTACCGCTTTGTCCAGGGTGATctcctt                                                                                                                                                                    >  2:617518/1‑293 (MQ=11)
                                                                         cGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCCGTGCCGTCATGCGCGGCACCACCTTAATCAGTGCATCGACTTCGGCGTTTGCGACCGCAGCCAGACTGACACCGCCCAGCGTCCCGGCCTTAGGTTTCATCAGCGTCACCTGTTCGATAACCTGCTCACCACTCTTGCCCGAATTGTCCCGGGAAATCACATTTTCTTTGtta                                                                                                                                                          >  2:755096/1‑297 (MQ=11)
                                                                            ccGGGAGGTGTTGCGGGACCGTCAGTAAACCGAGGTCGGGGGCGAAAACCCGACCGCCTTGCCGGCCGGCGGCACGAGGTCAGGCCGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCAGCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCAGCGACGTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGGCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATg                                                                                                                                                       <  2:143763/298‑1 (MQ=11)
                                                                              gagagaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGt                                                                                                                                                     >  1:691076/1‑298 (MQ=21)
                                                                               agagaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGtt                                                                                                                                                    <  1:8561/298‑1 (MQ=21)
                                                                                gagacgtgtCGGGAGCGTCGCGGCCACGAGGTGGGGGACAACAGGCCGGCCACGTTCCCGGCCGGCGCCACAGGGTCAGGGCGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCGTCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACGTCGGAGTTTGCGACCGCAGCCAGACGGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGttt                                                                                                                                                   <  2:84902/290‑1 (MQ=11)
                                                                                 agaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGtttt                                                                                                                                                  >  1:397387/1‑298 (MQ=18)
                                                                                  gaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTtc                                                                                                                                                 <  1:554471/298‑1 (MQ=18)
                                                                                  gaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTtc                                                                                                                                                 >  1:158932/1‑298 (MQ=18)
                                                                                               cGTCACCGACCCGAGTTCGGCGACAAAAACCCGCCCCCCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATc                                                                                                                                   <  2:77780/298‑1 (MQ=255)
                                                                                                             gtccGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGCCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTCCATCCACCTCAGAGTTTTCGACCGCAACCCAACTGACACCGCCCAGCGCCCCCGCCTCAGGTTCACTCCGCCGTACCCGCTCCATAACCTGCTCACCCCGCCTGCACGGATACTGCCGGGTAAGCCCTTTTTCTTTTTTCATGATCTTCTCATCCTCGCTTCGGGGTTAACCgcta                                                                                                                     >  2:669204/4‑295 (MQ=2)
                                                                                                              ttCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCa                                                                                                                    <  1:702058/298‑1 (MQ=18)
                                                                                                              ttCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGGTTTCTCACTTCTGAATCGGGGTTAACCGGTCa                                                                                                                    >  1:144566/1‑298 (MQ=18)
                                                                                                                    gACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGAGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGGTTTCTCACTTCTGCATCGGGGTTAACCGGTCAGCCAgg                                                                                                              >  1:206804/1‑298 (MQ=18)
                                                                                                                         aaaaCCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGAc                                                                                                         >  1:420857/1‑298 (MQ=18)
                                                                                                                                      ttACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGCCCCGGCATTAGGTTTCATCAGCGTGCCCGGTTCGCTAACCTCCTCACCACGCTTGACCGGATCGTCCAGCGTAATCACATTTCCTTTGTTCATGGTTCTCTCCCCTCTGAACCTGGGTTACCCGGTCAGCCACGCTGACCGGATGAAAATcc                                                                                            >  2:259675/1‑297 (MQ=11)
                                                                                                                                                            cAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGgtgt                                                                      <  1:275238/298‑1 (MQ=18)
                                                                                                                                                              ggCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTg                                                                    >  1:431845/1‑298 (MQ=18)
                                                                                                                                                               gCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGc                                                                   <  1:242168/298‑1 (MQ=18)
                                                                                                                                                                  gTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGTAAATCACAGGCCGATATTGCGGCGGTGTTGCTcc                                                                >  1:411616/1‑298 (MQ=18)
                                                                                                                                                                    tCCAGGGCGGCGACCGCCTGCTCGGTCAGACTCGGGGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAg                                                              <  2:72035/298‑1 (MQ=12)
                                                                                                                                                                        gTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGCCCGGATGAAAATCACAGGCCGATATTGCGGCCGTGTTGCTCCACCCgg                                                          >  2:813625/1‑298 (MQ=17)
                                                                                                                                                                                        cTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACCGGACGCTATTGCGGCGGTGGTGGCCCAGCCGGTCAGCCCCGTTccact                                          >  2:274428/1‑294 (MQ=11)
                                                                                                                                                                                         tCGGGCAGCATGGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGCCTTCGGAGGTTGCGACCGCAGCCAGACTGGCACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACGTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCtt                                         <  2:648256/298‑1 (MQ=9)
                                                                                                                                                                                            gTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGCGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGCTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCCCAGGCCGATATTGCGGCGGTGTTGCCCCCGCCCGTCCACTCCGTTCACCTtctc                                      >  2:60901/1‑298 (MQ=255)
                                                                                                                                                                                              cAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCaa                                    <  1:396948/298‑1 (MQ=17)
                                                                                                                                                                                                                         caCCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCa         <  1:348909/298‑1 (MQ=14)
                                                                                                                                                                                                                                aTCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCtt  >  1:339247/1‑298 (MQ=14)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
CGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTT  >  CP006584/2250951‑2251471

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: