Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,733,753 C→T 100% E253K (GAG→AAG)  aroF ← 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,733,7530CT89.6% 718.9 / 58.5 231E253K (GAG→AAG) aroF3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base T (116/90);  total (128/103)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 237 total aligned reads displayed.

TCATCAGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGGGCGG  >  CP009273/2733606‑2733899
                                                                                                                                                   |                                                                                                                                                  
tcatcaGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCtttt                                                                                                                                                  >  2:3545969/1‑150 (MQ=255)
 catcaGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTtttt                                                                                                                                                 >  2:1625830/1‑150 (MQ=255)
 catcaGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGAGTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCttt                                                                                                                                                 <  2:2642990/150‑1 (MQ=255)
   tcaGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTtttttt                                                                                                                                               <  2:642144/150‑1 (MQ=255)
    caGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTc                                                                                                                                              >  1:1164967/1‑150 (MQ=255)
     aGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTca                                                                                                                                             <  2:1256205/150‑1 (MQ=255)
      gACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTcac                                                                                                                                            <  1:3103404/150‑1 (MQ=255)
      gACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTcac                                                                                                                                            <  1:3103885/150‑1 (MQ=255)
      gACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTcac                                                                                                                                            >  2:3278198/1‑150 (MQ=255)
       aCCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTcaca                                                                                                                                           >  1:1766839/1‑150 (MQ=255)
       aCCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTcaca                                                                                                                                           >  1:4020334/1‑150 (MQ=255)
       aCCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTcaca                                                                                                                                           >  1:4020183/1‑150 (MQ=255)
         cAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACAtt                                                                                                                                         >  1:3709006/1‑150 (MQ=255)
         cAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACAtt                                                                                                                                         >  1:878516/1‑150 (MQ=255)
           ataatTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTgc                                                                                                                                       >  2:744449/1‑150 (MQ=255)
                     gATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCACCATCCg                                                                                                                             <  2:3829655/150‑1 (MQ=255)
                     gATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCg                                                                                                                             >  2:3229559/1‑150 (MQ=255)
                     gATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCg                                                                                                                             >  2:3229833/1‑150 (MQ=255)
                       ttGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCa                                                                                                                           >  2:1283829/1‑150 (MQ=255)
                       ttGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCa                                                                                                                           <  2:4293863/150‑1 (MQ=255)
                       ttGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCa                                                                                                                           <  2:4295296/150‑1 (MQ=255)
                        tGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCTGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAg                                                                                                                          <  1:3834415/150‑1 (MQ=255)
                           cATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGc                                                                                                                       >  2:1346707/1‑150 (MQ=255)
                              cTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCtat                                                                                                                    >  1:3603018/1‑150 (MQ=255)
                              cTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCtat                                                                                                                    <  2:192675/150‑1 (MQ=255)
                              cTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCtat                                                                                                                    <  2:191491/150‑1 (MQ=255)
                                ttgatttgaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAg                                                                                                                  <  2:3499794/150‑1 (MQ=255)
                                  gatttgaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGtt                                                                                                                <  2:1084266/150‑1 (MQ=255)
                                   atttgaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTc                                                                                                               <  1:1708031/150‑1 (MQ=255)
                                    tttgaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCg                                                                                                              <  1:3096709/150‑1 (MQ=255)
                                      tgaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGgc                                                                                                            <  1:2940910/150‑1 (MQ=255)
                                       gaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGgcg                                                                                                           >  1:3133413/1‑150 (MQ=255)
                                       gaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGgcg                                                                                                           <  2:65230/150‑1 (MQ=255)
                                        aGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATCGTTCGgcgc                                                                                                          >  2:4143934/1‑150 (MQ=255)
                                        aGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGgcgc                                                                                                          >  2:1864162/1‑150 (MQ=255)
                                         gCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCt                                                                                                         <  2:2234614/150‑1 (MQ=255)
                                               aCGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTacca                                                                                                   <  1:3870923/150‑1 (MQ=255)
                                                   aTTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACcgc                                                                                               >  1:346651/1‑150 (MQ=255)
                                                         ccACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGa                                                                                         <  2:2048984/150‑1 (MQ=255)
                                                          cACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGAt                                                                                        >  1:328028/1‑150 (MQ=255)
                                                           aCCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATc                                                                                       <  1:2642499/150‑1 (MQ=255)
                                                            ccGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATca                                                                                      >  1:2273427/1‑150 (MQ=255)
                                                            ccGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATca                                                                                      >  1:1094481/1‑150 (MQ=255)
                                                            ccGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATca                                                                                      >  1:2272494/1‑150 (MQ=255)
                                                               cAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACAt                                                                                   <  1:2125033/150‑1 (MQ=255)
                                                                aGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATg                                                                                  >  1:1169540/1‑150 (MQ=255)
                                                                 ggCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATgg                                                                                 >  1:4493203/1‑150 (MQ=255)
                                                                       cGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTcc                                                                           <  1:541792/150‑1 (MQ=255)
                                                                         gCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCgg                                                                         <  2:1963460/150‑1 (MQ=255)
                                                                           gATAATCTTTATTGCAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGAt                                                                       <  2:2706007/150‑1 (MQ=255)
                                                                                tCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTcccc                                                                  >  1:1460186/1‑150 (MQ=255)
                                                                                   ttATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGa                                                               >  1:984440/1‑150 (MQ=255)
                                                                                   ttATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGa                                                               >  1:3261512/1‑150 (MQ=255)
                                                                                   ttATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGa                                                               >  1:3261536/1‑150 (MQ=255)
                                                                                    tATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAg                                                              >  1:772501/1‑150 (MQ=255)
                                                                                                cGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacg                                                  <  2:4470036/150‑1 (MQ=255)
                                                                                                 gTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgc                                                 >  1:3668028/1‑150 (MQ=255)
                                                                                                  tGGCTGCAATCTACCATCAGAGACGGGAGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATATTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgca                                                <  1:3640390/150‑1 (MQ=255)
                                                                                                   ggCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgcaa                                               >  1:2673956/1‑150 (MQ=255)
                                                                                                   ggCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgcaa                                               >  2:12682/1‑150 (MQ=255)
                                                                                                   ggCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgcaa                                               >  2:3414464/1‑150 (MQ=255)
                                                                                                   ggCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTGGTAgcaacgcaa                                               >  2:2319678/1‑150 (MQ=255)
                                                                                                     cTGCAATCTACCATCAGAGACGGGCGCAGTCCCTCCTGTTCCATCAATTTTTCACATTGCGCAACATCCTCAGGGCTATAGATCGGCGCTTTACCACCGCGCAGGATCACAAGGCCGACCGGATTCCCCAGAGTTTGTAGCAACGCAAcc                                             >  1:4316535/1‑150 (MQ=255)
                                                                                                      tGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcct                                            <  1:2956021/150‑1 (MQ=255)
                                                                                                      tGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcct                                            <  1:1573830/150‑1 (MQ=255)
                                                                                                      tGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcct                                            <  1:4154299/150‑1 (MQ=255)
                                                                                                      tGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCCCCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcct                                            <  1:4473535/150‑1 (MQ=255)
                                                                                                        cAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgc                                          <  2:2569532/150‑1 (MQ=255)
                                                                                                          aTCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccc                                        <  1:2508661/150‑1 (MQ=255)
                                                                                                          aTCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccc                                        <  2:528833/150‑1 (MQ=255)
                                                                                                              aCCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACctgccctgcc                                    <  1:2402938/150‑1 (MQ=255)
                                                                                                                cATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCcctgcctg                                  >  2:412301/1‑150 (MQ=255)
                                                                                                                  tCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGt                                >  2:1207220/1‑150 (MQ=255)
                                                                                                                   cAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGtt                               >  1:3590070/1‑150 (MQ=255)
                                                                                                                     gagaCGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTaa                             <  1:3235660/150‑1 (MQ=255)
                                                                                                                      agaCGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAAt                            >  2:417782/1‑150 (MQ=255)
                                                                                                                      agaCGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAAt                            >  1:1144996/1‑150 (MQ=255)
                                                                                                                      agaCGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAAt                            <  1:568928/150‑1 (MQ=255)
                                                                                                                         cGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGTCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAAGGCAACCTGCCCTGCCTGGTTAATGcc                         <  2:2901168/150‑1 (MQ=255)
                                                                                                                         cGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGcc                         <  2:73658/150‑1 (MQ=255)
                                                                                                                         cGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGcc                         >  2:3847825/1‑150 (MQ=255)
                                                                                                                          gggCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCca                        >  2:4046589/1‑150 (MQ=255)
                                                                                                                          gggCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCca                        >  1:4800/1‑150 (MQ=255)
                                                                                                                          gggCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCca                        >  2:1718205/1‑150 (MQ=255)
                                                                                                                          gggCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCca                        >  2:1718471/1‑150 (MQ=255)
                                                                                                                          gggCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCca                        >  2:1795262/1‑150 (MQ=255)
                                                                                                                               cAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACaaa                   <  1:3342453/150‑1 (MQ=255)
                                                                                                                               cAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACaaa                   >  1:3170237/1‑150 (MQ=255)
                                                                                                                                aGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACaaaa                  >  2:1530748/1‑150 (MQ=255)
                                                                                                                                    ccGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGt              >  1:1883408/1‑150 (MQ=255)
                                                                                                                                    ccGCCTGTTCCATCTCTTTTTCACAGTGCGCACCATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGt              <  1:3524271/150‑1 (MQ=255)
                                                                                                                                      gCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGc            >  1:698958/1‑150 (MQ=255)
                                                                                                                                         tGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGc         >  2:3741627/1‑150 (MQ=255)
                                                                                                                                         tGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGc         >  2:3612276/1‑150 (MQ=255)
                                                                                                                                         tGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGc         >  2:606685/1‑150 (MQ=255)
                                                                                                                                          gTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCt        >  2:288025/1‑150 (MQ=255)
                                                                                                                                            tCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTgg      <  1:4421317/150‑1 (MQ=255)
                                                                                                                                            tCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTgg      >  1:135711/1‑150 (MQ=255)
                                                                                                                                             ccATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTggg     >  2:821224/1‑150 (MQ=255)
                                                                                                                                                tcttTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGggcgg  <  1:4289286/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCATCAGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGGGCGG  >  CP009273/2733606‑2733899

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: