Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,041,662 G→T 11.2% F580L (TTC→TTA gfcD ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,041,6620GT11.2% 804.8 / 49.6 222F580L (TTC→TTAgfcDputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (98/99);  new base T (11/14);  total (109/113)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 226 total aligned reads displayed.

CTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTCCAG  >  CP009273/1041513‑1041811
                                                                                                                                                     |                                                                                                                                                     
cTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACg                                                                                                                                                       <  1:210529/150‑1 (MQ=255)
 tGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGa                                                                                                                                                      <  1:3357181/150‑1 (MQ=255)
  gATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGaa                                                                                                                                                     >  1:527741/1‑150 (MQ=255)
   aTTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGaaa                                                                                                                                                    >  2:2337114/1‑150 (MQ=255)
   aTTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGaaa                                                                                                                                                    >  2:346615/1‑150 (MQ=255)
   aTTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGGTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGaaa                                                                                                                                                    >  1:296865/1‑150 (MQ=255)
    ttGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAAc                                                                                                                                                   >  2:741352/1‑150 (MQ=255)
    ttGCCGCACATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACTTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAAc                                                                                                                                                   <  2:1346056/150‑1 (MQ=255)
     tGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAAcc                                                                                                                                                  >  1:3078295/1‑150 (MQ=255)
        cgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGt                                                                                                                                               >  2:3597051/1‑150 (MQ=255)
           ccATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTgg                                                                                                                                            >  2:2810288/1‑150 (MQ=255)
            cATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTgga                                                                                                                                           <  1:4231217/150‑2 (MQ=255)
              tACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGt                                                                                                                                         >  2:4327587/1‑150 (MQ=255)
                 cccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGa                                                                                                                                      <  2:854149/150‑1 (MQ=255)
                    tACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGaata                                                                                                                                   <  2:1787073/150‑1 (MQ=255)
                     aCCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGaataa                                                                                                                                  >  2:4376114/1‑150 (MQ=255)
                      ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGaataat                                                                                                                                 <  1:2694744/150‑1 (MQ=255)
                       cGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATc                                                                                                                                <  2:3505097/150‑1 (MQ=255)
                           aCGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTg                                                                                                                            >  2:1831622/1‑150 (MQ=255)
                               cgcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAAc                                                                                                                        >  1:1008710/1‑150 (MQ=255)
                               cgcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAAc                                                                                                                        >  1:901096/1‑150 (MQ=255)
                               cgcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAAc                                                                                                                        <  2:109685/150‑1 (MQ=255)
                                gcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACg                                                                                                                       <  2:1590620/150‑1 (MQ=255)
                                  tGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGca                                                                                                                     <  2:264159/150‑1 (MQ=255)
                                       aaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCAtg                                                                                                                >  2:496350/1‑150 (MQ=255)
                                        aGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCAtgt                                                                                                               <  2:42802/150‑1 (MQ=255)
                                          cGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCAtgttg                                                                                                             <  1:4466550/150‑1 (MQ=255)
                                          cGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCAtgttg                                                                                                             <  2:3342480/150‑1 (MQ=255)
                                                  gACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAGGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTc                                                                                                     <  2:2462185/150‑1 (MQ=255)
                                                  gACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTc                                                                                                     >  1:1994456/1‑150 (MQ=255)
                                                  gACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTc                                                                                                     <  2:2462563/150‑1 (MQ=255)
                                                    cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCAc                                                                                                   <  1:2652922/150‑1 (MQ=255)
                                                    cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTCGTCCCAGTCAc                                                                                                   >  2:2582828/1‑150 (MQ=255)
                                                     gTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACg                                                                                                  >  2:4118371/1‑150 (MQ=255)
                                                       cGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGtt                                                                                                <  2:1687531/150‑1 (MQ=255)
                                                       cGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGtt                                                                                                >  2:3102230/1‑150 (MQ=255)
                                                        gATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGAGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTg                                                                                               >  1:1083004/1‑150 (MQ=255)
                                                        gATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTg                                                                                               >  2:3482854/1‑150 (MQ=255)
                                                        gATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTg                                                                                               >  1:3107192/1‑150 (MQ=255)
                                                        gATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTg                                                                                               <  1:1514594/150‑1 (MQ=255)
                                                           cGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCtt                                                                                            <  1:1706282/150‑1 (MQ=255)
                                                            gTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTc                                                                                           <  1:1227663/150‑1 (MQ=255)
                                                             ttGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCa                                                                                          <  2:1739074/150‑1 (MQ=255)
                                                              tGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCAc                                                                                         <  1:2957458/150‑1 (MQ=255)
                                                              tGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCAc                                                                                         <  1:3170046/150‑1 (MQ=255)
                                                               gcgcCTTAATATTTTGCCAGATATAGACCAACGCTAATGTTAATAAGTACGCCAATGATCTACGGCGGGTTCCAGTAATCCGAAACGAAACCTGATGTCGTGTCATAATCGCTGAAACGCAACATGTTGTCCCAGTCACGTTGCTTCACg                                                                                        >  2:661910/1‑150 (MQ=255)
                                                                cgcCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGt                                                                                       >  2:4546037/1‑150 (MQ=255)
                                                                cgcCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGt                                                                                       <  2:1149296/150‑1 (MQ=255)
                                                                      tATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTaa                                                                                 >  1:4434878/1‑150 (MQ=255)
                                                                        tCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACg                                                                               >  1:2901464/1‑150 (MQ=255)
                                                                        tCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACg                                                                               <  1:915884/150‑1 (MQ=255)
                                                                           tttGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGACTAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTcc                                                                            >  1:3308245/1‑150 (MQ=255)
                                                                           tttGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTcc                                                                            <  1:2764580/150‑1 (MQ=255)
                                                                              gCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCcaca                                                                         >  2:3341093/1‑150 (MQ=255)
                                                                                cAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACAcc                                                                       >  1:1959576/1‑150 (MQ=255)
                                                                                 aGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACAccc                                                                      >  1:3443115/1‑150 (MQ=255)
                                                                                       tGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGcc                                                                >  1:517851/1‑150 (MQ=255)
                                                                                         cccAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCca                                                              <  1:243710/150‑1 (MQ=255)
                                                                                          ccAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCcag                                                             <  1:4250841/150‑1 (MQ=255)
                                                                                           cAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCcagc                                                            <  1:1907376/150‑1 (MQ=255)
                                                                                           cAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCcagc                                                            <  2:739966/150‑1 (MQ=255)
                                                                                             aCGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCcagcag                                                          <  1:1264420/150‑1 (MQ=255)
                                                                                               gCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGc                                                        >  1:2800405/1‑150 (MQ=255)
                                                                                                 tAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCAt                                                      >  2:1955441/1‑150 (MQ=255)
                                                                                                     tttCATCAGTACGCCTTTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAg                                                  <  2:74407/150‑1 (MQ=255)
                                                                                                      ttCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGc                                                 >  1:1925835/1‑150 (MQ=255)
                                                                                                              tACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATa                                         <  2:2929096/150‑1 (MQ=255)
                                                                                                              tACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATa                                         <  2:2929533/150‑1 (MQ=255)
                                                                                                              tACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATa                                         <  2:3349735/150‑1 (MQ=255)
                                                                                                               aCGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATaa                                        >  1:3112710/1‑150 (MQ=255)
                                                                                                               aCGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATaa                                        >  2:3773401/1‑150 (MQ=255)
                                                                                                                cGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATaaa                                       >  2:1124811/1‑150 (MQ=255)
                                                                                                                  ccATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagc                                     <  2:2670178/150‑1 (MQ=255)
                                                                                                                   cATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagca                                    >  1:681799/1‑150 (MQ=255)
                                                                                                                   cATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagca                                    >  1:2406405/1‑150 (MQ=255)
                                                                                                                   cATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagca                                    <  2:517857/150‑1 (MQ=255)
                                                                                                                    aTTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagcag                                   >  1:3862074/1‑150 (MQ=255)
                                                                                                                     ttGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagcagc                                  >  2:909731/1‑150 (MQ=255)
                                                                                                                      tGAGCGTCTGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCAAGTTGCTTCACGTAGTTAACTTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCt                                 >  2:2154402/1‑150 (MQ=255)
                                                                                                                       gAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTc                                >  2:2061723/1‑150 (MQ=255)
                                                                                                                        aGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCg                               <  2:527790/150‑1 (MQ=255)
                                                                                                                             cggcggGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAAcc                          <  1:3706929/150‑1 (MQ=255)
                                                                                                                             cggcggGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAAcc                          <  1:329622/150‑1 (MQ=255)
                                                                                                                               gcggGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGa                        <  2:1877057/150‑1 (MQ=255)
                                                                                                                                cggGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGac                       <  1:2254153/150‑1 (MQ=255)
                                                                                                                                 gggTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCGGCTCGGAACCGaca                      <  1:3228853/150‑1 (MQ=255)
                                                                                                                                 gggATCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGaca                      >  2:2647076/1‑150 (MQ=255)
                                                                                                                                   gTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACAcc                    <  1:4206028/150‑1 (MQ=255)
                                                                                                                                    ttCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGACCCGACACCg                   <  2:512526/150‑1 (MQ=255)
                                                                                                                                    ttCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCg                   >  1:2503908/1‑150 (MQ=255)
                                                                                                                                      ccAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGc                 >  1:3423556/1‑150 (MQ=255)
                                                                                                                                       cAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCg                <  2:2761317/150‑1 (MQ=255)
                                                                                                                                       cAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCg                <  1:983435/150‑1 (MQ=255)
                                                                                                                                       cAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCg                <  1:983398/150‑1 (MQ=255)
                                                                                                                                           aaGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACa            <  2:3700964/150‑1 (MQ=255)
                                                                                                                                           aaGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACa            >  2:811079/1‑150 (MQ=255)
                                                                                                                                            aGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACAt           >  2:2751218/1‑150 (MQ=255)
                                                                                                                                                  cACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTc     >  2:3929380/1‑150 (MQ=255)
                                                                                                                                                   aCGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTcc    <  2:164019/150‑1 (MQ=255)
                                                                                                                                                     gAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTCCAg  >  1:4115181/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTCCAG  >  CP009273/1041513‑1041811

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