Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 395,884 G→T 27.7% R166S (CGT→AGT)  ddlA ← D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273395,8840GT27.7% 533.2 / 170.1 238R166S (CGT→AGT) ddlAD‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base G (98/74);  new base T (35/31);  total (133/105)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 244 total aligned reads displayed.

ATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCA  >  CP009273/395735‑396033
                                                                                                                                                     |                                                                                                                                                     
aTCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGCGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACg                                                                                                                                                       <  2:4201876/150‑1 (MQ=255)
   gaactcgaacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGAtt                                                                                                                                                    <  2:331946/150‑1 (MQ=255)
      ctcgaacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGgc                                                                                                                                                 <  2:3827219/150‑1 (MQ=255)
      ctcgaacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGgc                                                                                                                                                 >  1:213203/1‑150 (MQ=255)
       tcgaacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGgcg                                                                                                                                                >  2:2313835/1‑150 (MQ=255)
        cgaacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGgcgc                                                                                                                                               <  2:1537516/150‑1 (MQ=255)
          aacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGgcgcgc                                                                                                                                             >  2:192535/1‑150 (MQ=255)
          aacGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGgcgcgc                                                                                                                                             <  1:3570873/150‑1 (MQ=255)
            cGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGt                                                                                                                                           >  1:2599500/1‑150 (MQ=255)
            cGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGt                                                                                                                                           >  2:2499305/1‑150 (MQ=255)
             gCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTc                                                                                                                                          >  1:1209396/1‑150 (MQ=255)
              ccAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCa                                                                                                                                         <  1:1375049/150‑1 (MQ=255)
                aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAgg                                                                                                                                       >  1:3240606/1‑150 (MQ=255)
                aGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAgg                                                                                                                                       >  1:499091/1‑150 (MQ=255)
                 gATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAggg                                                                                                                                      >  1:1339271/1‑150 (MQ=255)
                 gATCGACGGCAATTGAGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAggg                                                                                                                                      <  1:1580173/150‑1 (MQ=255)
                  aTCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGt                                                                                                                                     >  2:4075497/1‑150 (MQ=255)
                      aCGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGtaata                                                                                                                                 >  1:1547138/1‑150 (MQ=255)
                      aCGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGtaata                                                                                                                                 >  2:2429142/1‑150 (MQ=255)
                       cGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGtaataa                                                                                                                                <  2:4496406/150‑1 (MQ=255)
                          cAATTGCGTACTGTTCTTCACTGGTTAATTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCACAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATg                                                                                                                             >  1:3290388/1‑150 (MQ=255)
                           aaTTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATgg                                                                                                                            <  1:2066865/150‑1 (MQ=255)
                           aaTTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATgg                                                                                                                            >  1:1147674/1‑150 (MQ=255)
                           aaTTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATgg                                                                                                                            <  1:4494289/150‑1 (MQ=255)
                            aTTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGgc                                                                                                                           >  1:4583277/1‑150 (MQ=255)
                            aTTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGgc                                                                                                                           >  1:3509705/1‑150 (MQ=255)
                            aTTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGgc                                                                                                                           <  2:235965/150‑1 (MQ=255)
                                 gTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAAt                                                                                                                      <  2:824672/150‑1 (MQ=255)
                                 gTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAAt                                                                                                                      >  2:2769759/1‑150 (MQ=255)
                                 gTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAAt                                                                                                                      >  1:561986/1‑150 (MQ=255)
                                     tGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTc                                                                                                                  >  1:3307112/1‑150 (MQ=255)
                                     tGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTc                                                                                                                  >  1:1484417/1‑150 (MQ=255)
                                       ttcttcACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCTGGGTAATAAATGGCGCAATGTTCAg                                                                                                                <  2:1590780/150‑1 (MQ=255)
                                          ttcACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGccc                                                                                                             <  2:3323140/150‑1 (MQ=255)
                                          ttcACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGccc                                                                                                             <  2:1731174/150‑1 (MQ=255)
                                          ttcACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGccc                                                                                                             >  1:36400/1‑150 (MQ=255)
                                           tcACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCg                                                                                                            <  1:747787/150‑1 (MQ=255)
                                           tcACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCg                                                                                                            <  1:2807834/150‑1 (MQ=255)
                                                ggTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATc                                                                                                       >  1:2894232/1‑150 (MQ=255)
                                                ggTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATc                                                                                                       <  1:1281209/150‑1 (MQ=255)
                                                   tACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATcgcg                                                                                                    <  1:3193903/150‑1 (MQ=255)
                                                         gCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACCAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAg                                                                                              <  1:2888374/150‑1 (MQ=255)
                                                           tGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGAc                                                                                            <  1:2700665/150‑1 (MQ=255)
                                                           tGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGAc                                                                                            <  1:550763/150‑1 (MQ=255)
                                                                ccAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTg                                                                                       >  1:3640512/1‑150 (MQ=255)
                                                                   aCAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCCCTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGtt                                                                                    <  2:1484630/150‑2 (MQ=255)
                                                                     agaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGAc                                                                                  >  1:3752863/1‑150 (MQ=255)
                                                                     agaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGAc                                                                                  <  2:135758/150‑1 (MQ=255)
                                                                     agaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTAc                                                                                  >  2:2843192/1‑150 (MQ=255)
                                                                      gaagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACa                                                                                 >  2:3900301/1‑150 (MQ=255)
                                                                       aagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACAt                                                                                >  2:3998637/1‑150 (MQ=255)
                                                                       aagaGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACAt                                                                                >  2:1669221/1‑150 (MQ=255)
                                                                        agagCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATc                                                                               >  2:2560075/1‑150 (MQ=255)
                                                                         gagCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCt                                                                              >  2:151036/1‑150 (MQ=255)
                                                                          agCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCtt                                                                             >  2:129120/1‑150 (MQ=255)
                                                                            cccTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTg                                                                           <  1:1197046/150‑1 (MQ=255)
                                                                              cTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTc                                                                         >  1:449374/1‑150 (MQ=255)
                                                                                gATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCa                                                                       >  2:2474545/1‑150 (MQ=255)
                                                                                      ccGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAgg                                                                 <  2:67795/150‑1 (MQ=255)
                                                                                       cGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGgc                                                                >  2:2175275/1‑150 (MQ=255)
                                                                                        ggTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGgct                                                               >  2:3657096/1‑150 (MQ=255)
                                                                                        ggTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGgct                                                               <  2:1340012/150‑1 (MQ=255)
                                                                                           tttACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCTGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGgctgct                                                            <  2:1010357/150‑1 (MQ=255)
                                                                                           tttACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGgctgct                                                            <  2:1010384/150‑1 (MQ=255)
                                                                                            ttACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGctgctg                                                           <  1:1868152/150‑1 (MQ=255)
                                                                                                    aGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGa                                                   <  2:3787142/150‑1 (MQ=255)
                                                                                                       ggTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACa                                                <  2:2611427/150‑1 (MQ=255)
                                                                                                       ggTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACa                                                >  1:806795/1‑150 (MQ=255)
                                                                                                       ggTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACa                                                <  1:3578480/150‑1 (MQ=255)
                                                                                                         tAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATc                                              <  1:1880579/150‑1 (MQ=255)
                                                                                                          aaCCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCa                                             >  2:2691321/1‑150 (MQ=255)
                                                                                                            ccccAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGa                                           >  2:4372684/1‑150 (MQ=255)
                                                                                                             cccAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGaa                                          <  2:759520/150‑1 (MQ=255)
                                                                                                             cccAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGaa                                          >  1:3775658/1‑150 (MQ=255)
                                                                                                               cAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAAcc                                        <  2:4507137/150‑1 (MQ=255)
                                                                                                                aGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCt                                       >  2:1843497/1‑150 (MQ=255)
                                                                                                                aGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCt                                       >  1:2155057/1‑150 (MQ=255)
                                                                                                                 gTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTa                                      <  1:1171220/150‑1 (MQ=255)
                                                                                                                   ttAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACa                                    >  2:3494152/1‑150 (MQ=255)
                                                                                                                          ccACTTCGGCAAAACTAATGTTGTGACTATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTa                             <  1:3937500/150‑1 (MQ=255)
                                                                                                                            aCTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTaaa                           >  1:2274665/1‑150 (MQ=255)
                                                                                                                               tCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTg                        <  1:629335/150‑1 (MQ=255)
                                                                                                                                   cAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGa                    <  2:36400/150‑1 (MQ=255)
                                                                                                                                    aaaaCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGAc                   <  1:2429142/150‑1 (MQ=255)
                                                                                                                                       aCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCg                <  1:2886287/150‑1 (MQ=255)
                                                                                                                                       aCTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCg                >  1:1783293/1‑150 (MQ=255)
                                                                                                                                         tAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCgca              <  2:3196117/150‑1 (MQ=255)
                                                                                                                                         tAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCgca              <  2:613724/150‑1 (MQ=255)
                                                                                                                                            tgttgtGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCgcagca           <  1:2287242/150‑1 (MQ=255)
                                                                                                                                               tgtgACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTc        <  2:1881847/150‑1 (MQ=255)
                                                                                                                                               tgtgACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTc        <  2:3103599/150‑1 (MQ=255)
                                                                                                                                                gtgACGATATGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTcc       >  2:778085/1‑150 (MQ=255)
                                                                                                                                                 tgACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTccc      >  1:1389196/1‑150 (MQ=255)
                                                                                                                                                  gACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCt     >  2:2491263/1‑150 (MQ=255)
                                                                                                                                                   aCGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTg    <  2:2984292/150‑1 (MQ=255)
                                                                                                                                                   aCGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTg    <  1:1794183/150‑1 (MQ=255)
                                                                                                                                                   aCGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTg    >  2:572989/1‑150 (MQ=255)
                                                                                                                                                   aCGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTg    <  1:1240730/150‑1 (MQ=255)
                                                                                                                                                    cGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGc   >  1:4054493/1‑150 (MQ=255)
                                                                                                                                                     gATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTTACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCa  >  1:3986546/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCAAAACTAATGTTGTGACGATTGGCGCGCGTCAGGGTAATAAATGGCGCAATGTTCAGCCCGGCATCGCGTAGCAGACGTTTGGTGACATCTTTGTCCATACAGGCTGCTGAAGCCAGAACATCAGAACCTACAAACGGTAAATTGGCGACCCGCAGCATTCCCTGCA  >  CP009273/395735‑396033

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