Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,712,562 C→A 75.9% D31Y (GAT→TAT)  mliC ← inhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,712,5620CA75.9% 450.0 / 228.0 299D31Y (GAT→TAT) mliCinhibitor of c‑type lysozyme, membrane‑bound; putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base A (111/116);  total (145/154)
Fisher's exact test for biased strand distribution p-value = 8.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 299 total aligned reads displayed.

AGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTG  >  CP009273/1712415‑1712708
                                                                                                                                                   |                                                                                                                                                  
aGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCgg                                                                                                                                                  >  1:1919357/1‑150 (MQ=255)
aGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAgg                                                                                                                                                  <  2:4175373/150‑1 (MQ=255)
  aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGtt                                                                                                                                                >  2:1459810/1‑150 (MQ=255)
  aTTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGtt                                                                                                                                                <  1:5309993/150‑1 (MQ=255)
   ttCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGttt                                                                                                                                               <  2:2017707/150‑1 (MQ=255)
      cGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCa                                                                                                                                            >  1:2354720/1‑150 (MQ=255)
        tCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATg                                                                                                                                          >  2:3256381/1‑150 (MQ=255)
         cAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATgc                                                                                                                                         >  1:2617120/1‑150 (MQ=255)
             ttatcgcgcgCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGGATGCGCTc                                                                                                                                     <  1:3128042/146‑1 (MQ=255)
             gtAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTc                                                                                                                                     <  1:450621/150‑1 (MQ=255)
               aaCGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCaa                                                                                                                                   >  1:4777462/1‑150 (MQ=255)
               aaCGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCaa                                                                                                                                   >  1:1722329/1‑150 (MQ=255)
                      gcCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTg                                                                                                                            >  1:4125599/1‑150 (MQ=255)
                      gcCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTg                                                                                                                            >  2:3000407/1‑150 (MQ=255)
                      gcCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTg                                                                                                                            <  2:4140647/150‑1 (MQ=255)
                        cAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGAt                                                                                                                          <  2:1836725/150‑1 (MQ=255)
                           agGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGAGTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTCTTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCGGATTaa                                                                                                                       >  1:4810982/1‑150 (MQ=255)
                            ggCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAATGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                      >  2:1248938/1‑150 (MQ=255)
                            ggCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTaaa                                                                                                                      <  2:4198817/150‑1 (MQ=255)
                               tGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGgc                                                                                                                   >  1:2084402/1‑150 (MQ=255)
                                    tGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTAc                                                                                                              <  1:1557745/150‑1 (MQ=255)
                                       ccTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAATGCGCTACAGc                                                                                                           >  1:4640509/1‑150 (MQ=255)
                                       ccTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGc                                                                                                           >  1:3937400/1‑150 (MQ=255)
                                       ccTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGc                                                                                                           >  1:2230514/1‑150 (MQ=255)
                                         tGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCg                                                                                                         >  1:1395288/1‑150 (MQ=255)
                                          gTTTGAGATGCNGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCgg                                                                                                        >  1:2759479/1‑150 (MQ=255)
                                            ttGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgag                                                                                                      >  1:221033/1‑150 (MQ=255)
                                            ttGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgag                                                                                                      >  2:4301468/1‑150 (MQ=255)
                                             tGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgaga                                                                                                     >  1:3240814/1‑150 (MQ=255)
                                             tGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgaga                                                                                                     >  1:1139322/1‑150 (MQ=255)
                                             tGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATGAAAGGCGCTACAGCCGgaga                                                                                                     >  2:4463126/1‑150 (MQ=255)
                                              gagaTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGgagag                                                                                                    <  2:3377696/150‑1 (MQ=255)
                                                 aTGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAg                                                                                                 <  1:1661083/150‑1 (MQ=255)
                                                  tGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAgg                                                                                                <  1:2502831/150‑1 (MQ=255)
                                                  tGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAgg                                                                                                <  1:847541/150‑1 (MQ=255)
                                                     agtagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGAcc                                                                                             <  1:5064618/150‑1 (MQ=255)
                                                       tagtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCgg                                                                                           >  2:366506/1‑150 (MQ=255)
                                                        agtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGc                                                                                          <  1:3618097/150‑1 (MQ=255)
                                                         gtTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCa                                                                                         >  1:3009566/1‑150 (MQ=255)
                                                           tGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCaaa                                                                                       <  2:5321663/150‑1 (MQ=255)
                                                            gATTATCGTAAACAAAACTGACCTCCTGGCGCGTATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCGGGAGAGAAGGACCGGCaaaa                                                                                      <  2:948340/150‑1 (MQ=255)
                                                              ttATCGTAAACATTACTGACCTCCTGGCGCGGATTTTTCAGTCTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCTACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATg                                                                                    <  2:5160794/150‑1 (MQ=255)
                                                              ttATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATg                                                                                    <  2:4622451/150‑1 (MQ=255)
                                                              ttATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATg                                                                                    <  2:3192377/150‑1 (MQ=255)
                                                              ttATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATg                                                                                    <  2:4798748/150‑1 (MQ=255)
                                                                 tCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATa                                                                                 >  2:1548845/1‑150 (MQ=255)
                                                                  cGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATaa                                                                                <  1:594230/150‑1 (MQ=255)
                                                                   gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAATATGATAAt                                                                               >  1:1296067/1‑150 (MQ=255)
                                                                   gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAt                                                                               >  1:5136161/1‑150 (MQ=255)
                                                                   gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAt                                                                               <  2:2483066/150‑1 (MQ=255)
                                                                   gTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAt                                                                               <  1:351775/150‑1 (MQ=255)
                                                                    tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAtt                                                                              >  1:4607015/1‑150 (MQ=255)
                                                                    tAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAtt                                                                              <  1:3816803/150‑1 (MQ=255)
                                                                     aaaCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTa                                                                             >  1:3077504/1‑150 (MQ=255)
                                                                      aaCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTaa                                                                            >  1:4289229/1‑150 (MQ=255)
                                                                       aCAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAAc                                                                           <  2:2401632/150‑1 (MQ=255)
                                                                        cAAAACTGACCTCCTGGCGCGGATTATTCAGTTTTACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCAAGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACa                                                                          >  1:1991157/1‑150 (MQ=255)
                                                                          aaaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGt                                                                        <  1:4098692/150‑1 (MQ=255)
                                                                           aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtt                                                                       <  2:4711791/150‑1 (MQ=255)
                                                                           aaCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtt                                                                       <  1:4886762/150‑1 (MQ=255)
                                                                            aCTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACAATGTTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTTAAGGCGCTACAGCCGGAGAGAAGGACCGTCAAAATGATAATTTACAGttt                                                                      <  2:2865702/150‑1 (MQ=255)
                                                                               gACCTCCTGGCGCTGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                   <  2:1004901/150‑1 (MQ=255)
                                                                               gACCTCCTGGCGCTGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGtttttt                                                                   <  2:3913408/150‑1 (MQ=255)
                                                                                aCCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTc                                                                  >  2:1656483/1‑150 (MQ=255)
                                                                                  ctcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCAt                                                                <  2:2647971/150‑1 (MQ=255)
                                                                                   tcctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATa                                                               <  2:2196368/150‑1 (MQ=255)
                                                                                    cctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAg                                                              <  2:4230620/150‑1 (MQ=255)
                                                                                     ctGGCGCGTATTATTCAGTTTGTCCGTCAACGGTTTTTCATCACTCTGGTATTCGAGCGTATAGGTTTGCATGCGCTCAAACTGCTGATTAAAGGAGCTACAGCCGGAGAGAAGTACCGGCAAAATGATAATTAACAGTTTTTTCATAGt                                                             <  1:4931947/150‑1 (MQ=255)
                                                                                     ctGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGt                                                             <  2:3700473/150‑1 (MQ=255)
                                                                                         cgcgGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCata                                                         <  1:4361507/150‑1 (MQ=255)
                                                                                            ggATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATcc                                                      <  1:892628/150‑1 (MQ=255)
                                                                                             gATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATccc                                                     <  1:3406416/150‑1 (MQ=255)
                                                                                             gATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATccc                                                     <  1:1208686/150‑1 (MQ=255)
                                                                                               ttattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGa                                                   <  1:4010484/150‑1 (MQ=255)
                                                                                                tattCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGaa                                                  <  2:3078861/150‑1 (MQ=255)
                                                                                                 attCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAg                                                 >  2:1933620/1‑150 (MQ=255)
                                                                                                     aGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACtt                                             >  1:2117116/1‑150 (MQ=255)
                                                                                                       tttGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTc                                           >  2:1521212/1‑150 (MQ=255)
                                                                                                        ttGACCGTCAACGGTTTTTCATCACACTGTTATTCCAGCGTATAGGTTTTCATGCGCTCAACCAGCTGATTAAATGCGCTACAGCCGGAGAGAAGGACAGTCAAAATGATAATTAACAGTTTTTTCATAGTAATATCACGAAGACTTTca                                          >  1:229628/1‑149 (MQ=255)
                                                                                                         tGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCt                                         >  2:3478777/1‑150 (MQ=255)
                                                                                                            ccGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGt                                      <  2:4782487/150‑1 (MQ=255)
                                                                                                                cAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGAGAATTAACAGTTTTTTCATCGTCATATCCCGAAGACTTTCCTGGTCTgg                                  >  2:3403990/1‑150 (MQ=255)
                                                                                                                   cGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAgg                               <  1:5236318/150‑1 (MQ=255)
                                                                                                                     gTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGc                             >  1:3743955/1‑150 (MQ=255)
                                                                                                                     gTTTTTCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGc                             >  2:4736765/1‑150 (MQ=255)
                                                                                                                        tttCATCACACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAAt                          >  2:742275/1‑150 (MQ=255)
                                                                                                                          tcatcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATAc                        <  2:1395288/150‑1 (MQ=255)
                                                                                                                          tcatcaCACTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATAc                        <  1:4245785/150‑1 (MQ=255)
                                                                                                                             tcaCACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGcc                     <  1:1035287/150‑1 (MQ=255)
                                                                                                                               acacTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTATAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGATTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccct                   <  2:707286/150‑1 (MQ=255)
                                                                                                                                 acTGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                 >  1:2530304/1‑150 (MQ=255)
                                                                                                                                 acTGGTATTCCAGCGGATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAAGTAACAGTTTTTTCATAGTCGTATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctcc                 >  1:751654/1‑150 (MQ=255)
                                                                                                                                   tGGTATTCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAAGGATAAGTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGccctccct               >  1:4659852/1‑150 (MQ=255)
                                                                                                                                       aTTCCAGCGTATCGGTTTGCATGCGCTCCACCAGCTGATTAACGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATTATTAACAGTTTTTTCATAGTCATATCCCGAAGTCTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACg           <  1:2706677/150‑1 (MQ=255)
                                                                                                                                         tCCAGCGTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGtt         >  2:4434195/1‑150 (MQ=255)
                                                                                                                                               gTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt   >  2:1035949/1‑150 (MQ=255)
                                                                                                                                               gTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt   >  2:4968073/1‑150 (MQ=255)
                                                                                                                                               gTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt   >  2:5069824/1‑150 (MQ=255)
                                                                                                                                               gTATAGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgt   >  1:2905198/1‑150 (MQ=255)
                                                                                                                                                tATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAgtg  <  2:4793236/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGATTCCGTCAGTGTAACGCGCGCCAGAGGCTGAAATGCCCTGTTTGAGATGCAGTAGTTGATTATCGTAAACAAAACTGACCTCCTGGCGCGGATTATTCAGTTTGACCGTCAACGGTTTTTCATCACACTGGTATTCCAGCGTATCGGTTTGCATGCGCTCAACCAGCTGATTAAAGGCGCTACAGCCGGAGAGAAGGACCGGCAAAATGATAATTAACAGTTTTTTCATAGTCATATCCCGAAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTG  >  CP009273/1712415‑1712708

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