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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,184,4710AT73.5% 541.4 / 229.8 308I519I (ATT→ATAthiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base T (101/125);  total (142/166)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 314 total aligned reads displayed.

TTTCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCGGCGC  >  CP009273/4184323‑4184620
                                                                                                                                                    |                                                                                                                                                     
  tCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATc                                                                                                                                                    <  2:85347/150‑1 (MQ=255)
  tCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATc                                                                                                                                                    >  2:4019912/1‑150 (MQ=255)
   cATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCt                                                                                                                                                   >  2:392937/1‑150 (MQ=255)
     tCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTta                                                                                                                                                 <  2:2125945/150‑1 (MQ=255)
      cgggATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTtat                                                                                                                                                <  1:5192389/147‑1 (MQ=255)
      cGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTtat                                                                                                                                                <  2:4329470/150‑1 (MQ=255)
       gTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTtata                                                                                                                                               >  1:2467670/1‑150 (MQ=255)
        tGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAg                                                                                                                                              >  1:3919007/1‑150 (MQ=255)
         gATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAgg                                                                                                                                             <  2:248275/150‑1 (MQ=255)
          aTAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGt                                                                                                                                            >  1:1651604/1‑150 (MQ=255)
            aaGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGa                                                                                                                                          <  1:4121116/150‑1 (MQ=255)
            aaGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGa                                                                                                                                          >  2:1113986/1‑150 (MQ=255)
            aaGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGa                                                                                                                                          >  2:5133661/1‑150 (MQ=255)
             aGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGAt                                                                                                                                         <  2:4612502/150‑1 (MQ=255)
              gcgcGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATa                                                                                                                                        <  2:1350277/150‑1 (MQ=255)
               cgcgGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATaa                                                                                                                                       <  2:5203109/150‑1 (MQ=255)
                 cgggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGc                                                                                                                                     >  1:2347650/1‑150 (MQ=255)
                  gggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGcc                                                                                                                                    <  1:3023402/150‑1 (MQ=255)
                   ggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGccc                                                                                                                                   <  2:3369748/150‑1 (MQ=255)
                     cggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCt                                                                                                                                 <  1:275825/150‑1 (MQ=255)
                     cggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCt                                                                                                                                 <  1:1744925/150‑1 (MQ=255)
                      ggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTg                                                                                                                                <  2:138656/150‑1 (MQ=255)
                      ggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATGAGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTCTCTTATAGGTGATAAGCCCCTg                                                                                                                                >  2:2825486/1‑150 (MQ=255)
                       gTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGc                                                                                                                               >  2:4558587/1‑150 (MQ=255)
                        tAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTCTAGGTGATAAGCCCCTGCt                                                                                                                              >  2:809809/1‑150 (MQ=255)
                        tAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCt                                                                                                                              >  1:328458/1‑150 (MQ=255)
                          aaCGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTa                                                                                                                            <  1:951731/150‑1 (MQ=255)
                            cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAAc                                                                                                                          >  2:3710341/1‑150 (MQ=255)
                            cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAAc                                                                                                                          <  1:3107352/150‑1 (MQ=255)
                            cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAAc                                                                                                                          <  2:884098/150‑1 (MQ=255)
                              ggTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAt                                                                                                                        <  2:3228574/150‑1 (MQ=255)
                                 cGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAtctt                                                                                                                     <  1:3799819/150‑1 (MQ=255)
                                  gAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAtcttc                                                                                                                    <  1:456335/150‑1 (MQ=255)
                                     ggCCAGATTACACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCttt                                                                                                                 <  1:160837/150‑1 (MQ=255)
                                     ggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCttt                                                                                                                 >  1:1478463/1‑150 (MQ=255)
                                       ccAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAt                                                                                                               >  1:236908/1‑150 (MQ=255)
                                        cAGATTACACTGGTCTTCACAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                              <  2:617541/150‑1 (MQ=255)
                                        cAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGTTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                              >  1:3484607/1‑150 (MQ=255)
                                        cAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                              >  2:1707526/1‑150 (MQ=255)
                                         aGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTg                                                                                                             <  2:4483690/150‑1 (MQ=255)
                                         aGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTg                                                                                                             <  2:1881612/150‑1 (MQ=255)
                                           aTTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTggg                                                                                                           >  1:1015775/1‑150 (MQ=255)
                                           aTTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTggg                                                                                                           >  2:4979945/1‑150 (MQ=255)
                                            ttAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGc                                                                                                          <  2:641937/150‑1 (MQ=255)
                                             tAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCa                                                                                                         <  2:4722527/150‑1 (MQ=255)
                                              aaaCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAg                                                                                                        <  1:4019912/150‑1 (MQ=255)
                                              aaaCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAg                                                                                                        <  1:2054601/150‑1 (MQ=255)
                                                   ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCa                                                                                                   <  2:3097540/150‑1 (MQ=255)
                                                   ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCa                                                                                                   >  1:2528617/1‑150 (MQ=255)
                                                   ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCa                                                                                                   <  2:1298049/150‑1 (MQ=255)
                                                   ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCAGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCa                                                                                                   >  2:1864873/1‑150 (MQ=255)
                                                    gTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAg                                                                                                  >  2:745568/1‑150 (MQ=255)
                                                          ccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctc                                                                                            >  1:3765416/1‑150 (MQ=255)
                                                          ccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctc                                                                                            <  1:3056474/150‑1 (MQ=255)
                                                            aGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtt                                                                                          <  1:975794/150‑1 (MQ=255)
                                                             gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGTCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCttt                                                                                         <  1:493260/150‑1 (MQ=255)
                                                               gAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTg                                                                                       <  1:4327406/150‑1 (MQ=255)
                                                                aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTgg                                                                                      >  1:5125080/1‑150 (MQ=255)
                                                                 aaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAAATTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTAAACATCTTCTTTATTGTGCAGACCCAGATGCTCtttagtc                                                                                     >  1:4311952/1‑146 (MQ=255)
                                                                  aTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCg                                                                                    >  1:193763/1‑150 (MQ=255)
                                                                    tCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGtt                                                                                  <  2:3811923/150‑1 (MQ=255)
                                                                     cGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTa                                                                                 <  2:3427889/150‑1 (MQ=255)
                                                                      gAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTAc                                                                                >  2:5237311/1‑150 (MQ=255)
                                                                        aGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGt                                                                              >  1:4392286/1‑150 (MQ=255)
                                                                         gCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTa                                                                             <  1:1114144/150‑1 (MQ=255)
                                                                          cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                            <  1:4238483/150‑1 (MQ=255)
                                                                          cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                            >  2:3183995/1‑150 (MQ=255)
                                                                          cGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                            <  1:4476279/150‑1 (MQ=255)
                                                                           gggCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAAc                                                                           >  2:2184075/1‑150 (MQ=255)
                                                                              cTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACaga                                                                        >  2:2743215/1‑150 (MQ=255)
                                                                                  cGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCAt                                                                    <  2:4414582/150‑1 (MQ=255)
                                                                                   gACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATc                                                                   <  1:3355330/150‑1 (MQ=255)
                                                                                    aCATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcg                                                                  <  1:1551875/150‑1 (MQ=255)
                                                                                      aTGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcgcg                                                                >  2:4295494/1‑150 (MQ=255)
                                                                                          cGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGc                                                            >  2:797392/1‑150 (MQ=255)
                                                                                           gTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCGTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGcc                                                           >  2:1764685/1‑150 (MQ=255)
                                                                                           gTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGcc                                                           >  1:5319832/1‑150 (MQ=255)
                                                                                            ttATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCa                                                          <  2:4864453/150‑1 (MQ=255)
                                                                                                cgcgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAAcc                                                      <  2:2977610/150‑1 (MQ=255)
                                                                                                    aaTTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGcc                                                  <  2:2039843/150‑1 (MQ=255)
                                                                                                     aTTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCa                                                 <  1:4816718/150‑1 (MQ=255)
                                                                                                      tttGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCaa                                                <  1:3239666/150‑1 (MQ=255)
                                                                                                         gcgcgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatca                                             <  2:4898679/150‑1 (MQ=255)
                                                                                                           gcgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatcatc                                           <  2:2458670/150‑1 (MQ=255)
                                                                                                              cccGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcgcc                                        <  2:4392286/150‑1 (MQ=255)
                                                                                                              cccGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcgcc                                        <  2:3266831/150‑1 (MQ=255)
                                                                                                               ccGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcgccg                                       >  2:2267367/1‑150 (MQ=255)
                                                                                                                 ggATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCa                                     <  1:3100480/150‑1 (MQ=255)
                                                                                                                    tGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCa                                  <  2:1558464/150‑1 (MQ=255)
                                                                                                                       ccTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                               >  2:2582738/1‑150 (MQ=255)
                                                                                                                       ccTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                               >  1:2915311/1‑150 (MQ=255)
                                                                                                                           tCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCg                           <  2:3598281/150‑1 (MQ=255)
                                                                                                                                   tCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAg                   >  1:2663940/1‑150 (MQ=255)
                                                                                                                                   tCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAg                   >  2:1208786/1‑150 (MQ=255)
                                                                                                                                      gcggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTgg                <  1:224439/150‑1 (MQ=255)
                                                                                                                                       cggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGt               <  1:4750190/150‑1 (MQ=255)
                                                                                                                                       cggcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGt               <  2:183080/150‑1 (MQ=255)
                                                                                                                                        ggcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTc              >  1:807368/1‑150 (MQ=255)
                                                                                                                                         gcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCa             <  1:3688315/150‑1 (MQ=255)
                                                                                                                                           gTGGGCAGCTATCGTAGAGGTGATAAGCCCCTGCTTTACTTCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGTGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATa           <  1:2825486/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCGGCGC  >  CP009273/4184323‑4184620

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: