Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,894,388 C→A 94.3% V146V (GTC→GTA yqcE → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,894,3880CA94.3% 605.4 / 26.3 212V146V (GTC→GTAyqcEMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base C (3/9);  new base A (94/106);  total (97/115)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 212 total aligned reads displayed.

GTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGAT  >  CP009273/2894244‑2894532
                                                                                                                                                |                                                                                                                                                
 tgGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGgtg                                                                                                                                              >  2:1531784/1‑148 (MQ=255)
  gggTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGgtgt                                                                                                                                             >  1:1957786/1‑148 (MQ=255)
   ggTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTa                                                                                                                                            >  1:553307/1‑148 (MQ=255)
      aaTGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATg                                                                                                                                         <  2:1161414/148‑1 (MQ=255)
        tGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTc                                                                                                                                       >  2:874563/1‑148 (MQ=255)
           tctGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCt                                                                                                                                    >  2:945007/1‑148 (MQ=255)
           tctGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCt                                                                                                                                    >  2:2535268/1‑148 (MQ=255)
               tATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCg                                                                                                                                >  1:3176502/1‑148 (MQ=255)
                aTTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                                                                                                               >  2:1069925/1‑148 (MQ=255)
                 ttCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGgt                                                                                                                              >  2:1504066/1‑148 (MQ=255)
                   cAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGgtgt                                                                                                                            <  1:2683218/148‑1 (MQ=255)
                     gATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGttt                                                                                                                          >  1:1262084/1‑148 (MQ=255)
                      aTCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTa                                                                                                                         <  1:1092504/148‑1 (MQ=255)
                      aTCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTa                                                                                                                         >  2:1790656/1‑148 (MQ=255)
                      aTCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTa                                                                                                                         >  1:156879/1‑148 (MQ=255)
                       tCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTAc                                                                                                                        >  1:53123/1‑148 (MQ=255)
                       tCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTAc                                                                                                                        >  2:1658747/1‑148 (MQ=255)
                       tCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTAc                                                                                                                        <  1:3441009/148‑1 (MQ=255)
                        cGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATAGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGATGGCGGTGTTTAcc                                                                                                                       >  2:466700/1‑148 (MQ=255)
                                 gATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTc                                                                                                              <  1:1845935/148‑1 (MQ=255)
                                 gATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTc                                                                                                              <  1:580400/148‑1 (MQ=255)
                                 gATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTc                                                                                                              <  2:1632734/148‑1 (MQ=255)
                                  aTAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATAATAGCGAGCAAGGGACAATTATGGGCTGGATTGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCt                                                                                                             >  1:696019/1‑148 (MQ=255)
                                   tAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtt                                                                                                            <  2:2083739/148‑1 (MQ=255)
                                   tAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtt                                                                                                            >  1:3305534/1‑148 (MQ=255)
                                   tAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtt                                                                                                            <  2:627961/148‑1 (MQ=255)
                                    aaCGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCttt                                                                                                           <  2:2140220/148‑1 (MQ=255)
                                     acgacgATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtttt                                                                                                          >  2:2346003/1‑148 (MQ=255)
                                     acgacgATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtttt                                                                                                          >  1:215236/1‑148 (MQ=255)
                                      cgacgaTTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTtc                                                                                                         <  2:2407783/148‑1 (MQ=255)
                                          gaTTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGc                                                                                                     <  2:3482790/148‑1 (MQ=255)
                                            ttttAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCtt                                                                                                   <  1:2906117/148‑1 (MQ=255)
                                            ttttAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCtt                                                                                                   <  2:80797/148‑1 (MQ=255)
                                              ttAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTg                                                                                                 <  1:889106/148‑1 (MQ=255)
                                                   gCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTATCGCTTTGCACCg                                                                                            >  1:973280/1‑148 (MQ=255)
                                                     tgtgGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGa                                                                                          >  2:3110965/1‑148 (MQ=255)
                                                          tCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACa                                                                                     <  2:2877862/148‑1 (MQ=255)
                                                              tgtCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCAc                                                                                 <  1:2842321/148‑1 (MQ=255)
                                                              tgtCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCAc                                                                                 <  1:1128848/148‑1 (MQ=255)
                                                                 cGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAg                                                                              >  1:1467308/1‑148 (MQ=255)
                                                                  gATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGc                                                                             <  2:1044562/148‑1 (MQ=255)
                                                                     tAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTg                                                                          <  2:1054248/148‑1 (MQ=255)
                                                                        aGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGaaa                                                                       <  1:361014/148‑1 (MQ=255)
                                                                        aGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGaaa                                                                       >  2:1376691/1‑148 (MQ=255)
                                                                         gCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGaaaa                                                                      >  1:1972102/1‑148 (MQ=255)
                                                                           tGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACg                                                                    <  1:136046/148‑1 (MQ=255)
                                                                            gCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACgg                                                                   <  1:3216151/148‑1 (MQ=255)
                                                                               tCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGtaa                                                                <  2:2916147/148‑3 (MQ=255)
                                                                                 gTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCAtt                                                              <  2:2856132/148‑1 (MQ=255)
                                                                                   tGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTat                                                            >  1:353352/1‑148 (MQ=255)
                                                                                   tGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTat                                                            >  2:2226758/1‑148 (MQ=255)
                                                                                   tGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTat                                                            >  1:2534983/1‑148 (MQ=255)
                                                                                    gCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatc                                                           <  1:1825346/148‑1 (MQ=255)
                                                                                      ttGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcat                                                         <  2:1704464/148‑1 (MQ=255)
                                                                                      ttGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcat                                                         >  1:1493401/1‑148 (MQ=255)
                                                                                       tGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcatc                                                        >  2:1983723/1‑148 (MQ=255)
                                                                                       tGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcatc                                                        >  2:1807374/1‑148 (MQ=255)
                                                                                        ggCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCt                                                       <  1:2591070/148‑1 (MQ=255)
                                                                                         gCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCta                                                      >  1:2526161/1‑148 (MQ=255)
                                                                                             tCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAgt                                                  >  2:2011268/1‑148 (MQ=255)
                                                                                                 aGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTgg                                              >  1:3043917/1‑148 (MQ=255)
                                                                                                 aGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTgg                                              <  2:2798883/148‑1 (MQ=255)
                                                                                                  gCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTTGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGt                                             <  2:1577348/148‑1 (MQ=255)
                                                                                                  gCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGt                                             <  2:14180/148‑1 (MQ=255)
                                                                                                     aGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTa                                          >  1:1625147/1‑148 (MQ=255)
                                                                                                      gCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTAc                                         <  1:20229/148‑1 (MQ=255)
                                                                                                       cAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACa                                        >  1:240837/1‑148 (MQ=255)
                                                                                                       cAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACa                                        <  1:1068506/148‑1 (MQ=255)
                                                                                                        aaGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACAt                                       <  1:95421/148‑1 (MQ=255)
                                                                                                         aGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATc                                      <  2:2954662/148‑1 (MQ=255)
                                                                                                         aGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATc                                      >  2:1272070/1‑148 (MQ=255)
                                                                                                         aGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATc                                      >  2:3021048/1‑148 (MQ=255)
                                                                                                         aGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGGAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATc                                      <  2:915556/148‑1 (MQ=255)
                                                                                                            gAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCttg                                   <  1:2991030/148‑1 (MQ=255)
                                                                                                            gAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCttg                                   <  1:3327710/148‑1 (MQ=255)
                                                                                                             aaaaTTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCAACGGATGACAGCACCAGCCTGAAAACTGTCATTATCATCTATAGTGTTTTGTACATCttgt                                  >  2:1094919/1‑148 (MQ=255)
                                                                                                               aaTTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCttgttg                                <  2:2528492/148‑1 (MQ=255)
                                                                                                                 ttATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTggg                              >  2:3274992/1‑148 (MQ=255)
                                                                                                                     gggCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGAtt                          <  2:1370808/148‑1 (MQ=255)
                                                                                                                     gggCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGAtt                          <  1:1531784/148‑1 (MQ=255)
                                                                                                                     gggCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGAtt                          <  1:2201178/148‑1 (MQ=255)
                                                                                                                     gggCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGAtt                          <  1:1678468/148‑1 (MQ=255)
                                                                                                                        cTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCtg                       <  1:1045898/148‑1 (MQ=255)
                                                                                                                        cTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCtg                       <  1:3487410/148‑1 (MQ=255)
                                                                                                                          ggatgtaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCtgtg                     <  2:2088472/148‑1 (MQ=255)
                                                                                                                          ggatgtaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCtgtg                     <  2:315857/148‑1 (MQ=255)
                                                                                                                           tatggaAGGGCTGCGCGGCGTATGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATATTGTGGTTTACATCTTGTTGGGTATTCTGTGc                    <  2:646187/147‑1 (MQ=255)
                                                                                                                           gatggaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGc                    <  2:773549/148‑1 (MQ=255)
                                                                                                                           gatggaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGc                    <  2:2994494/148‑1 (MQ=255)
                                                                                                                            atggaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCt                   >  2:1908191/1‑148 (MQ=255)
                                                                                                                               gaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCTGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGt                <  1:1796010/148‑1 (MQ=255)
                                                                                                                               gaAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGt                <  1:1504696/148‑1 (MQ=255)
                                                                                                                                  gggCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGtttt             <  2:3337253/148‑1 (MQ=255)
                                                                                                                                  gggCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGtttt             <  1:1376691/148‑1 (MQ=255)
                                                                                                                                  gggCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGtttt             <  1:1830919/148‑1 (MQ=255)
                                                                                                                                   ggCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGttttt            >  2:2765696/1‑148 (MQ=255)
                                                                                                                                       gcgcGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGtt        <  2:3107100/148‑1 (MQ=255)
                                                                                                                                           ggcgTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCg    <  2:2645283/148‑1 (MQ=255)
                                                                                                                                            gcgTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGa   <  1:1327054/148‑1 (MQ=255)
                                                                                                                                            gcgTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGa   <  1:1747322/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTATAGTGTGGTTTACATCTTGTTGGGGATTCTGTGCTGGTTTTTTGTTAGCGAT  >  CP009273/2894244‑2894532

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: