Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,184,471 A→T 42.1% I519I (ATT→ATA thiC ← phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,184,4710AT42.1% 130.3 / 289.4 216I519I (ATT→ATAthiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (58/67);  new base T (28/63);  total (86/130)
Fisher's exact test for biased strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TTCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCGG  >  CP009273/4184324‑4184617
                                                                                                                                                   |                                                                                                                                                  
 tCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCaa                                                                                                                                                   >  2:323165/1‑148 (MQ=255)
 tCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCaa                                                                                                                                                   >  2:3387537/1‑148 (MQ=255)
   aTCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATc                                                                                                                                                 <  1:2472124/148‑1 (MQ=255)
    tCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAACCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCt                                                                                                                                                <  1:2809108/148‑1 (MQ=255)
    tCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCt                                                                                                                                                >  2:1184356/1‑148 (MQ=255)
    tCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCt                                                                                                                                                <  1:2050342/148‑1 (MQ=255)
       tGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTtat                                                                                                                                             >  1:2489635/1‑148 (MQ=255)
       tGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTtat                                                                                                                                             <  1:316663/148‑1 (MQ=255)
        gATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTtata                                                                                                                                            <  2:890954/148‑1 (MQ=255)
           aaGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGt                                                                                                                                         >  1:662206/1‑148 (MQ=255)
            aGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTg                                                                                                                                        <  1:918859/148‑1 (MQ=255)
             gcgcGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGa                                                                                                                                       <  2:71058/148‑1 (MQ=255)
              cgcgGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGAt                                                                                                                                      >  1:3433211/1‑148 (MQ=255)
                  ggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGc                                                                                                                                  <  2:3345870/148‑1 (MQ=255)
                       tAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTg                                                                                                                             <  2:3089788/148‑1 (MQ=255)
                        aaaCGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGc                                                                                                                            <  1:3303765/148‑1 (MQ=255)
                        aaaCGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGc                                                                                                                            >  2:2681892/1‑148 (MQ=255)
                           cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTa                                                                                                                         <  2:3594647/148‑1 (MQ=255)
                           cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTa                                                                                                                         <  2:3569755/148‑1 (MQ=255)
                               tCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTGATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAt                                                                                                                     <  2:2131404/148‑1 (MQ=255)
                                cGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAtc                                                                                                                    <  1:1460682/148‑1 (MQ=255)
                                   gggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtcttc                                                                                                                 >  2:1418542/1‑148 (MQ=255)
                                     gCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATcttctt                                                                                                               <  1:2688001/148‑1 (MQ=255)
                                      ccAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCttt                                                                                                              <  2:3655793/148‑1 (MQ=255)
                                      ccAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCttt                                                                                                              <  2:2122256/148‑1 (MQ=255)
                                      ccAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCttt                                                                                                              <  2:1289630/148‑1 (MQ=255)
                                        aGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAt                                                                                                            >  1:548380/1‑148 (MQ=255)
                                        aGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAt                                                                                                            <  2:1731014/148‑1 (MQ=255)
                                         gATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGTCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                           <  2:1543691/148‑1 (MQ=255)
                                         gATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                           <  2:2785734/148‑1 (MQ=255)
                                         gATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                           >  1:170674/1‑148 (MQ=255)
                                         gATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                           <  2:3170482/148‑1 (MQ=255)
                                          aTTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTg                                                                                                          <  1:221726/148‑1 (MQ=255)
                                          aTTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTg                                                                                                          <  2:450596/148‑1 (MQ=255)
                                          aTTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTg                                                                                                          >  1:1121784/1‑148 (MQ=255)
                                           ttAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTgg                                                                                                         >  1:1641675/1‑148 (MQ=255)
                                            tAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTggg                                                                                                        >  2:3369003/1‑148 (MQ=255)
                                             aaaCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGc                                                                                                       <  1:2449771/148‑1 (MQ=255)
                                             aaaCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGc                                                                                                       <  1:1053955/148‑1 (MQ=255)
                                               aCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAg                                                                                                     >  1:2652300/1‑148 (MQ=255)
                                                  ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAcc                                                                                                  >  2:1004978/1‑148 (MQ=255)
                                                  ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAcc                                                                                                  <  1:157811/148‑1 (MQ=255)
                                                   gTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAccc                                                                                                 >  1:2159459/1‑148 (MQ=255)
                                                   gTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAccc                                                                                                 <  1:512374/148‑1 (MQ=255)
                                                      ttCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGa                                                                                              <  2:3726794/148‑1 (MQ=255)
                                                       tCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGAt                                                                                             >  2:2141088/1‑148 (MQ=255)
                                                        cccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATg                                                                                            <  2:2967283/148‑1 (MQ=255)
                                                         ccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGc                                                                                           <  1:1627002/148‑1 (MQ=255)
                                                          cAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGct                                                                                          >  1:1074736/1‑148 (MQ=255)
                                                          cAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGct                                                                                          <  2:2693715/148‑1 (MQ=255)
                                                            gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                        <  1:3057705/148‑1 (MQ=255)
                                                            gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                        >  1:2853173/1‑148 (MQ=255)
                                                             cGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtt                                                                                       <  2:1843521/148‑1 (MQ=255)
                                                               aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                     <  2:1310321/148‑1 (MQ=255)
                                                               aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                     >  2:2348883/1‑148 (MQ=255)
                                                               aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                     >  2:1271389/1‑148 (MQ=255)
                                                               aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                     >  2:2316421/1‑148 (MQ=255)
                                                                aaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTg                                                                                    >  1:114711/1‑148 (MQ=255)
                                                                 aTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTgg                                                                                   >  1:2355441/1‑148 (MQ=255)
                                                                    cGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGt                                                                                <  1:1955231/148‑1 (MQ=255)
                                                                     gAAGCGGGCTTTCGACATGGCGGTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGtt                                                                               <  2:3606552/148‑1 (MQ=255)
                                                                        gCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACg                                                                            <  2:1148625/148‑1 (MQ=255)
                                                                        gCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACg                                                                            >  2:3492892/1‑148 (MQ=255)
                                                                         cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTTATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGt                                                                           <  1:1815325/148‑1 (MQ=255)
                                                                         cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGt                                                                           >  2:1321334/1‑148 (MQ=255)
                                                                         cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGt                                                                           >  1:1247027/1‑148 (MQ=255)
                                                                         cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGt                                                                           >  1:1247413/1‑148 (MQ=255)
                                                                           ggCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                         >  2:1385847/1‑148 (MQ=255)
                                                                           ggCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                         >  1:1002695/1‑148 (MQ=255)
                                                                            gCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAAc                                                                        <  2:739264/148‑1 (MQ=255)
                                                                            gCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAAc                                                                        <  2:2654732/148‑1 (MQ=255)
                                                                               ttCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACaga                                                                     <  1:3492892/148‑1 (MQ=255)
                                                                                     aTGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcg                                                               <  1:3007769/148‑1 (MQ=255)
                                                                                        gCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCgcgc                                                            <  1:3244037/148‑1 (MQ=255)
                                                                                               cgcgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCaaa                                                     >  1:564870/1‑148 (MQ=255)
                                                                                               cgcgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGACCCTGCTTAACATCTTCTTTATTGTGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCaaa                                                     >  1:1195499/1‑148 (MQ=255)
                                                                                                   aaTTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAg                                                 <  1:818519/148‑1 (MQ=255)
                                                                                                     tttGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGcc                                               <  1:3679042/148‑1 (MQ=255)
                                                                                                      ttGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCa                                              >  1:800488/1‑148 (MQ=255)
                                                                                                       tGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGATATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCaa                                             <  2:3174928/148‑1 (MQ=255)
                                                                                                             cccGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcg                                       <  2:548380/148‑1 (MQ=255)
                                                                                                                  aTGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCa                                  >  1:2039737/1‑148 (MQ=255)
                                                                                                                     cccTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCa                               <  2:170674/148‑1 (MQ=255)
                                                                                                                        tttCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                            >  1:1357225/1‑148 (MQ=255)
                                                                                                                        tttCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                            >  1:3035687/1‑148 (MQ=255)
                                                                                                                        tttCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                            <  2:2489635/148‑1 (MQ=255)
                                                                                                                        tttCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                            <  2:1170871/148‑1 (MQ=255)
                                                                                                                         ttCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATcc                           <  1:2113019/148‑1 (MQ=255)
                                                                                                                                    aGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAg                <  2:1040134/148‑1 (MQ=255)
                                                                                                                                      cggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCATAACATCTTCTTTATTGGGCAGACCCAGATGATCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGGGAAGTg              >  2:196299/1‑148 (MQ=255)
                                                                                                                                       ggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTgg             <  1:365866/148‑1 (MQ=255)
                                                                                                                                       ggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTgg             <  1:413798/148‑1 (MQ=255)
                                                                                                                                       ggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTgg             <  1:2337293/148‑1 (MQ=255)
                                                                                                                                       ggcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTgg             <  2:2355441/148‑1 (MQ=255)
                                                                                                                                        gcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGt            <  2:1247027/148‑1 (MQ=255)
                                                                                                                                        gcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGt            <  2:1247413/148‑1 (MQ=255)
                                                                                                                                         cgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTc           <  2:771326/148‑1 (MQ=255)
                                                                                                                                         cgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTc           <  2:28971/148‑1 (MQ=255)
                                                                                                                                          gTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCa          <  2:2400467/148‑1 (MQ=255)
                                                                                                                                               cagctATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAgtc     <  2:120474/148‑3 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCGG  >  CP009273/4184324‑4184617

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: