Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,144,697 C→A 7.3% A27S (GCA→TCA)  argE ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,144,6970CA 7.3% 1034.6 / 36.2 279A27S (GCA→TCA) argEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (108/150);  new base A (10/11);  total (118/161)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTT  >  CP009273/4144548‑4144846
                                                                                                                                                     |                                                                                                                                                     
 aacaaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGcc                                                                                                                                                      >  1:904670/1‑150 (MQ=255)
 aacaaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGcc                                                                                                                                                      >  1:1303473/1‑150 (MQ=255)
   caaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCctc                                                                                                                                                    <  2:1510893/150‑1 (MQ=255)
     aGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCtt                                                                                                                                                  <  1:3612568/150‑1 (MQ=255)
     aGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCtt                                                                                                                                                  <  1:3300718/150‑1 (MQ=255)
     aGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCtt                                                                                                                                                  <  1:1341112/150‑1 (MQ=255)
       ccgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTcc                                                                                                                                                <  1:1067828/150‑1 (MQ=255)
       ccgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTcc                                                                                                                                                <  2:3834087/150‑1 (MQ=255)
        cgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTTCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCg                                                                                                                                               >  1:824066/1‑150 (MQ=255)
        cgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCg                                                                                                                                               >  2:1854211/1‑150 (MQ=255)
        cgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCg                                                                                                                                               >  1:2647694/1‑150 (MQ=255)
          ccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTg                                                                                                                                             <  1:3260038/150‑1 (MQ=255)
          ccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTg                                                                                                                                             <  1:984405/150‑1 (MQ=255)
                  cTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTAt                                                                                                                                     <  1:2087297/150‑1 (MQ=255)
                     tCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGa                                                                                                                                  >  2:105496/1‑150 (MQ=255)
                           aCTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGgtgt                                                                                                                            <  2:679683/150‑1 (MQ=255)
                              tGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGc                                                                                                                         <  2:2564894/150‑1 (MQ=255)
                                 cAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAAt                                                                                                                      <  1:4006455/150‑1 (MQ=255)
                                  aGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATc                                                                                                                     <  1:3238379/150‑1 (MQ=255)
                                   gCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCa                                                                                                                    >  2:1148404/1‑150 (MQ=255)
                                       atTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAgc                                                                                                                <  2:2325620/150‑1 (MQ=255)
                                       atTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAgc                                                                                                                <  2:3743682/150‑1 (MQ=255)
                                       atTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAgc                                                                                                                >  1:966216/1‑150 (MQ=255)
                                          gAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGcgcg                                                                                                             >  1:2133804/1‑150 (MQ=255)
                                            aTTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGt                                                                                                           <  1:3531856/150‑1 (MQ=255)
                                             tttGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTa                                                                                                          <  1:1982372/150‑1 (MQ=255)
                                               tgttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTaaa                                                                                                        <  2:3667391/150‑1 (MQ=255)
                                                 ttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAtc                                                                                                      <  1:4415813/150‑1 (MQ=255)
                                                    cGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCg                                                                                                   >  2:568700/1‑150 (MQ=255)
                                                      aGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGAt                                                                                                 <  2:4305285/150‑1 (MQ=255)
                                                       gTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATa                                                                                                <  2:1822893/150‑1 (MQ=255)
                                                       gTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATa                                                                                                <  2:2180760/150‑1 (MQ=255)
                                                         tCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTACCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATaaa                                                                                              <  1:4450203/150‑1 (MQ=255)
                                                          ccTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAAt                                                                                             <  2:1653811/150‑1 (MQ=255)
                                                            tGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATgg                                                                                           >  1:4141923/1‑150 (MQ=255)
                                                               aaCAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCgg                                                                                        <  2:629513/150‑1 (MQ=255)
                                                               aaCAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCgg                                                                                        >  1:2294190/1‑150 (MQ=255)
                                                                  aGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTaa                                                                                     >  1:2201458/1‑150 (MQ=255)
                                                                        cACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGtt                                                                               >  1:527721/1‑150 (MQ=255)
                                                                          cTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGtttt                                                                             >  1:84137/1‑150 (MQ=255)
                                                                           ttCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGttttt                                                                            >  1:4285347/1‑150 (MQ=255)
                                                                             ccACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCa                                                                          >  2:3905739/1‑150 (MQ=255)
                                                                              cacaTTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCAt                                                                         >  2:4241619/1‑150 (MQ=255)
                                                                               acaTTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCAtt                                                                        <  1:422531/150‑1 (MQ=255)
                                                                                caTTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCACGATAAATGGCGGTAATTTGTTTTTCAttg                                                                       >  2:927939/1‑150 (MQ=255)
                                                                                      aGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGaca                                                                 <  2:2775455/150‑1 (MQ=255)
                                                                                         ccAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGacacac                                                              <  2:2993174/150‑1 (MQ=255)
                                                                                          cAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACAcc                                                             <  1:958331/150‑1 (MQ=255)
                                                                                           aaaTCTTTAAACCAGTCCGACAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACct                                                            >  1:3681290/1‑150 (MQ=255)
                                                                                            aaTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACctc                                                           >  1:4330759/1‑150 (MQ=255)
                                                                                             aTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACctct                                                          >  1:1434939/1‑150 (MQ=255)
                                                                                              tCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTg                                                         <  1:4263646/150‑1 (MQ=255)
                                                                                               cTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTgg                                                        >  2:3277205/1‑150 (MQ=255)
                                                                                               cTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTgg                                                        <  2:125845/150‑1 (MQ=255)
                                                                                                tttAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTAATTGTTGACACACCTCTGGt                                                       >  2:2473519/1‑150 (MQ=255)
                                                                                                 ttAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTc                                                      <  2:3507326/150‑1 (MQ=255)
                                                                                                   aaaCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCAt                                                    >  2:2405858/1‑150 (MQ=255)
                                                                                                   aaaCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCAt                                                    <  1:2859107/150‑1 (MQ=255)
                                                                                                    aaCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATg                                                   <  2:1085279/150‑1 (MQ=255)
                                                                                                        aGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAg                                               <  1:1007231/150‑1 (MQ=255)
                                                                                                        aGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACATCTGGTCATGATAg                                               <  1:61628/150‑1 (MQ=255)
                                                                                                         gTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGt                                              <  1:1749710/150‑1 (MQ=255)
                                                                                                         gTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGt                                              <  1:1226705/150‑1 (MQ=255)
                                                                                                            cGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTGATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATc                                           >  1:795472/1‑150 (MQ=255)
                                                                                                             gCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCa                                          >  1:2980952/1‑150 (MQ=255)
                                                                                                             gCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCa                                          <  1:739555/150‑1 (MQ=255)
                                                                                                             gCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCGCTGGTCATGATAGTATCa                                          >  1:3873465/1‑150 (MQ=255)
                                                                                                               cagcagAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAat                                        >  2:898269/1‑150 (MQ=255)
                                                                                                                agcagAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAata                                       >  1:1565693/1‑150 (MQ=255)
                                                                                                                 gcagAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAatat                                      <  1:189218/150‑1 (MQ=255)
                                                                                                                 gcagAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAatat                                      <  1:1087045/150‑1 (MQ=255)
                                                                                                                   agagTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTc                                    >  1:190692/1‑150 (MQ=255)
                                                                                                                   agagTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTc                                    >  2:2053362/1‑150 (MQ=255)
                                                                                                                   agagTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACAACTCTGGTAATGATAGTATCAATATTc                                    >  1:1293059/1‑150 (MQ=255)
                                                                                                                    gagTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCa                                   >  2:1748156/1‑150 (MQ=255)
                                                                                                                     agTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCAt                                  <  2:2439065/150‑1 (MQ=255)
                                                                                                                     agTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCAt                                  <  2:136309/150‑1 (MQ=255)
                                                                                                                           tAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACGCTCCTCTGGTCATGATAGTATCAATATTCATGCAGTa                            <  1:380559/150‑1 (MQ=255)
                                                                                                                            aaaTCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTAt                           <  2:4435621/150‑1 (MQ=255)
                                                                                                                            aaaTCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTAt                           >  1:3986342/1‑150 (MQ=255)
                                                                                                                              aTCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTAttt                         >  1:2696317/1‑150 (MQ=255)
                                                                                                                               tCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTa                        <  2:3682441/150‑1 (MQ=255)
                                                                                                                                cTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTAt                       >  1:1446475/1‑150 (MQ=255)
                                                                                                                                 tGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATg                      <  2:3808521/150‑1 (MQ=255)
                                                                                                                                  gCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGa                     <  1:2798557/150‑1 (MQ=255)
                                                                                                                                   cATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGaa                    <  2:2294190/150‑1 (MQ=255)
                                                                                                                                    atagCTTTGATCGAGTGCCTCTTACGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAATATTTATGAAt                   <  2:2916065/147‑1 (MQ=255)
                                                                                                                                    aTTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAAt                   >  2:2658301/1‑150 (MQ=255)
                                                                                                                                     ttGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATa                  <  2:1094810/150‑1 (MQ=255)
                                                                                                                                        cTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaaaa               <  2:79051/150‑1 (MQ=255)
                                                                                                                                           tGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATCGTATCAATATTCATGCAGTATTTATGAATATAAATa            >  2:4103007/1‑150 (MQ=255)
                                                                                                                                            gATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATac           >  2:2132397/1‑150 (MQ=255)
                                                                                                                                             aTCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATaca          >  1:3987977/1‑150 (MQ=255)
                                                                                                                                             aTCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAATATTTATGAATAAAAATaca          <  1:2385504/150‑1 (MQ=255)
                                                                                                                                               cGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACt        <  2:264364/150‑1 (MQ=255)
                                                                                                                                               cGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACt        >  2:3450216/1‑150 (MQ=255)
                                                                                                                                                 aGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTaa      >  2:2361390/1‑150 (MQ=255)
                                                                                                                                                  gTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAAc     <  2:3057376/150‑1 (MQ=255)
                                                                                                                                                     ccTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGtt  <  2:516165/150‑1 (MQ=255)
                                                                                                                                                     ccTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGtt  <  2:1386864/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTT  >  CP009273/4144548‑4144846

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: