Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,135 T→C 100% Y365C (TAC→TGC)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,1350TC99.6% 932.9 / ‑4.8 233Y365C (TAC→TGC) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (123/108);  total (125/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 233 total aligned reads displayed.

GCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTC  >  CP009273/2304988‑2305281
                                                                                                                                                   |                                                                                                                                                  
    cgcgGGTAAAGCTTATTTCCATACTTGATTATCCTCATGCGAACGGTCTCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCa                                                                                                                                              <  1:1185361/150‑1 (MQ=255)
    cgcgGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCa                                                                                                                                              <  2:3675625/150‑1 (MQ=255)
     gcgGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAg                                                                                                                                             >  1:4451/1‑150 (MQ=255)
     gcgGGAAAAGCTTATTTCCATACGTGATTATCCTCAGGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAg                                                                                                                                             <  1:3496576/150‑1 (MQ=255)
       gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAc                                                                                                                                           <  2:1879752/150‑1 (MQ=255)
       gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAc                                                                                                                                           <  1:3814307/150‑1 (MQ=255)
       gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAc                                                                                                                                           <  2:1879774/150‑1 (MQ=255)
        ggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAcc                                                                                                                                          <  2:3891652/150‑1 (MQ=255)
         gTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGAccc                                                                                                                                         >  2:2945565/1‑150 (MQ=255)
          tAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCa                                                                                                                                        >  1:2033683/1‑150 (MQ=255)
          tAAAGCTTATTTCAATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACAAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCa                                                                                                                                        <  1:827960/150‑1 (MQ=255)
              gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGATTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGc                                                                                                                                    <  1:4055356/150‑1 (MQ=255)
              gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGATTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGc                                                                                                                                    <  1:2690298/150‑1 (MQ=255)
                 tATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTAc                                                                                                                                 <  1:234013/150‑1 (MQ=255)
                     tCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATg                                                                                                                             <  1:3597733/150‑1 (MQ=255)
                      ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGt                                                                                                                            <  1:3001605/150‑1 (MQ=255)
                      ccATACGTGATTATCCACATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTTTACGCTTAAAACACTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGt                                                                                                                            <  1:2091734/150‑1 (MQ=255)
                       cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGtt                                                                                                                           <  2:4191919/150‑1 (MQ=255)
                       cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGtt                                                                                                                           <  1:3422667/150‑1 (MQ=255)
                          aCGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAACAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCGGCGGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATc                                                                                                                        <  2:4512226/150‑1 (MQ=255)
                          aCGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATc                                                                                                                        <  2:225204/150‑1 (MQ=255)
                           cGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCa                                                                                                                       <  1:4330529/150‑1 (MQ=255)
                            gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAg                                                                                                                      <  1:1710128/150‑1 (MQ=255)
                            gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAg                                                                                                                      <  2:1987709/150‑1 (MQ=255)
                               aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTCTCAgggt                                                                                                                   >  2:1356263/1‑150 (MQ=255)
                                  ctCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACTCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTGTGTTCTATATAAATCGAGATTATAACTGTCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGa                                                                                                                >  1:3862918/2‑150 (MQ=255)
                                  aTCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGa                                                                                                                >  1:4604121/1‑150 (MQ=255)
                                  aTCCTCACGCGAACGGTCGCAAGAGTACACCAAACAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGCGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGa                                                                                                                <  2:2715442/150‑1 (MQ=255)
                                     cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGAGGTTATCAGTGTTGATGc                                                                                                             >  2:323188/1‑150 (MQ=255)
                                      tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGAGGTTATCAGTGTTGATGcc                                                                                                            >  2:4138422/1‑150 (MQ=255)
                                         tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAATGGGCCCGCAGGCCCGTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCCGc                                                                                                         >  1:3069850/1‑150 (MQ=255)
                                         tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGc                                                                                                         <  1:2281337/150‑1 (MQ=255)
                                          gCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCg                                                                                                        <  2:2869195/150‑1 (MQ=255)
                                          gCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCg                                                                                                        <  2:761028/150‑1 (MQ=255)
                                            gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTc                                                                                                      >  2:2062371/1‑150 (MQ=255)
                                             aaCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCCGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCa                                                                                                     >  2:240714/1‑150 (MQ=255)
                                              aCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCAc                                                                                                    >  2:4609461/1‑150 (MQ=255)
                                                ggTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                                                                  <  2:1685584/150‑1 (MQ=255)
                                                    gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAACAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGa                                                                                              <  2:4285010/150‑1 (MQ=255)
                                                     cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                                                             >  2:1465048/1‑150 (MQ=255)
                                                     cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                                                             >  2:2184065/1‑150 (MQ=255)
                                                        gagTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTg                                                                                          >  2:120166/1‑150 (MQ=255)
                                                           tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                                                       <  1:3202410/150‑1 (MQ=255)
                                                              acCAAAAACCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAACTCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtc                                                                                    <  1:2197736/150‑1 (MQ=255)
                                                              acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtc                                                                                    >  1:1424998/1‑150 (MQ=255)
                                                                cAAAAAACTGTGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgt                                                                                  <  1:2167508/150‑1 (MQ=255)
                                                                cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgt                                                                                  <  1:3324795/150‑1 (MQ=255)
                                                                    aaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcag                                                                              >  1:1364707/1‑150 (MQ=255)
                                                                     aaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCCGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagc                                                                             >  2:166078/1‑150 (MQ=255)
                                                                      aCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                                                            >  2:215257/1‑150 (MQ=255)
                                                                       cTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagcag                                                                           >  1:2414141/1‑150 (MQ=255)
                                                                             ttGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAt                                                                     >  2:792330/1‑150 (MQ=255)
                                                                               gTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAttt                                                                   >  2:3404356/1‑150 (MQ=255)
                                                                                 aCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTg                                                                 >  1:3340453/1‑150 (MQ=255)
                                                                                       aaaaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCa                                                           >  2:978645/1‑150 (MQ=255)
                                                                                         aaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAAc                                                         >  1:613105/1‑150 (MQ=255)
                                                                                             aTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAg                                                     >  2:851051/1‑150 (MQ=255)
                                                                                              tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg                                                    <  2:55962/150‑1 (MQ=255)
                                                                                               gAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt                                                   >  2:2309189/1‑150 (MQ=255)
                                                                                                  aaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa                                                >  2:2986669/1‑150 (MQ=255)
                                                                                                    aaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACa                                              >  2:444167/1‑150 (MQ=255)
                                                                                                       ggCCCTCAGGCCCATTGTTCGATATAAATCTAGATAAGAACAGGAAAACCAGACCCAGAGCTACGATGTTATCAGTGTAGATGCCAGCGTCACGAGTGACCTGGTTTTCGTCCAGCAGTTTGATATTGTAGTCAACGTAGGTGGACATGt                                           >  2:699558/1‑150 (MQ=255)
                                                                                                         cccGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttt                                         <  2:3733131/150‑1 (MQ=255)
                                                                                                          ccGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGtttt                                        <  1:4292759/150‑1 (MQ=255)
                                                                                                          ccGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGtttt                                        <  1:1604734/150‑1 (MQ=255)
                                                                                                          ccGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGtttt                                        <  2:50533/150‑1 (MQ=255)
                                                                                                           cGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                                       <  1:1099268/150‑1 (MQ=255)
                                                                                                           cGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                                       >  2:311352/1‑150 (MQ=255)
                                                                                                           cGCAGGCCATTTGTTCTATATAAATCGAGATTAGACCTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACTTGttttt                                       <  1:2585232/150‑1 (MQ=39)
                                                                                                              aGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgtt                                    >  2:4167267/1‑150 (MQ=255)
                                                                                                              aGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgtt                                    >  1:3153665/1‑150 (MQ=255)
                                                                                                               ggCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                   <  2:2075738/150‑1 (MQ=255)
                                                                                                                 cccTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                 >  2:3531596/1‑150 (MQ=255)
                                                                                                                  ccTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                >  1:2504232/1‑150 (MQ=255)
                                                                                                                   cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt                               <  2:3305297/150‑1 (MQ=255)
                                                                                                                   cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt                               <  1:3895796/150‑1 (MQ=255)
                                                                                                                   cTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt                               <  1:4057322/150‑1 (MQ=255)
                                                                                                                     ttGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                             >  2:3552001/1‑150 (MQ=255)
                                                                                                                      tGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                            >  1:1196640/1‑150 (MQ=255)
                                                                                                                      tGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGACATCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                            >  1:1196208/1‑150 (MQ=255)
                                                                                                                       gTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAgtagta                           >  2:291799/1‑150 (MQ=255)
                                                                                                                        ttCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCGACGTAGGTGGACATGTTTTTGTTGAAGtagtag                          <  2:959195/150‑1 (MQ=255)
                                                                                                                         tCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                         <  1:962679/150‑1 (MQ=255)
                                                                                                                           gATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggta                       <  2:2364922/150‑1 (MQ=255)
                                                                                                                              ataaaTCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCa                    <  1:2517573/147‑1 (MQ=255)
                                                                                                                               taaaTCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAc                   >  2:361647/3‑150 (MQ=255)
                                                                                                                                 aaTCTAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACca                 >  1:968206/1‑150 (MQ=255)
                                                                                                                                  aTCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaa                <  2:1829606/150‑1 (MQ=255)
                                                                                                                                   tCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaac               <  2:3775058/150‑1 (MQ=255)
                                                                                                                                     gagaTTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacat             <  1:1020277/150‑1 (MQ=255)
                                                                                                                                     gagaTTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacat             >  1:3630433/1‑150 (MQ=255)
                                                                                                                                     gagaTTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacat             >  2:2741839/1‑150 (MQ=255)
                                                                                                                                       gaTTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacatca           >  1:1118188/1‑150 (MQ=255)
                                                                                                                                           aGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacatcaacat       <  2:4267102/150‑1 (MQ=255)
                                                                                                                                           aGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacatcaacat       <  1:2934982/150‑1 (MQ=255)
                                                                                                                                           aGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacatcaacat       <  1:4088377/150‑1 (MQ=255)
                                                                                                                                             aaCTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACatat     >  1:636841/1‑150 (MQ=255)
                                                                                                                                             aaCTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACatat     >  2:1647515/1‑150 (MQ=255)
                                                                                                                                               cTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt   >  1:1769204/1‑150 (MQ=255)
                                                                                                                                                tGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTc  >  1:3279729/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTC  >  CP009273/2304988‑2305281

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: