Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,836,192 T→A 100% I67I (ATA→ATT ynjH ← DUF1496 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,836,1920TA97.5% 660.9 / 4.1 203I67I (ATA→ATTynjHDUF1496 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base A (94/104);  total (95/108)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 203 total aligned reads displayed.

GCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGG  >  CP009273/1836046‑1836336
                                                                                                                                                  |                                                                                                                                                
 ctggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTc                                                                                                                                              >  1:1543459/1‑150 (MQ=255)
 ctggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCACCAGAATTc                                                                                                                                              >  2:3534801/1‑150 (MQ=255)
  tggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTcc                                                                                                                                             >  1:4368779/1‑150 (MQ=255)
   ggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCt                                                                                                                                            >  2:1976080/1‑150 (MQ=255)
    gtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCtt                                                                                                                                           <  1:3170833/150‑1 (MQ=255)
     tctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTc                                                                                                                                          <  2:3684101/150‑1 (MQ=255)
      ctggtGCTCACCTGAAGAGTCGCTGCACTATCCGCGTGCCCCTGCTGACTTTCCATTATTAGAGGCGTTTATGGCTTTAAGCGCCGCCAGAACAGCGGATTAGTGCTAATGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCATTCg                                                                                                                                         <  2:1327315/150‑1 (MQ=255)
       tggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATTCCTTCgg                                                                                                                                        <  1:1172431/150‑1 (MQ=255)
        ggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTCTTAGAGGCGTTTATGGCTTTACCCGCCGCCAGACCACCGGCTTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCCGCAGAATTCCATCGGc                                                                                                                                       >  2:749772/1‑150 (MQ=255)
          tGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCtt                                                                                                                                     <  2:2204079/150‑1 (MQ=255)
           gCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCttt                                                                                                                                    <  2:812840/150‑1 (MQ=255)
            cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                                   <  1:1464810/150‑1 (MQ=255)
            cTCACCTGAAGAGGCGCTGCAATATCCACTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                                   <  1:1927484/150‑1 (MQ=255)
               aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAtt                                                                                                                                <  1:775584/150‑1 (MQ=255)
               aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAtt                                                                                                                                <  1:775696/150‑1 (MQ=255)
                  tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTAcc                                                                                                                             <  1:5000856/150‑1 (MQ=255)
                  tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTAcc                                                                                                                             <  2:1563598/150‑1 (MQ=255)
                  tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTAcc                                                                                                                             >  2:2149574/1‑150 (MQ=255)
                  tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTAcc                                                                                                                             <  1:1084529/150‑1 (MQ=255)
                   gAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCg                                                                                                                            >  1:2587401/1‑150 (MQ=255)
                    aaGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTTCGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCGTTAATTACCGc                                                                                                                           >  1:2327048/1‑150 (MQ=255)
                       agGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCCATGCAGCAGAATTCCTTAGGCTTTAATTACCGcccc                                                                                                                        >  1:2625083/1‑150 (MQ=255)
                       agGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGcccc                                                                                                                        <  2:4858958/150‑1 (MQ=255)
                        ggCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCt                                                                                                                       <  1:1673697/150‑1 (MQ=255)
                          cgcTGCAATATCCGCTGTCCACTGCTGACATTCCATTATTAGATGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGTTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTc                                                                                                                     <  2:1634742/150‑1 (MQ=255)
                          cgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTc                                                                                                                     <  1:3327967/150‑1 (MQ=255)
                           gcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCa                                                                                                                    <  1:1039821/150‑1 (MQ=255)
                              gNAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTTCGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAgta                                                                                                                 >  1:183044/1‑148 (MQ=255)
                                   atCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                                                                            >  2:2830138/1‑150 (MQ=255)
                                      cGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTg                                                                                                         >  2:4057041/1‑150 (MQ=255)
                                         tGGCCCCTGCTGACATTCCATTATGAGAGGCGTTTATGGCTTTTCGCGCCGCCAGACCAGCGGATTAGTGCTGAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTGATTTACCGCCCCTTCCGCATAATGTTGAtg                                                                                                      >  2:767321/1‑149 (MQ=255)
                                          ggCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATcc                                                                                                     <  1:4632081/150‑1 (MQ=255)
                                           gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATTGCTTTACGCCCCCCCAGACAACCGGATTAGTTCTAAGGGTTTTGTCAACACGCTTGCCTTGCACCAGAATTCCTTCGGCGTTAAATAACGCCCCTTCAGAAAAAATTTTATACt                                                                                                    >  1:3465169/1‑150 (MQ=255)
                                           gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCATAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATcca                                                                                                    <  1:2816855/150‑2 (MQ=255)
                                           gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCt                                                                                                    <  1:670475/150‑1 (MQ=255)
                                           gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCGTCAGAATAATTTTGATCCt                                                                                                    >  2:1460657/1‑150 (MQ=255)
                                            ccccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTg                                                                                                   <  2:3512387/150‑1 (MQ=255)
                                              ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGGCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGAt                                                                                                 <  2:3958706/150‑1 (MQ=255)
                                              ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGAt                                                                                                 <  2:3201164/150‑1 (MQ=255)
                                                tgctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATaa                                                                                               <  2:1084986/150‑1 (MQ=255)
                                                  ctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAAc                                                                                             <  2:4115759/150‑1 (MQ=255)
                                                          ccATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTAGGCTTTAATTAACGCCCCATCAGAATAATTTTGATCCTGATAAAAGCAGCACt                                                                                     >  1:3098626/1‑150 (MQ=255)
                                                            attattAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTTCAGCAGAATTACTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTTATCCAGATAAAAACAGCAcata                                                                                   >  1:5077686/1‑147 (MQ=255)
                                                             ttattaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGTGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAg                                                                                  <  2:680155/150‑1 (MQ=255)
                                                             ttattaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAg                                                                                  <  1:3740025/150‑1 (MQ=255)
                                                                  agagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAg                                                                             >  2:1730599/1‑150 (MQ=255)
                                                                    agGCGTTTATGGCTTTTCGCGCCGCCAGACCAGCGGATTTGTGCTATGGGTTTTGTCATCACGCTGGCATTGCTGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAAAGCAGCACTGAGTTCAggg                                                                           <  1:3013258/150‑1 (MQ=255)
                                                                      gCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCt                                                                         >  2:1888246/1‑150 (MQ=255)
                                                                       cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTg                                                                        >  2:1173429/1‑150 (MQ=255)
                                                                       cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTg                                                                        >  1:1268292/1‑150 (MQ=255)
                                                                          ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgcgc                                                                     >  2:167222/1‑150 (MQ=255)
                                                                              ggCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGAc                                                                 >  1:508228/1‑150 (MQ=255)
                                                                                cTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGNGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                               >  1:2680974/1‑150 (MQ=255)
                                                                                cTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                               >  1:91696/1‑150 (MQ=255)
                                                                                cTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                               >  1:1458227/1‑150 (MQ=255)
                                                                                cTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                               <  2:774356/150‑1 (MQ=255)
                                                                                   tACGCGCCGCCAGACTAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                            <  2:2387173/150‑1 (MQ=255)
                                                                                   tACGCGCCGCCAGACCAGCGGATTAGTGCTACGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                            <  2:4270756/150‑1 (MQ=255)
                                                                                   tACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGccc                                                            <  2:3983347/150‑1 (MQ=255)
                                                                                      gcgcCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGcc                                                         >  1:1155499/1‑150 (MQ=255)
                                                                                        gccgccAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctg                                                       >  1:3094292/1‑150 (MQ=255)
                                                                                          cgccAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCGGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaa                                                     <  2:1004745/150‑1 (MQ=255)
                                                                                            ccAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaact                                                   <  1:3815099/150‑1 (MQ=255)
                                                                                             cAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaactg                                                  >  1:4102599/1‑150 (MQ=255)
                                                                                              aGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaactga                                                 >  1:4141348/1‑150 (MQ=255)
                                                                                               gACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaactgaa                                                >  2:2364580/1‑150 (MQ=255)
                                                                                                  cAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAacac                                             >  1:4658696/1‑150 (MQ=255)
                                                                                                       gATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTg                                        <  1:1571314/150‑1 (MQ=255)
                                                                                                        aTTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTgg                                       >  1:4902695/1‑150 (MQ=255)
                                                                                                        aTTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTgg                                       <  1:4885142/150‑1 (MQ=255)
                                                                                                        aTTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTgg                                       <  2:4816078/150‑1 (MQ=255)
                                                                                                          tAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCg                                     >  2:383314/1‑150 (MQ=255)
                                                                                                            gTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGt                                   <  1:1592562/150‑1 (MQ=255)
                                                                                                            gTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGt                                   >  1:1218536/1‑150 (MQ=255)
                                                                                                              gCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTAc                                 >  2:3928339/1‑150 (MQ=255)
                                                                                                                tAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGt                               >  1:2367746/1‑150 (MQ=255)
                                                                                                                    ggTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTAc                           >  2:4560304/1‑150 (MQ=255)
                                                                                                                        ttGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTcc                       <  1:938359/150‑1 (MQ=255)
                                                                                                                          gTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCAc                     >  2:4978642/1‑150 (MQ=255)
                                                                                                                           tcatcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACg                    <  1:2082273/150‑1 (MQ=255)
                                                                                                                             atcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTc                  <  1:1990567/150‑1 (MQ=255)
                                                                                                                             atcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTc                  >  1:2932411/1‑150 (MQ=255)
                                                                                                                              tgaCGCTGGCAGTGCAGCAGATTTCCTTCGGCTTTTATTACCGACGCTTCAGAATTATTTTGCTCCTGATAAACGCAGCACTGAGTACAGGGCTGCGATGACTGCCCGCCTGAACTGATCACTTCTGGCGGTTCGTTTACCTACACGTcc                 <  2:3230584/148‑1 (MQ=255)
                                                                                                                              tcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTCCCTCCACGTcc                 >  2:284998/1‑150 (MQ=255)
                                                                                                                                   cTGGCATTGCAGCAGAATTCCTTCGGCTTTACTTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACg            <  2:2693060/150‑1 (MQ=255)
                                                                                                                                   cTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACg            <  1:1722235/150‑1 (MQ=255)
                                                                                                                                    tGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGCATAATTTTGATCCTGATAAACGCAGCACTGAGTCCAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGa           >  2:903011/1‑150 (MQ=255)
                                                                                                                                     agCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGAt          >  1:4033474/2‑150 (MQ=255)
                                                                                                                                      gCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGCCCTGACCACTTCTGGCGGTACGTTTACCGCCCCGTCCGGACGATa         >  2:3368476/1‑150 (MQ=255)
                                                                                                                                      gCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATa         <  2:2757563/150‑1 (MQ=255)
                                                                                                                                      gCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATa         >  2:534018/1‑150 (MQ=255)
                                                                                                                                       cATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATaa        <  1:1375901/150‑1 (MQ=255)
                                                                                                                                       cATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATaa        >  2:1803802/1‑150 (MQ=255)
                                                                                                                                       cATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATaa        >  1:3670562/1‑150 (MQ=255)
                                                                                                                                          tGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGc     <  2:3094292/150‑1 (MQ=255)
                                                                                                                                          tGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGc     >  1:4533012/1‑150 (MQ=255)
                                                                                                                                            cagcagAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCgg   <  1:2149574/150‑1 (MQ=255)
                                                                                                                                            cagcagAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCgg   <  1:4971516/150‑1 (MQ=255)
                                                                                                                                            cagcagAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCgg   >  1:3142595/1‑150 (MQ=255)
                                                                                                                                             agcagAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCggg  >  1:1093644/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGG  >  CP009273/1836046‑1836336

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