Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,153,615 G→C 5.9% S171C (TCC→TGC)  U712_16195 ← putative xanthine dehydrogenase subunit E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,153,6150GC 5.9% 544.0 / 4.5 245S171C (TCC→TGC) U712_16195putative xanthine dehydrogenase subunit E
Reads supporting (aligned to +/- strand):  ref base G (105/84);  new base C (3/22);  total (109/136)
Fisher's exact test for biased strand distribution p-value = 2.95e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

TGGCGTCTTTGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACAATGCTTTA  >  CP006881/3153468‑3153764
                                                                                                                                                   |                                                                                                                                                      
   cGTCTTTGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTcc                                                                                                                                                   <  1:5795723/150‑1 (MQ=255)
     tCTTTGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTTCGGGTTTTAAATGGCGGAATCCCg                                                                                                                                                 >  1:2305560/1‑150 (MQ=255)
       tttGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTGCGGCTTTTAAAAGCCGGACTCCCGTc                                                                                                                                               >  1:2193560/1‑150 (MQ=255)
       tttGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCCTc                                                                                                                                               >  2:4906436/1‑150 (MQ=255)
        ttGCTTTGTGTTCACTAATCTTTTTTTACTCCAGAAGCTGTCATGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCt                                                                                                                                              <  2:3931345/150‑1 (MQ=255)
        ttGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTCCCCCCAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGCCCGGCTTTTAAAAGCCGAACTCCCGTCt                                                                                                                                              <  2:4706015/150‑1 (MQ=255)
          gCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTAc                                                                                                                                            <  2:4578187/150‑1 (MQ=255)
          gCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTACGGCTTTTAAAAGCCGGACTCCCGTCTAc                                                                                                                                            >  2:2634339/1‑150 (MQ=255)
          gCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCAACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCAAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTAc                                                                                                                                            <  2:118432/150‑1 (MQ=255)
            tATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTTCGGCTTTTAATAGCCGGACTCCCGTCTAccc                                                                                                                                          >  2:4471461/1‑150 (MQ=255)
              tgtgTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTAccccc                                                                                                                                        <  1:722253/150‑1 (MQ=255)
              gctgtTCACTAATTTTTTTGTAATCCAGAACCGGTCAAGTTCTACGAGTTTTATCGTGTGTTTTCCCACAACGCTTGTGCGTCTTGATCGGGCTGACTGATTGAGGCTAAAAAAGTCCGGCTATAAAAATCCTGACTCCCGTCTAccccc                                                                                                                                        <  2:4952152/148‑1 (MQ=255)
               gtgtTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTCCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCg                                                                                                                                       <  2:1678478/150‑1 (MQ=255)
                   tCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGTGGCTAAAAAAGTTCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCTTTCa                                                                                                                                  >  1:367155/1‑150 (MQ=255)
                   tCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCa                                                                                                                                  <  2:2444856/150‑1 (MQ=255)
                    cACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCCCTCaa                                                                                                                                 >  1:5080864/1‑150 (MQ=255)
                      cTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTACCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGCCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAAct                                                                                                                               <  1:4844707/150‑1 (MQ=255)
                        aaTCTTTTTGTAATCCAGAAGCTGTCCCCGTTTANCCGATTTATGGTCGGTTTGCCCGCAACGCGTGTGGGTGTTGATTGCGCTGTCTGACTGAGGCTAAAAAAGTCCGACTTTTAAAAGCCGGTCTCCCCACCACCCCCGCTCAActct                                                                                                                             >  1:4672122/1‑150 (MQ=255)
                         ctCTTTTTGTAATCAAGAAGCTGCCAAGTTCTACGAGATTTAACGTGTGTTTTCCCACAACGCTTCTGCTTCTTGATCGCGCTGTCTGACTGAGGCAAAAAAAGTCCGGCTTTTAAAAGCCGTACTCCCGTCTACCCCCGTTCAActctc                                                                                                                            <  1:3534142/149‑1 (MQ=255)
                         aTCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGTCTGAGGCTAAAAAAGTCCGGCTTTTAAAAACCGCAATAGCCTCTACCCCCATTCAActctc                                                                                                                            >  1:2619621/1‑150 (MQ=255)
                          ccTTTTTGTAACCAAAAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTCCCCACACCGTTTCTGCTTCTTGATCGCGCTGTCTGTTTGAGGCTAAAAAATTCCGCCTTTTAAAACCCGTACTCCCGTCTACCCCCGTTCAActctct                                                                                                                           <  2:1844103/149‑1 (MQ=255)
                           cTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTACGGCTTTTTAAAGCCGGACTCCCGGCTACCCCCCTTCAACCCtctt                                                                                                                          >  1:1100090/1‑150 (MQ=255)
                            tttttGTAATCCAGAAGCTGTCAAGTTCTACGAGTTTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCGGGACTCCCGTCTACCCCAGTTCAACTCtcttc                                                                                                                         <  1:3936154/150‑1 (MQ=255)
                            tttttGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGGCTCCCGTCTACCCCCGTTCAACTCtcttc                                                                                                                         >  1:672558/1‑150 (MQ=255)
                            tttttGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCCTTCAACCCtcttc                                                                                                                         >  1:4356827/1‑150 (MQ=255)
                                tGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTCCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAAt                                                                                                                     <  2:2486938/150‑1 (MQ=255)
                                tGTAACCCAGAAGCAGCCAAGTTCTCCGAGANTTACTGTGTGTTTACTCAAAACGCTTCTGCGCCTTGCTCGCGCTGTCTGACTGAGGCAAAAAAAGTCCGCTTTTTAAAAGCCTAACTTCCGTCTACCCCCGTTCAACTCTCTTCTAAt                                                                                                                     <  1:2096949/150‑1 (MQ=255)
                                  tAACCCAGAACCTCCCACGTTCAACAATCTTCATCGTCTTTTTTCCCCCTCCGCTTCTCCTTCTTGACCGCGCTGTCTGACTGAGCGTAAAAAAGCCCGTCTTTTAAAACCCCCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATcc                                                                                                                   <  1:5544809/150‑1 (MQ=255)
                                   acaccAGAAGCTCACAAGTTCTAGGAGATTTATCGTGTGTTTCCCCCCACCGCTCCTGCGTCTTGATCGCGCTGTCTTACTGAGGCAAAAAAAGTCCGGCTTTTAAAAGCCCAACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCt                                                                                                                  <  1:1793427/147‑1 (MQ=255)
                                   aaTCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCTGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCt                                                                                                                  <  2:2767965/150‑1 (MQ=255)
                                   aaTCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCt                                                                                                                  <  2:4004000/150‑1 (MQ=255)
                                     tCCAGAAGCTGTCAAGTCCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGc                                                                                                                <  2:592713/150‑1 (MQ=255)
                                        aGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAg                                                                                                             >  2:4668948/1‑150 (MQ=255)
                                         gAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTGCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGc                                                                                                            >  2:2801141/1‑150 (MQ=255)
                                         gAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGc                                                                                                            >  2:3079349/1‑150 (MQ=255)
                                          aaGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGGGCTGTCTGACTGAGGCAAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCt                                                                                                           <  2:5181516/150‑1 (MQ=255)
                                          aaGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCt                                                                                                           <  1:1373410/150‑1 (MQ=255)
                                          aaGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCAAAAAAAGTCCGGCTTTTAAAAGCCCGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCt                                                                                                           <  1:4748629/150‑1 (MQ=255)
                                           aGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTg                                                                                                          <  2:5359372/150‑1 (MQ=255)
                                                  aaGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCAt                                                                                                   <  2:3834840/150‑1 (MQ=255)
                                                   aGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGCGTTTTAAAAGCCGGGCTCTCGTCAACGCCCCTTCAACCCCCTTCTAATACTGCAAGCTGCCCGCaga                                                                                                  >  1:1602288/1‑148 (MQ=255)
                                                     ttCTAGGAGTTTTATCTTGTGTTTTCCCACAACTCTTCTGCGTTTTGATCCCGCTGTCTGACTGAGGCAAAAAAAGTCCGGCTTTTAAAAGCCGAACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAAt                                                                                                <  1:325072/150‑1 (MQ=255)
                                                     ttCTAGGAGATTTATCGTGTTTTTTCCCACAACGCTTCTGCTTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAAt                                                                                                <  1:4576365/150‑1 (MQ=255)
                                                      tCTACGAGTGTTATCGTGTGTTTTCCAACAACGTTTCTGCGTCTTGATCCCGCTGTCTGATTGAGGTTAAACAAGTCCGCCCTTTAAAAGCCTTACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATc                                                                                               <  1:917017/150‑1 (MQ=255)
                                                       cTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTACGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCCTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAAAcc                                                                                              >  2:1613223/1‑150 (MQ=255)
                                                        tACGAGATTTATCGTGTGTTTTCCCACACCGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATccc                                                                                             <  2:5316672/150‑1 (MQ=255)
                                                        tACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGACGGACTCCCCTCTTCCCCCGCTTCAACCTCTTCCAATCCCTCAAGCTAGCCCCATAATccc                                                                                             >  1:3381962/1‑150 (MQ=255)
                                                        tACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTACGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGGAAGCTGACCGCATAATccc                                                                                             >  2:3029845/1‑150 (MQ=255)
                                                            agaTTTACTGTGTGTTTTCCCCCAACTCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCAGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCa                                                                                         <  2:4856085/150‑1 (MQ=255)
                                                             gaTTTATCTGGTGTTTTCCCACACCGCTTCAGCGTCTTGATCGAGCTGTCTGTCTGAGCCTAAAAAAGTCGGGCTTTTAAAACCCAGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCAt                                                                                        <  2:3923683/150‑1 (MQ=255)
                                                               tttATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAATTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATAc                                                                                      <  2:777863/150‑1 (MQ=255)
                                                               tttATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATAc                                                                                      >  2:4344651/1‑150 (MQ=255)
                                                                     gtgtttTTTCCCACAACGCTTTTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTa                                                                                <  1:2828458/150‑1 (MQ=255)
                                                                     gtgtctTTTCCCACAAGGCTTCTGCGTCTTGATCGCCCTGTTTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCAGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTa                                                                                <  2:4319195/150‑1 (MQ=255)
                                                                      tgtgtTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTac                                                                               <  2:4612423/150‑1 (MQ=255)
                                                                             tCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTTTTCTAATCCTTCAAGCTGGCCACATAATCCCGCCATACCCCGTACACCGGca                                                                        >  2:2089752/1‑150 (MQ=255)
                                                                                   accGCTTCTGTGTTTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAAtt                                                                  <  2:189689/148‑1 (MQ=255)
                                                                                     cGCTTCTGCGCCTTGATCGTGCTGTCCACCTGAGGCAAAAAAAGTCCGGCTTTTAAAACCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAACCCCGCCATACCCGGTACACCGGCACAAATTcc                                                                <  2:5620871/150‑1 (MQ=255)
                                                                                      gCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATAACCGGGACACCGGCACAAATTccc                                                               >  1:3602729/1‑150 (MQ=255)
                                                                                      gCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTccc                                                               <  2:3473487/150‑1 (MQ=255)
                                                                                       cTTCTGCTTCTTGATCGCGCTGTCTGATTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCCTACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTcccc                                                              <  1:4996896/150‑1 (MQ=255)
                                                                                        ttCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTTCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCCCCATACCCGGTACACCCGCACAAATTCCCCg                                                             >  1:5155032/1‑150 (MQ=255)
                                                                                        ttCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCTGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCg                                                             <  2:4642929/150‑1 (MQ=255)
                                                                                        ttCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCg                                                             <  2:5429655/150‑1 (MQ=255)
                                                                                         tCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGc                                                            <  2:2290129/150‑1 (MQ=255)
                                                                                               tCTTGATCGCGCTGTCTGACTGAGGCTAAATAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGccc                                                      <  2:2195556/150‑1 (MQ=255)
                                                                                                cTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCct                                                     <  2:1051737/150‑1 (MQ=255)
                                                                                                   gATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAAATTCCCCCCAGCcctcct                                                  >  1:3391476/1‑150 (MQ=255)
                                                                                                       gcgcTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCta                                              <  2:2697119/150‑1 (MQ=255)
                                                                                                        cgcTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCtat                                             <  2:3749115/150‑1 (MQ=255)
                                                                                                          cTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACATCGGCACAAATTCCCCGCCAGCCCCTCTTCTatat                                           >  1:3730035/1‑150 (MQ=255)
                                                                                                           tgcctgactgaGGCTAAAAAAGTCCGGCTTTTAAAAGCCGTACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATc                                          <  2:5711731/147‑1 (MQ=255)
                                                                                                           tGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATc                                          <  1:2468595/150‑1 (MQ=255)
                                                                                                             actgactgaGCCTAAAAAACTGGGGCTTTTAAAACCCCTACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCtt                                        <  2:5265081/149‑1 (MQ=255)
                                                                                                                 actgaGGCTAAAAAAGTCCGGCTTTTAGAACCCCTACTTCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTtctc                                    <  2:862161/150‑1 (MQ=255)
                                                                                                                  ctgaGGCTAAAAAAGTGCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGAACACCGGCACAAAATCCCCGCCAGCCCTTCTTCTATATCTTTGTCa                                   >  1:5450199/1‑150 (MQ=255)
                                                                                                                  ctgaGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCa                                   <  1:4190560/150‑1 (MQ=255)
                                                                                                                    gaGGCTAAAAAAGCCAGCCTTTTAAAAGCCCCACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGa                                 <  1:3289609/150‑1 (MQ=255)
                                                                                                                       ggaaaaaaaCTCCGGCTTTTAAAAGCCCGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAgg                              <  1:84326/147‑1 (MQ=255)
                                                                                                                       gCTAAAAAAGTCCGGCTTTTAAAAGCCGAACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAgg                              <  1:4656073/150‑1 (MQ=255)
                                                                                                                        cTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGt                             <  2:2038989/150‑1 (MQ=255)
                                                                                                                        cTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGt                             >  2:4791380/1‑150 (MQ=255)
                                                                                                                          aaaaaaGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTg                           <  2:5638051/150‑1 (MQ=255)
                                                                                                                             caaGTCCGGCTTTTAAAAGCCGCACTTCCGTCTACCCCTGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAAACCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAgg                        <  2:2305560/149‑1 (MQ=255)
                                                                                                                               aGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAgggg                      >  2:2902928/1‑150 (MQ=255)
                                                                                                                                gTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGGAGGTTGGGGGGt                     >  2:3017772/1‑150 (MQ=255)
                                                                                                                                gTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGt                     >  2:737125/1‑150 (MQ=255)
                                                                                                                                 tCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGtt                    >  2:252649/1‑150 (MQ=255)
                                                                                                                                 tCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGtt                    >  1:1001053/1‑150 (MQ=255)
                                                                                                                                   cGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTc                  >  2:3628000/1‑150 (MQ=255)
                                                                                                                                     gCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTccc                <  1:4325129/150‑1 (MQ=255)
                                                                                                                                        tttAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGaa             >  2:2248133/1‑150 (MQ=255)
                                                                                                                                         ttAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGaaa            >  2:3223336/1‑150 (MQ=255)
                                                                                                                                           aaaaGCCGTACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACa          <  2:1100090/150‑1 (MQ=255)
                                                                                                                                           aaaaGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACa          <  2:451604/150‑1 (MQ=255)
                                                                                                                                            aaaGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGTCCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACaa         >  2:800221/1‑150 (MQ=255)
                                                                                                                                               gCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACAATGc      <  2:448216/150‑1 (MQ=255)
                                                                                                                                                ccGTACTCCCGTCTACCCCCGTTCAACTCTTTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACAATGCt     <  1:1672351/150‑1 (MQ=255)
                                                                                                                                                 cGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACAATGCtt    >  1:2731300/1‑150 (MQ=255)
                                                                                                                                                   gACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACAATGCTTTa  <  2:3388586/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TGGCGTCTTTGCTATGTGTTCACTAATCTTTTTGTAATCCAGAAGCTGTCAAGTTCTACGAGATTTATCGTGTGTTTTCCCACAACGCTTCTGCGTCTTGATCGCGCTGTCTGACTGAGGCTAAAAAAGTCCGGCTTTTAAAAGCCGGACTCCCGTCTACCCCCGTTCAACTCTCTTCTAATCCTGCAAGCTGACCGCATAATCCCGCCATACCCGGTACACCGGCACAAATTCCCCGCCAGCCCTTCTTCTATATCTTTGTCAGAAGGTTGAGGGGTTTCCCGAAACAATGCTTTA  >  CP006881/3153468‑3153764

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