Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,973,968 T→A 54.0% G9G (GGA→GGT ygeA ← Asp/Glu_racemase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,973,9680TA54.0% 55.5 / inf 266G9G (GGA→GGTygeAAsp/Glu_racemase family protein
Reads supporting (aligned to +/- strand):  ref base T (66/55);  new base A (81/63);  total (147/119)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 283 total aligned reads displayed.

TTATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCG  >  CP009273/2973819‑2974117
                                                                                                                                                     |                                                                                                                                                     
ttATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCAtt                                                                                                                                                       >  2:1312548/1‑150 (MQ=255)
   gcccATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATtcct                                                                                                                                                    <  1:220964/149‑1 (MQ=255)
    cccATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAGGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATtcctc                                                                                                                                                   >  2:696411/1‑150 (MQ=255)
     ccATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATacctcc                                                                                                                                                  >  1:241206/1‑150 (MQ=255)
      cagtCCCCGCGACGCTGGCACTATTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTccc                                                                                                                                                 <  2:1783692/147‑1 (MQ=255)
      cATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTccc                                                                                                                                                 <  1:3043329/150‑1 (MQ=255)
       aTTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCa                                                                                                                                                <  1:947807/150‑1 (MQ=255)
         tCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGACTGGTGGATTCCCAGCTCATACCTCCCAGc                                                                                                                                              >  2:3861919/1‑150 (MQ=255)
            ccGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCaaa                                                                                                                                           <  1:4649044/150‑1 (MQ=255)
              gcgACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCACAcc                                                                                                                                         >  2:5014066/1‑150 (MQ=255)
              gcgACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAAcc                                                                                                                                         >  2:5159784/1‑150 (MQ=255)
               cgACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCa                                                                                                                                        >  2:4897494/1‑150 (MQ=255)
               cgACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCa                                                                                                                                        >  2:1859370/1‑150 (MQ=255)
                     tGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTTAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAAAGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGAAAACCAATTGtt                                                                                                                                  <  2:1672385/150‑1 (MQ=255)
                     tGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGtt                                                                                                                                  >  1:3826546/1‑150 (MQ=255)
                     tGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGtt                                                                                                                                  >  1:5231896/1‑150 (MQ=255)
                        cACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTc                                                                                                                               <  1:180389/150‑1 (MQ=255)
                         aCTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCa                                                                                                                              >  1:2571917/1‑150 (MQ=255)
                          ctctTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCAt                                                                                                                             <  1:2042226/150‑1 (MQ=255)
                          ctctTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCAt                                                                                                                             >  1:4253572/1‑150 (MQ=255)
                           tcttctATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCAtt                                                                                                                            >  1:389872/1‑150 (MQ=255)
                             ttctATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTcc                                                                                                                          <  2:2349101/150‑1 (MQ=255)
                              tctATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCct                                                                                                                         <  2:926929/150‑1 (MQ=255)
                               ctATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCctc                                                                                                                        >  2:2632821/1‑150 (MQ=255)
                               ctATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCctc                                                                                                                        >  1:1606741/1‑150 (MQ=255)
                                tATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCctct                                                                                                                       <  2:3620134/150‑1 (MQ=255)
                                tATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCctct                                                                                                                       >  2:4947540/1‑150 (MQ=255)
                                 aTTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCctctc                                                                                                                      >  2:1218169/1‑150 (MQ=255)
                                  tttCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCACACCAATTGTTTTCATTCctctct                                                                                                                     >  2:464551/1‑150 (MQ=255)
                                  tttCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCctctct                                                                                                                     <  1:2364988/150‑1 (MQ=255)
                                   ttCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCtctctc                                                                                                                    >  1:877592/1‑150 (MQ=255)
                                    tCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTcc                                                                                                                   <  2:1472573/150‑1 (MQ=255)
                                      aTGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTg                                                                                                                 >  2:3135436/1‑150 (MQ=255)
                                          aaaTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCGGAtgt                                                                                                             >  1:1013886/1‑150 (MQ=255)
                                          aaaTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGAtgt                                                                                                             >  1:1333926/1‑150 (MQ=255)
                                            aTCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGAtgtgt                                                                                                           >  2:3458629/1‑150 (MQ=255)
                                             tCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGAtgtgtg                                                                                                          <  1:990365/150‑1 (MQ=255)
                                             tCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGAtgtgtg                                                                                                          <  1:2908904/150‑1 (MQ=255)
                                              cGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATgtgtgt                                                                                                         >  2:329874/1‑150 (MQ=255)
                                               gACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGtt                                                                                                        >  1:883856/1‑150 (MQ=255)
                                               gACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGtt                                                                                                        >  1:3878806/1‑150 (MQ=255)
                                               gACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGtt                                                                                                        >  2:1287990/1‑150 (MQ=255)
                                               gACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGtt                                                                                                        <  1:1292019/150‑1 (MQ=255)
                                                aCGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTa                                                                                                       <  1:3630669/150‑1 (MQ=255)
                                                 cGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTAc                                                                                                      >  2:1588104/1‑150 (MQ=255)
                                                 cGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTAc                                                                                                      >  2:4518081/1‑150 (MQ=255)
                                                  gCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTTTTTTCATTCCTCTCTACTGATGTGTGTTAca                                                                                                     >  1:5119163/1‑149 (MQ=255)
                                                   cTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCg                                                                                                    >  1:4434268/1‑150 (MQ=255)
                                                     aTGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGca                                                                                                  >  1:5291451/1‑150 (MQ=255)
                                                      tGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGcaa                                                                                                 >  2:3722979/1‑150 (MQ=255)
                                                             gcaCTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCt                                                                                          <  1:420625/150‑1 (MQ=255)
                                                               aCTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCAGACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGa                                                                                        <  1:4981639/150‑1 (MQ=255)
                                                                 ttGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAAt                                                                                      <  1:371939/150‑1 (MQ=255)
                                                                 ttGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAAt                                                                                      <  2:1744830/150‑1 (MQ=255)
                                                                      cAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGa                                                                                 >  2:5299461/1‑150 (MQ=255)
                                                                      cAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGa                                                                                 <  2:4254004/150‑1 (MQ=255)
                                                                       agagTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGaa                                                                                >  1:457576/1‑150 (MQ=255)
                                                                             aaGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGAc                                                                          <  1:235502/150‑1 (MQ=255)
                                                                              aGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACa                                                                         >  1:3202386/1‑150 (MQ=255)
                                                                              aGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACa                                                                         >  2:5239894/1‑150 (MQ=255)
                                                                                cccACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACaaa                                                                       <  2:302451/150‑1 (MQ=255)
                                                                                 ccaccaAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAAc                                                                      <  1:4103112/150‑1 (MQ=255)
                                                                                 ccaccaAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAAc                                                                      <  1:5212203/150‑1 (MQ=255)
                                                                                 ccaccaAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAAc                                                                      <  1:1144847/150‑1 (MQ=255)
                                                                                      aaGCCTCTGTTTAATGCATTCATTTATCAAACGATTGTATGGATTGGTGGTTGCCCCGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGAAACATGGCTGAATAAGGAACGTGACAAACGCCTc                                                                 <  2:4469526/150‑1 (MQ=255)
                                                                                      aaGCCGCTGTTTAATTCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTc                                                                 <  2:1640550/150‑1 (MQ=255)
                                                                                      aaGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTc                                                                 <  2:1640684/150‑1 (MQ=255)
                                                                                        gCCGATGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCaa                                                               <  2:74935/150‑1 (MQ=255)
                                                                                         ccGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAAc                                                              <  2:660527/150‑1 (MQ=255)
                                                                                           gCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACgg                                                            >  1:1543273/1‑150 (MQ=255)
                                                                                           gCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACgg                                                            <  2:4709493/150‑1 (MQ=255)
                                                                                              gTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCa                                                         <  2:1594364/150‑1 (MQ=255)
                                                                                              gTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCa                                                         >  1:3103289/1‑150 (MQ=255)
                                                                                                  aaTGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTg                                                     >  1:3990123/1‑150 (MQ=255)
                                                                                                    tGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGcc                                                   >  2:2986237/1‑150 (MQ=255)
                                                                                                       cTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCaa                                                >  1:645590/1‑150 (MQ=255)
                                                                                                       cTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCaa                                                >  2:3149319/1‑150 (MQ=255)
                                                                                                         tCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAtt                                              <  2:491214/150‑2 (MQ=255)
                                                                                                          cATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAttt                                             >  2:3381941/1‑148 (MQ=255)
                                                                                                          aaTTTATAAAACGATTGTAAGTATTGGTGTATTCCAAGCACATACCTCCCAGCAAACAAATTGTTTTAATTCCTCTCTCGTTTTGTGTGTAACCGCAACAAGTCTTAATAAGTAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAttt                                             <  1:4942208/149‑3 (MQ=255)
                                                                                                           aTTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAAtttc                                            >  2:4598565/1‑150 (MQ=255)
                                                                                                             ttATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTa                                          >  1:3536841/1‑150 (MQ=255)
                                                                                                              tATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTAt                                         >  1:678985/1‑150 (MQ=255)
                                                                                                                   aaCGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACg                                    <  1:3458629/150‑1 (MQ=255)
                                                                                                                     cGATAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGaa                                  >  2:419682/1‑150 (MQ=255)
                                                                                                                       aTAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGaaaa                                >  1:4714047/1‑150 (MQ=255)
                                                                                                                       aTAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTTTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACATGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGaaaa                                <  1:1167209/150‑1 (MQ=255)
                                                                                                                       aTAGTAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGaaaa                                >  1:4143801/1‑150 (MQ=255)
                                                                                                                          gTAAGGAATGGTGGATTCCCAGCTCAGTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAg                             <  2:2536045/150‑1 (MQ=255)
                                                                                                                           tAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAgg                            <  2:4315309/150‑1 (MQ=255)
                                                                                                                           tAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAgg                            >  2:3742946/1‑150 (MQ=255)
                                                                                                                           tAAGGAATGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAgg                            >  1:5191764/1‑150 (MQ=255)
                                                                                                                             aGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCATATCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGc                          <  1:3654827/150‑1 (MQ=255)
                                                                                                                               gattggtgtATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTTTGTTACCGCAACAAGGCTGAATAAGGAACGAGAAAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGCCg                        <  1:2016262/147‑1 (MQ=255)
                                                                                                                               gAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGCCg                        >  1:3550097/1‑150 (MQ=255)
                                                                                                                                 aTGGTGGATTCCCAGCTCATACCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGCCGGa                      <  2:2695025/150‑1 (MQ=255)
                                                                                                                                        aTTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGCCGGATGTACAg               <  2:1544542/150‑1 (MQ=255)
                                                                                                                                              aGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATc         >  1:4528203/1‑150 (MQ=255)
                                                                                                                                               gCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATTTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc        <  2:877592/150‑1 (MQ=255)
                                                                                                                                                    ttCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTTATGTGTGTTTCCGCTACAAGGATGAATAATGAACGAGACAAACGCCTCTACGGCCAAGTGCCCTATTTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGTccc   <  2:3878884/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATCCCATTCCCCGCGACGCTGGCACTCTTCTATTTCATGAAAATCGACGCTATGTAGCAGCACTTGCGCAGAGTGAAGCCCACCAAGCCGCTGTTTAATGCCTTCATTTATCAAACGATAGTAAGGAATGGTGGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCG  >  CP009273/2973819‑2974117

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: