Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 288,249 T→A 100% T137T (ACA→ACT)  yagK ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273288,2490TA99.2% 773.9 / ‑2.8 234T137T (ACA→ACT) yagKCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (105/127);  total (105/129)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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 236 total aligned reads displayed.

GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACG  >  CP009273/288104‑288392
                                                                                                                                                 |                                                                                                                                               
   aGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGat                                                                                                                                            >  1:1027683/1‑148 (MQ=255)
   aGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGat                                                                                                                                            >  1:1740838/1‑148 (MQ=255)
      ccAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatcat                                                                                                                                         >  2:852428/1‑148 (MQ=255)
      ccAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatcat                                                                                                                                         <  2:189887/148‑1 (MQ=255)
           cGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCctc                                                                                                                                    <  2:1740838/148‑1 (MQ=255)
           cGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCctc                                                                                                                                    <  1:1561820/148‑1 (MQ=255)
              gTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTaa                                                                                                                                 >  1:2809189/1‑148 (MQ=255)
               ttAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAg                                                                                                                                <  2:1839935/148‑1 (MQ=255)
                tAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAgg                                                                                                                               <  1:1043080/148‑1 (MQ=255)
                   cAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTa                                                                                                                            <  2:1311752/148‑1 (MQ=255)
                    aGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTAt                                                                                                                           <  1:2226761/148‑1 (MQ=255)
                        ccTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                       <  2:1142231/148‑1 (MQ=255)
                         cTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTcc                                                                                                                      <  2:3354145/148‑1 (MQ=255)
                          tGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCt                                                                                                                     <  1:617793/148‑1 (MQ=255)
                            gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGa                                                                                                                   >  2:2954789/1‑148 (MQ=255)
                            gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGa                                                                                                                   <  2:645393/148‑1 (MQ=255)
                              aGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATc                                                                                                                 >  1:890282/1‑148 (MQ=255)
                                 tCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTa                                                                                                              <  2:1676457/148‑1 (MQ=255)
                                     ctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATcc                                                                                                          >  2:3132172/1‑148 (MQ=255)
                                      tgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATccc                                                                                                         <  1:1504103/148‑1 (MQ=255)
                                      gggAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGGTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATccc                                                                                                         <  2:841404/147‑1 (MQ=255)
                                       ggAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATcccc                                                                                                        >  2:2268391/1‑148 (MQ=255)
                                       ggAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATcccc                                                                                                        >  1:1426243/1‑148 (MQ=255)
                                        gAATTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCg                                                                                                       >  1:1907002/1‑148 (MQ=255)
                                        gAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCg                                                                                                       >  2:2019416/1‑148 (MQ=255)
                                           gTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAgg                                                                                                    >  2:2830064/1‑148 (MQ=255)
                                             cAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtg                                                                                                  <  2:1027669/148‑1 (MQ=255)
                                              aGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtgg                                                                                                 <  2:2740315/148‑1 (MQ=255)
                                                tACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgta                                                                                               >  2:769651/1‑148 (MQ=255)
                                                    aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                           <  2:3173222/148‑1 (MQ=255)
                                                    aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                           <  2:142855/148‑1 (MQ=255)
                                                    aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTaa                                                                                           <  2:2238611/148‑1 (MQ=255)
                                                        gtgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGt                                                                                       >  1:801179/1‑148 (MQ=255)
                                                        gtgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGt                                                                                       <  1:1068056/148‑1 (MQ=255)
                                                         tgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTc                                                                                      <  2:1582941/148‑1 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                     <  1:2951838/148‑1 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                     <  2:3505331/148‑1 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                     <  1:668851/148‑1 (MQ=255)
                                                          gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                     <  1:2364630/148‑1 (MQ=255)
                                                             tAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttg                                                                                  <  1:519236/148‑1 (MQ=255)
                                                               gCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgtt                                                                                <  1:2512343/148‑1 (MQ=255)
                                                                 acaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGa                                                                              <  2:2533037/148‑1 (MQ=255)
                                                                 acaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGa                                                                              <  2:3293049/148‑1 (MQ=255)
                                                                   aTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaac                                                                            >  2:3530442/1‑148 (MQ=255)
                                                                     aCTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaa                                                                          <  2:2822638/148‑1 (MQ=255)
                                                                       ttACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAacaaca                                                                        >  2:2474460/1‑148 (MQ=255)
                                                                        tACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAg                                                                       >  2:1037058/1‑148 (MQ=255)
                                                                        tACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAg                                                                       <  1:500509/148‑1 (MQ=255)
                                                                         aCAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAgg                                                                      <  1:2311671/148‑1 (MQ=255)
                                                                          cAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGc                                                                     <  2:358962/148‑1 (MQ=255)
                                                                           aGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCa                                                                    <  2:3411943/148‑1 (MQ=255)
                                                                            gTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaa                                                                   <  1:932307/148‑1 (MQ=255)
                                                                              tttCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAt                                                                 >  2:102772/1‑148 (MQ=255)
                                                                                  aGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGa                                                             <  1:1011144/148‑1 (MQ=255)
                                                                                   ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                            <  1:164809/148‑1 (MQ=255)
                                                                                    gAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATa                                                           >  2:3041961/1‑148 (MQ=255)
                                                                                     aaaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAg                                                          <  2:924867/148‑1 (MQ=255)
                                                                                      aaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGt                                                         <  1:2130580/148‑1 (MQ=255)
                                                                                       aaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTg                                                        >  1:1491283/1‑148 (MQ=255)
                                                                                        aTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                       <  1:2753825/148‑1 (MQ=255)
                                                                                        aTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                       >  2:639739/1‑148 (MQ=255)
                                                                                        aTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                       >  2:3108283/1‑148 (MQ=255)
                                                                                         tGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                      >  2:3183557/1‑148 (MQ=255)
                                                                                         tGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                      >  1:1666053/1‑148 (MQ=255)
                                                                                          ggACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCa                                                     <  2:263778/148‑1 (MQ=255)
                                                                                              cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGATCAACAGGCAAATGTGATAGTGGCATTTa                                                 <  2:500380/148‑1 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                <  1:2265202/148‑1 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                >  2:1218211/1‑148 (MQ=255)
                                                                                                gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAcc                                               <  1:3391456/148‑1 (MQ=255)
                                                                                                  ccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                             <  1:405699/148‑1 (MQ=255)
                                                                                                   cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGa                                            >  1:734157/1‑148 (MQ=255)
                                                                                                   cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGa                                            >  2:598297/1‑148 (MQ=255)
                                                                                                   cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGa                                            <  1:3533510/148‑1 (MQ=255)
                                                                                                    agGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                           >  2:2385992/1‑148 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                          <  2:890282/148‑1 (MQ=255)
                                                                                                      ggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                         >  1:2393249/1‑148 (MQ=255)
                                                                                                         gATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                      >  2:2089636/1‑148 (MQ=255)
                                                                                                           tCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAt                                    >  1:2522655/1‑148 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                              >  1:980892/1‑148 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                              <  1:1885131/148‑1 (MQ=255)
                                                                                                                  tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCttt                             <  1:3134351/148‑1 (MQ=255)
                                                                                                                  tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCttt                             <  1:3134128/148‑1 (MQ=255)
                                                                                                                  tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCttt                             <  1:2859697/148‑1 (MQ=255)
                                                                                                                      cAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGcc                         <  2:2858396/148‑1 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCa                       >  2:2776029/1‑148 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCa                       <  1:538189/148‑1 (MQ=255)
                                                                                                                         tctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                      >  2:1482687/1‑148 (MQ=255)
                                                                                                                           tcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCatat                    <  1:1245819/148‑1 (MQ=255)
                                                                                                                             aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataa                  <  2:398903/148‑1 (MQ=255)
                                                                                                                               ggCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataata                >  1:2610830/1‑148 (MQ=255)
                                                                                                                                  gcTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat             <  2:2605332/148‑1 (MQ=255)
                                                                                                                                  gcTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat             <  1:557059/148‑1 (MQ=255)
                                                                                                                                   cTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaata            <  2:1909753/148‑1 (MQ=255)
                                                                                                                                    tGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAataaataa           <  1:3093322/148‑1 (MQ=255)
                                                                                                                                     gTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAAtaaataaa          <  2:2270460/148‑1 (MQ=255)
                                                                                                                                      tACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAg         >  2:2053247/1‑148 (MQ=255)
                                                                                                                                          aGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCa     <  1:289271/148‑1 (MQ=255)
                                                                                                                                          aGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCa     >  2:1601804/1‑148 (MQ=255)
                                                                                                                                            gCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAAc   <  2:1308725/148‑1 (MQ=255)
                                                                                                                                             cTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACg  >  2:2448078/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACG  >  CP009273/288104‑288392

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: ‑