Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 2,251,586 T→C E104E (GAA→GAG LY180_10745 ← major tail tube protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065842,251,5860TC97.3% 60.9 / ‑4.0 37E104E (GAA→GAGLY180_10745major tail tube protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (17/19);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

TGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTT  >  CP006584/2251322‑2251876
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                   
tgtttttCTCACTTCGTCATGGGGGTTAACCGGTGACCCGGGGTGACCGGGTGAAAATCAACGGCCGTTCTTCGGGCGGGGTTGCTCAAGCCGGTCGACGCCGTTCACCTTCTCATTCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCgcgc                                                                                                                                                                                                                                                                    <  2:339247/295‑1 (MQ=1)
                        gttacccagttggCCGGGCTGACCGGCGGAAAAGCACAGGCCTGTTGTGCGGGGGTGTGTCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTGGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGGCCAGTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATc                                                                                                                                                                                                                                            <  2:617547/286‑1 (MQ=1)
                                                                     tttGCGGGGGTGTTCCTCCCGCCGTCTGACGCCGTTCACCTTGCCAATCATGTTGATGGTGGCGATTTCCACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAAcc                                                                                                                                                                                               <  2:54267/297‑1 (MQ=11)
                                                                      ttGGGGCGGGGTGCCTCCGGCCGGTCGACGCCGTTCACCTTCTCAATCAGGTGGATGGGGTCGATTTCGACCAGCTCCTGACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCg                                                                                                                                                                                              <  2:788516/298‑1 (MQ=11)
                                                                                     tCCAGGCGGTGGCCGCCTGTCACCTCTCCAAGCATGTTGCGGGTGTCGATGTCGACCAGCTCCCTGCCGTCCCTCGTCAGCCGGACATAGCTGCACGCCACGGAGATTTTCGACTGGTTGTCTTCTCTCTGTTTGCCCTCGCCGGTGCCGATTTCTTTCTGACGTCCACGCGTGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAACGCTCCCAGATAACCGAATCCGGGAAGCCAc                                                                                                                                                                               <  2:125179/298‑1 (MQ=11)
                                                                                               cGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCa                                                                                                                                                                    >  1:341060/1‑298 (MQ=14)
                                                                                                        tCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATCGGTGCAGACCACGGAGATTTTCGCCTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGCGTCCATTTCTTTCTGACGTCCACGCATGACCCCCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCCCACCCCCCCCCCCGGTTGCGCCGTCAACTTCCCAGATACCCGAACCAGGGCACCCACACGAGcaccgccccacttc                                                                                                                                                            >  2:111943/1‑285 (MQ=11)
                                                                                                        tCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGa                                                                                                                                                           <  1:274428/298‑1 (MQ=14)
                                                                                                              tctcAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAgcgc                                                                                                                                                     <  1:583573/298‑1 (MQ=14)
                                                                                                                       ttgggATGGTGTCGATTCCGACCGGCTCCTTTACCTGCAGCGTCAGGCGGAAATAGGTGCGGACCACGGGGGTTTTCGGCTCGGTGTCGTCTCCCTGTTTGCCCTCGGCGTGGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGGAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGtcatc                                                                                                                                            <  2:431845/294‑1 (MQ=11)
                                                                                                                             tGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGcc                                                                                                                                      >  1:733012/1‑298 (MQ=14)
                                                                                                                                 gtCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAgg                                                                                                                                  <  2:341060/298‑1 (MQ=14)
                                                                                                                                              agttcgttACCGGCCACCGCCAGCCGGAAATAGGTGCAGCCCACGGAGATTTTCGACTCGGTGGCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTgcgc                                                                                                                     <  2:708825/292‑1 (MQ=255)
                                                                                                                                                cTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGCTTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGGCGCGCTGGTCAGAGCCTGCAAAACGAAGCGGCACCTCCTACACCCCGTCTGCGGCGTGAAGCTCCCACATACCCGATTCCGGGACGCCACCGAGAGACGACTCCATTACCAGCTCATCGCCATCAAGGCCCAGGCCTACCgctcccctg                                                                                                                   >  2:428508/1‑290 (MQ=11)
                                                                                                                                                cTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGCTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGCTCACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACGCCGCATCTTCAGCAAGGCCGAGGtacacacgggggcgg                                                                                                                   >  2:214814/1‑284 (MQ=11)
                                                                                                                                                    ttACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGCTAACCGAATCCGGGAAGCCACCCAGCGACCCCTCCCTTGACAGCGCATCGTCATCCAGGCCGAGGTCTCCCGGTGCCCTGCCgc                                                                                                               >  2:310463/1‑297 (MQ=255)
                                                                                                                                                    ttACCGACCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCATGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACACCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACAGGGGCGCTGCCga                                                                                                               >  1:381074/1‑297 (MQ=9)
                                                                                                                                                      aCCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTc                                                                                                             >  1:411137/1‑298 (MQ=14)
                                                                                                                                                                gTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACAGGTGCGCTGCCGTTCATCCCCGCAc                                                                                                   >  1:123334/1‑298 (MQ=14)
                                                                                                                                                                 tCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCCCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGGTCATCCCCGCAcc                                                                                                  >  1:612707/1‑298 (MQ=14)
                                                                                                                                                                                 gCAGACCACGGAGATTTTCGACTCCGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGCCCACCTCGACGGCCCCCGTTTCGCCGCTATCGTCGCGCTGCTAAGAGCCTGCCAAACCAATCGGCACACCCTCCACCCCGGTTGAGGCGTACCACCCCCAGATAACCCCATCCGCGAACCAACCCCGCCCCCACCCCCTTGACACCCCACCGTCATCAAGCACGACGCCTACCGGCCCGACGCCGCTCATCCCCACGCCCCGCTAGTTTTCGacc                                                                                  >  2:59378/1‑296 (MQ=2)
                                                                                                                                                                                      ccACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACg                                                                             <  1:538862/298‑1 (MQ=14)
                                                                                                                                                                                       cACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACgg                                                                            <  1:437846/298‑1 (MQ=14)
                                                                                                                                                                                        aCGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACggg                                                                           <  1:669998/298‑1 (MQ=14)
                                                                                                                                                                                             gaTTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGTTCAGt                                                                      >  1:177457/1‑298 (MQ=14)
                                                                                                                                                                                              aTTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGtt                                                                     <  1:783961/298‑1 (MQ=14)
                                                                                                                                                                                              aTTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGtt                                                                     >  1:118238/1‑298 (MQ=14)
                                                                                                                                                                                                tttCGACTCGGTGGCTTCTCCCGGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGtttt                                                                   <  1:617456/298‑1 (MQ=14)
                                                                                                                                                                                                         ggTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTg                                                          <  1:259675/298‑1 (MQ=14)
                                                                                                                                                                                                              gtcttcCCTGGTTGGCCCGCCCGTTGTCAGTTTCTTTCTTAGGGCCCCCCAGGACCCCCTCGCCGGCCACCGTTTCGCCGGTCTCGTCGGGCTGGGAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCCGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGa                                                     <  2:804791/294‑1 (MQ=11)
                                                                                                                                                                                                                                  ccGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGCCGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGc                                 <  1:441314/298‑1 (MQ=14)
                                                                                                                                                                                                                                      tgtCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTa                             >  1:600882/1‑298 (MQ=14)
                                                                                                                                                                                                                                             ttctttctGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTc                      <  1:38928/298‑1 (MQ=14)
                                                                                                                                                                                                                                                      aCGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACg             >  1:119613/1‑298 (MQ=14)
                                                                                                                                                                                                                                                           cACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCAt        >  1:506013/1‑298 (MQ=14)
                                                                                                                                                                                                                                                                aTGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGACGCCACCGCGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATCGTTTTCGAGCTTCCGCGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTTATAGCTTGAGCGGTTCAGAAGGACGTTCATTAAc  >  2:144242/1‑297 (MQ=11)
                                                                                                                                                                                                                                                                 tGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACtt  <  1:414209/298‑1 (MQ=14)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                   
TGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTT  >  CP006584/2251322‑2251876

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