Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 2,265,396 A→G intergenic (‑149/‑25) LY180_10855 ← / → LY180_10860 precorrin‑8W decarboxylase/terminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065842,265,3960AG97.7% 73.4 / NA 45intergenic (‑149/‑25)LY180_10855/LY180_10860precorrin‑8W decarboxylase/terminase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (28/15);  minor base T (0/1);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCATTTCCTGAATATCCTGCGCACTGATGACACGCCCGTCACAGGTGTCACCCTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGGCACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCATGGATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAACCGCAGAAAACCGGCGGTGACTTCAAGGAAATTGACCTGCTCGGACGCC  >  CP006584/2265151‑2265688
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    
cgcccgtttCCTGAATGCGCGCCGCACTGATGACACGCCCGTGACAGGTGTCACCCTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCGCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCttt                                                                                                                                                                                                                                                  <  2:172748/292‑1 (MQ=1)
                  cccgcgcACATGATGACCCGCCCGGCACAGGGGTCGCCCTCAACGCCGATACGCAAGAATTGTGGGACTTTTTTTGCCAGTGTCAGGGGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGCCGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGGCACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCATGGATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGcgtcgt                                                                                                                                                                                                                                 <  2:445136/295‑1 (MQ=1)
                                                       ccTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGc                                                                                                                                                                                            >  1:443800/1‑298 (MQ=14)
                                                         tCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCg                                                                                                                                                                                          >  1:228376/1‑298 (MQ=14)
                                                             cGCCGTGACGAAAGGATGTTGGTGCTTTTTTTGCCGTGGGCGGGAGTCCTGAATAGTGAGGAGAGGGGTCATATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTGGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCAGCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGa                                                                                                                                                                                      <  2:895185/298‑1 (MQ=11)
                                                               actaTAAGAAAGAATGTTGAGACTTTTTTGGGCTATGTCCGGAATGCTGAAGAGTGATGAGAGGAGTAACCGATCGGCAGCATTTTCCCGACGATGCGCAGCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGAtt                                                                                                                                                                                    <  2:479560/295‑1 (MQ=11)
                                                                cGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGGACTGGCAGGGATTTTCCGTGCCGCAGATTg                                                                                                                                                                                   >  1:303618/1‑298 (MQ=14)
                                                                   tACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTgccg                                                                                                                                                                                >  1:239538/1‑298 (MQ=14)
                                                                    aCGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTgccgc                                                                                                                                                                               >  1:644118/1‑298 (MQ=14)
                                                                            atttgtacaCTTTTTTTGCCATTGTGAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGGAGCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGc                                                                                                                                                                       <  2:552405/290‑1 (MQ=9)
                                                                            aTTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGc                                                                                                                                                                       >  1:410111/1‑298 (MQ=14)
                                                                             ttttGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGCTGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGCCTCCGTAGCCTTGCAGCGTATTCATCACGGCGCGGTATTCATGACCATCCCCACAGACACCCCTCTTTTACAACACCCGCCGCGTCCGGCGGCGACGCTGTACTGGCAGGGACTTTCCGTTCCGCCGGTTGccaccagctcgca                                                                                                                                                                      >  2:305904/1‑291 (MQ=11)
                                                                                 gagaCTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATg                                                                                                                                                                  >  1:637556/1‑298 (MQ=14)
                                                                                    aCTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGaaa                                                                                                                                                               >  1:940448/1‑298 (MQ=14)
                                                                                        tttttGCCAGTGTCGGGGGTCCTGAAGGGTAGTTAGAGGAGTACCAGGTCGGGATCAGTTTCCCGACGATGCGCATCCTCCGCCATCGGTCCCGGGTGGCTTAGCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGcc                                                                                                                                                           <  2:439752/298‑1 (MQ=11)
                                                                                               cATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGAGTAAACGCCCGACGGt                                                                                                                                                    >  1:82007/1‑298 (MQ=9)
                                                                                                 ttGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGc                                                                                                                                                  >  1:882989/1‑298 (MQ=14)
                                                                                                  tGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCa                                                                                                                                                 <  1:538738/298‑1 (MQ=14)
                                                                                                  tGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCa                                                                                                                                                 <  1:546197/298‑1 (MQ=14)
                                                                                                          gTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGa                                                                                                                                         >  1:958879/1‑298 (MQ=14)
                                                                                                                          agagGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGCCCCGCGTCGTCAGGCGGCGCCGCCGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGCTGCAGATGAAACGCCCGACGGGGCCGTGCAGGAACCAACGCAGCGGCTg                                                                                                                         >  2:1017001/1‑298 (MQ=11)
                                                                                                                           gaggagTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACCGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACCGGCAGGGACTTTCCGTGCCGCAGCTTCCCCCCCCGCTGCAGATGAAACGCCCGGCGGTCCACAGCGCGCAACAGCGCGACGCCtcg                                                                                                                        >  2:673928/1‑296 (MQ=11)
                                                                                                                              gagTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTCCACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGGTGAAACCCCCGCCGGTGCCGAGCAGGAGACACACCGACGGCCGGGAc                                                                                                                     >  2:270846/1‑298 (MQ=11)
                                                                                                                                   acctGTCGGCCTCAGTGTCCCGGCGAGGCGCATCGCCGGCGAGTAGTGCCGGTTGGCTTAGCACTGACACTACGGGAGCTGAGCGAAGTGCGGGGCGCGGCGCGGTAGCCGAGCCGTGTAGTCATCACGGGGGGGTATTCATGACCATCACCACAGACACCACTCGTTTACACGACCCGCGTCGTCAGGCGGCGCGGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGTCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGt                                                                                                                <  2:155552/295‑1 (MQ=11)
                                                                                                                                           gCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCa                                                                                                        <  1:137273/298‑1 (MQ=14)
                                                                                                                                             aTCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATc                                                                                                      <  1:169442/298‑1 (MQ=14)
                                                                                                                                               cAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGc                                                                                                   >  1:544224/1‑298 (MQ=14)
                                                                                                                                                           gaTGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCGGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAAt                                                                                        <  1:503868/298‑1 (MQ=14)
                                                                                                                                                                  aTCCTCCGCCATCAGTCCCGGATGGCTTATCACTGAAACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCCTGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCCACAGCCGTGTCGAAAAGAGTCTg                                                                                 >  2:659928/1‑298 (MQ=9)
                                                                                                                                                                     ctccGCCATCCTACCCGGTTGGCTACTCACTGCAACACCCCCACCTCTCCGAATCGCTGGGCCCGACTCATTAGCCTTGCCTTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGGGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGaa                                                                              <  2:161188/298‑1 (MQ=11)
                                                                                                                                                                      tccGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAg                                                                             >  1:44985/1‑298 (MQ=14)
                                                                                                                                                                                    cGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGGAGCGCGGCCTACCCAg                                                               >  1:854043/1‑298 (MQ=9)
                                                                                                                                                                                     ggATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACAGTCAGGGACTTTCCGTGCCGCAGATTGCCGCCCTGTTGCAGACGACACGCCCGCCGGTGCCAAGCTGGGAACAGCCCGACGGCCGGGACCGCGTGGCCCCCACCCGCCGTGTCCAAAAGCGACAGGAAgccgcgcggcacccgc                                                              >  2:536370/1‑284 (MQ=11)
                                                                                                                                                                                      gATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCt                                                             <  1:433976/298‑1 (MQ=14)
                                                                                                                                                                                                cACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACACGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGGGTCGAAATGAGTCTGGAAGCGCGGCGGACCCAGCTAATTATCaaa                                                   >  1:127235/1‑298 (MQ=9)
                                                                                                                                                                                                 aCTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTCGCAGGGATTTTCCGTGCCGCAGATTGCCGCCCTGTCGCAGATGAAACGCCCGACGGTGCAGCGCTGGAAACAGCGCGCCGGCTGGGACAGCGATGCCCCCCTCAGCCCTGCCGAACTGAGTCTGGACGAGCGGCGCACCAACCCCATCATCAAAc                                                  >  2:961817/1‑298 (MQ=11)
                                                                                                                                                                                                  cTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAAcc                                                 <  1:12193/298‑1 (MQ=14)
                                                                                                                                                                                                         aacaGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAACCGCAGaaa                                          >  1:494958/1‑298 (MQ=14)
                                                                                                                                                                                                                            gcgGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGTTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTTCCGAAATGAGTCCGGAACCGCGGCTGACCCACCGCCCCCTCAACCCGCCCAAACCCGACGGTGAGTTAAAgg                       >  2:675593/1‑298 (MQ=11)
                                                                                                                                                                                                                               gggCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAACCGCAGAAAACCGGCGGTGACTTCAAGGaaa                    <  1:324089/298‑1 (MQ=14)
                                                                                                                                                                                                                                       cTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAACCGCAGAAAACCGGCGGTGACTTCAAGGAAATTGACCTg            >  1:366760/1‑298 (MQ=14)
                                                                                                                                                                                                                                            tAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATCTTGCAGATGACACGCCCTAAGGTGCAGAGCTGGACACAGCCCTACGGCCAGGACAGCGTTCCCCCCACAAGCCATGTCTAAATAAGTCTGGAAACCCGTCGGCCCAAGATCCTCCACCCACCGCAACAAACCGGCGACGACTTCAAAGAAATTCAccgccccgc       >  2:910475/1‑291 (MQ=11)
                                                                                                                                                                                                                                                 ttGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCaa                                                    >  2:169522/1‑248 (MQ=14)
                                                                                                                                                                                                                                                 ttGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCaa                                                    <  1:169522/248‑1 (MQ=14)
                                                                                                                                                                                                                                                 ttGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAACCGCAGAAAACCGGCGGTGACTTCAAGGAAATTGACCTGCTCGGACGcc  >  1:373613/1‑298 (MQ=14)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    
CGGCCATTTCCTGAATATCCTGCGCACTGATGACACGCCCGTCACAGGTGTCACCCTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGGCACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCATGGATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCAAACCGCAGAAAACCGGCGGTGACTTCAAGGAAATTGACCTGCTCGGACGCC  >  CP006584/2265151‑2265688

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