Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,577,095 G→T 37.1% V110V (GTG→GTT mrr → methylated adenine and cytosine restriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,577,0950GT37.1% 407.7 / 281.3 277V110V (GTG→GTTmrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base G (71/103);  new base T (37/65);  total (109/168)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 290 total aligned reads displayed.

TGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTC  >  CP009273/4576946‑4577243
                                                                                                                                                     |                                                                                                                                                    
    cgcATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATg                                                                                                                                                  <  2:3977955/150‑1 (MQ=255)
     gcATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGt                                                                                                                                                 <  2:4156756/150‑1 (MQ=255)
      cATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTc                                                                                                                                                <  1:3859295/150‑1 (MQ=255)
      cATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTc                                                                                                                                                >  1:2478746/1‑150 (MQ=255)
       aTGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCa                                                                                                                                               <  2:2674789/150‑1 (MQ=255)
       aTGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCa                                                                                                                                               <  1:3245106/150‑1 (MQ=255)
       aTGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCa                                                                                                                                               >  1:1838817/1‑150 (MQ=255)
       aTGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGACTTCGCTTTTGTTAATGTCa                                                                                                                                               >  2:3825753/1‑150 (MQ=255)
        tGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCaa                                                                                                                                              >  1:3072263/1‑150 (MQ=255)
        tGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCaa                                                                                                                                              >  2:3616754/1‑150 (MQ=255)
        tGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGGTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCaa                                                                                                                                              <  2:1354303/150‑1 (MQ=255)
         gACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCaaa                                                                                                                                             >  2:338643/1‑150 (MQ=255)
          aCCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAAc                                                                                                                                            <  1:1633878/150‑1 (MQ=255)
             gTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTa                                                                                                                                         <  2:1829911/150‑1 (MQ=255)
               ttAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAg                                                                                                                                       <  1:2303332/150‑1 (MQ=255)
               ttAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAg                                                                                                                                       <  2:1261508/150‑1 (MQ=255)
                 aaaaCGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTc                                                                                                                                     <  2:3871546/150‑1 (MQ=255)
                    aCGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACg                                                                                                                                  <  2:3998638/150‑1 (MQ=255)
                     cGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACgg                                                                                                                                 <  1:786316/150‑1 (MQ=255)
                      gTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGc                                                                                                                                >  1:1403830/1‑150 (MQ=255)
                        gCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCg                                                                                                                              <  2:1765745/150‑1 (MQ=255)
                          cGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGa                                                                                                                            >  1:2753120/1‑150 (MQ=255)
                                gTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGt                                                                                                                      <  1:2684296/150‑1 (MQ=255)
                                 tCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTc                                                                                                                     >  2:203410/1‑150 (MQ=255)
                                     aaaGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGAtt                                                                                                                 <  2:3865910/150‑1 (MQ=255)
                                      aaGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGAttt                                                                                                                <  2:851189/150‑1 (MQ=255)
                                       aGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTa                                                                                                               <  2:1147293/150‑1 (MQ=255)
                                       aGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTa                                                                                                               >  1:3698254/1‑150 (MQ=255)
                                        gTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAg                                                                                                              >  1:3047682/1‑150 (MQ=255)
                                          ttGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGAt                                                                                                            <  2:3729050/150‑1 (MQ=255)
                                            gTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATcc                                                                                                          <  1:1078692/150‑1 (MQ=255)
                                             tCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCg                                                                                                         >  1:4442732/1‑150 (MQ=255)
                                             tCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCg                                                                                                         <  2:2478746/150‑1 (MQ=255)
                                                cgTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGaaa                                                                                                      <  1:4270977/150‑1 (MQ=255)
                                                    gCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCg                                                                                                  <  1:1017212/150‑1 (MQ=255)
                                                    gCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCg                                                                                                  <  2:770896/150‑1 (MQ=255)
                                                       aaTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTTATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGTAAGCCGAct                                                                                               <  1:211794/150‑1 (MQ=255)
                                                       aaTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGAct                                                                                               >  1:1736604/1‑150 (MQ=255)
                                                       aaTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGAct                                                                                               <  2:4025671/150‑1 (MQ=255)
                                                              ccTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGa                                                                                        <  2:265961/150‑1 (MQ=255)
                                                              ccTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGa                                                                                        <  1:776197/150‑1 (MQ=255)
                                                                 gACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCa                                                                                     <  1:704708/150‑1 (MQ=255)
                                                                 gACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCa                                                                                     <  1:623287/150‑1 (MQ=255)
                                                                  aCTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCAt                                                                                    >  2:3464179/1‑150 (MQ=255)
                                                                  aCTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCAt                                                                                    >  2:3168417/1‑150 (MQ=255)
                                                                          gggTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACt                                                                            >  1:3460505/1‑150 (MQ=255)
                                                                          gggTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACt                                                                            >  2:3510614/1‑150 (MQ=255)
                                                                           ggTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACtt                                                                           <  2:4363158/150‑1 (MQ=255)
                                                                             ttttGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGc                                                                         <  2:2858917/150‑1 (MQ=255)
                                                                                tGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCaaa                                                                      >  2:887283/1‑150 (MQ=255)
                                                                                  aCTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAaga                                                                    <  1:3790030/150‑1 (MQ=255)
                                                                                   cTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAagag                                                                   <  1:1480182/150‑1 (MQ=255)
                                                                                      ggTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGagcag                                                                <  2:1907114/150‑1 (MQ=255)
                                                                                         tGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGccc                                                             <  1:887918/150‑1 (MQ=255)
                                                                                         tGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGccc                                                             <  1:3616754/150‑1 (MQ=255)
                                                                                            gTCTCATCCCCAGCCAATGACGGAGCAGGAGACGACCCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGa                                                          <  1:843851/150‑1 (MQ=255)
                                                                                            gTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGa                                                          <  1:4195816/150‑1 (MQ=255)
                                                                                              ctcATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACg                                                        >  2:111669/1‑150 (MQ=255)
                                                                                               tcATCCCCAGCCAATGAGGGTGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGa                                                       <  2:1545579/150‑1 (MQ=255)
                                                                                                cATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGAt                                                      >  1:3828361/1‑150 (MQ=255)
                                                                                                cATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGAt                                                      >  2:3209580/1‑150 (MQ=255)
                                                                                                    cccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCTATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGt                                                  <  2:2610448/150‑1 (MQ=255)
                                                                                                     ccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGtt                                                 <  1:2236604/150‑1 (MQ=255)
                                                                                                     ccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGtt                                                 <  1:440348/150‑1 (MQ=255)
                                                                                                      cAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTa                                                <  2:3221926/150‑1 (MQ=255)
                                                                                                       aGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAg                                               >  1:1312142/1‑150 (MQ=255)
                                                                                                       aGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAg                                               >  1:3257100/1‑150 (MQ=255)
                                                                                                          cAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATc                                            <  1:2159959/150‑1 (MQ=255)
                                                                                                           aaTGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCa                                           >  2:4399614/1‑150 (MQ=255)
                                                                                                            aTGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAACCTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGATAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAg                                          <  2:2712497/150‑1 (MQ=255)
                                                                                                                      aGGAGACGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAag                                >  2:4529642/1‑150 (MQ=255)
                                                                                                                         agaCGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGc                             <  1:1279653/150‑1 (MQ=255)
                                                                                                                          gaCGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCt                            <  2:889060/150‑1 (MQ=255)
                                                                                                                            cGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCTCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTc                          <  2:4334355/150‑1 (MQ=255)
                                                                                                                            cGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTc                          >  1:1826163/1‑150 (MQ=255)
                                                                                                                            cGAACCATCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTc                          >  1:4156336/1‑150 (MQ=255)
                                                                                                                              aaCCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGt                        <  1:1924425/150‑1 (MQ=255)
                                                                                                                              aaCCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGt                        <  1:4586177/150‑1 (MQ=255)
                                                                                                                               aCCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTg                       >  1:2896066/1‑150 (MQ=255)
                                                                                                                                ccATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGa                      <  2:2150593/150‑1 (MQ=255)
                                                                                                                                 cATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGAt                     >  1:1055036/1‑150 (MQ=255)
                                                                                                                                  aTCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATg                    <  2:3417764/150‑1 (MQ=255)
                                                                                                                                   tCTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGc                   >  2:657348/1‑150 (MQ=255)
                                                                                                                                   tCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGc                   <  1:4358863/150‑1 (MQ=255)
                                                                                                                                    cTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCg                  >  2:591011/1‑150 (MQ=255)
                                                                                                                                      ggCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGt                >  2:2967199/1‑150 (MQ=255)
                                                                                                                                        ccTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTgg              >  1:3981341/1‑150 (MQ=255)
                                                                                                                                        ccTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTgg              <  1:1900533/150‑1 (MQ=255)
                                                                                                                                         cTTCGCTTGTGTGAATGTCAACCTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCAGGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGGGGc             <  2:450164/150‑1 (MQ=255)
                                                                                                                                          ttCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCt            <  2:4033010/150‑1 (MQ=255)
                                                                                                                                          ttCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCt            >  1:2844243/1‑150 (MQ=255)
                                                                                                                                            cGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCtga          >  1:2944271/1‑150 (MQ=255)
                                                                                                                                             gCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCtgat         >  2:4602840/1‑150 (MQ=255)
                                                                                                                                                tttGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAg      >  2:1026281/1‑150 (MQ=255)
                                                                                                                                                 ttGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAgg     >  2:1200491/1‑150 (MQ=255)
                                                                                                                                                  tgtgAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGt    <  1:3292237/150‑1 (MQ=255)
                                                                                                                                                  tgtgAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGt    <  2:3840272/150‑1 (MQ=255)
                                                                                                                                                  tgtgAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGt    <  2:1598725/150‑1 (MQ=255)
                                                                                                                                                   gtgAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGtt   >  2:4519661/1‑150 (MQ=255)
                                                                                                                                                   gtgAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGtt   >  2:1631547/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTC  >  CP009273/4576946‑4577243

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: