Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,237,654 C→A 14.7% M149I (ATG→ATT uxaC ← uronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,237,6540CA14.7% 869.2 / 85.8 263M149I (ATG→ATTuxaCuronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (105/119);  new base A (24/15);  total (129/134)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

GTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTCCAGTGGTATAA  >  CP009273/3237507‑3237803
                                                                                                                                                   |                                                                                                                                                     
  tAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATct                                                                                                                                                   <  2:385755/150‑1 (MQ=255)
   aaaGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATctg                                                                                                                                                  >  1:941420/1‑150 (MQ=255)
   aaaGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTAATctg                                                                                                                                                  >  2:3617658/1‑150 (MQ=255)
    aaGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATctgc                                                                                                                                                 >  2:3742056/1‑150 (MQ=255)
      ggTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATctgctg                                                                                                                                               <  1:2543074/150‑1 (MQ=255)
       gTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCtgctgc                                                                                                                                              >  1:2717278/1‑150 (MQ=255)
       gTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCtgctgc                                                                                                                                              <  2:3287401/150‑1 (MQ=255)
       gTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCtgctgc                                                                                                                                              >  2:3525049/1‑150 (MQ=255)
         cGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCat                                                                                                                                            <  2:4588494/150‑1 (MQ=255)
            cTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCatgat                                                                                                                                         <  1:820505/150‑1 (MQ=255)
             tGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCatgatg                                                                                                                                        >  1:2821450/1‑150 (MQ=255)
             tGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCatgatg                                                                                                                                        >  1:4441444/1‑150 (MQ=255)
              gTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTATGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCCTCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGc                                                                                                                                       >  2:2260222/1‑150 (MQ=255)
                  gATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCcgcg                                                                                                                                   >  2:4499909/1‑150 (MQ=255)
                   aTGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCcgcgc                                                                                                                                  <  2:3417095/150‑1 (MQ=255)
                      ttGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCcgcgcgcg                                                                                                                               <  1:2679774/150‑1 (MQ=255)
                         aaGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGgag                                                                                                                            <  2:1276174/150‑1 (MQ=255)
                          aGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGgaga                                                                                                                           >  2:270076/1‑150 (MQ=255)
                            gCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAg                                                                                                                         <  2:3127107/150‑1 (MQ=255)
                              tttGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGTCGCGCGCGGAGAAGtt                                                                                                                       <  2:2277248/150‑1 (MQ=255)
                              tttGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGtt                                                                                                                       <  1:1676845/150‑1 (MQ=255)
                               ttGTCCGGTCGCCTGGTCGGCAGCACTTTGATGGTGATAGAGGCGTTTTTGGAGATATCTGCGTGATGCTCCAGTGAATAGATCGGGTCATCGTTGGTGCCGACCATTTTCACGTTCATCTGCTGCTTGTTGCCGCGGGTGGAGTAGTTa                                                                                                                      <  2:445122/150‑1 (MQ=255)
                                  tCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATcc                                                                                                                   >  1:1606741/1‑150 (MQ=255)
                                     ggACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGACTCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCTTGATGCCGCGCGCGGAGAAGTTATCCTgc                                                                                                                >  2:1939218/1‑150 (MQ=255)
                                     ggACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTgc                                                                                                                >  2:1153798/1‑150 (MQ=255)
                                      gACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTgcg                                                                                                               <  2:4308796/150‑1 (MQ=255)
                                        cGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGcc                                                                                                             <  2:264749/150‑1 (MQ=255)
                                        cGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGcc                                                                                                             <  1:3156159/150‑1 (MQ=255)
                                         gCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCca                                                                                                            >  1:37391/1‑150 (MQ=255)
                                         gCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCca                                                                                                            >  2:2136788/1‑150 (MQ=255)
                                         gCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCca                                                                                                            <  1:2060497/150‑1 (MQ=255)
                                         gCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCAGCGCca                                                                                                            >  1:3845292/1‑150 (MQ=255)
                                            aGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCcagca                                                                                                         >  2:3985943/1‑150 (MQ=255)
                                               tCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGtt                                                                                                      <  2:4148445/150‑1 (MQ=255)
                                                cGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTc                                                                                                     >  2:1464985/1‑150 (MQ=255)
                                                cGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTc                                                                                                     >  2:463508/1‑150 (MQ=255)
                                                cGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTc                                                                                                     <  1:4054080/150‑1 (MQ=255)
                                                cGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTc                                                                                                     >  2:4398880/1‑150 (MQ=255)
                                                    agcaCTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTa                                                                                                 <  1:3246284/150‑1 (MQ=255)
                                                      caCTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACa                                                                                               >  1:3868384/1‑150 (MQ=255)
                                                       aCTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACAt                                                                                              <  1:3820171/150‑1 (MQ=255)
                                                          ttGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCg                                                                                           >  2:1781440/1‑150 (MQ=255)
                                                          ttGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCg                                                                                           >  2:4112036/1‑150 (MQ=255)
                                                            gATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGtt                                                                                         <  1:3057983/150‑1 (MQ=255)
                                                            gATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGtt                                                                                         >  2:3258856/1‑150 (MQ=255)
                                                               ggTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCa                                                                                      >  2:1810160/1‑150 (MQ=255)
                                                                 tGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGa                                                                                    >  2:2638483/1‑150 (MQ=255)
                                                                 tGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGa                                                                                    >  1:1476756/1‑150 (MQ=255)
                                                                    aaGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGAttt                                                                                 >  2:1458055/1‑150 (MQ=255)
                                                                        gCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATc                                                                             <  1:2788834/150‑1 (MQ=255)
                                                                        gCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATc                                                                             <  1:882830/150‑1 (MQ=255)
                                                                        gCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATc                                                                             <  1:882569/150‑1 (MQ=255)
                                                                         ccGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCg                                                                            >  1:657846/1‑150 (MQ=255)
                                                                           gTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGc                                                                          >  2:3197167/1‑150 (MQ=255)
                                                                               ttGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGtt                                                                      <  2:3824790/150‑1 (MQ=255)
                                                                               ttGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGtt                                                                      >  2:799843/1‑150 (MQ=255)
                                                                                tGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTg                                                                     <  1:2878691/150‑1 (MQ=255)
                                                                                    gATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACgg                                                                 >  2:1595989/1‑150 (MQ=255)
                                                                                    gATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACgg                                                                 <  2:673379/150‑1 (MQ=255)
                                                                                     aTCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGga                                                                >  2:1577211/1‑150 (MQ=255)
                                                                                     aTCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGga                                                                <  1:4277868/150‑1 (MQ=255)
                                                                                       ctctGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGgaga                                                              >  1:1004481/1‑150 (MQ=255)
                                                                                         ctGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGaaa                                                            <  2:5441/150‑1 (MQ=255)
                                                                                           gCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGc                                                          >  2:944848/1‑150 (MQ=255)
                                                                                            cGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCa                                                         >  2:2671711/1‑150 (MQ=255)
                                                                                             gTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCaa                                                        <  1:3651237/150‑1 (MQ=255)
                                                                                             gTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCaa                                                        <  1:1803303/150‑1 (MQ=255)
                                                                                               gATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAAtt                                                      >  2:1702132/1‑150 (MQ=255)
                                                                                                aTGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAAttt                                                     <  2:3893387/150‑1 (MQ=255)
                                                                                                 tGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTg                                                    <  1:3577223/150‑1 (MQ=255)
                                                                                                  gCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGc                                                   <  2:4289098/150‑1 (MQ=255)
                                                                                                       agagAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTg                                              <  1:2020817/150‑1 (MQ=255)
                                                                                                       agagAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTg                                              >  2:3213784/1‑150 (MQ=255)
                                                                                                          gaATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATa                                           >  2:1661074/1‑150 (MQ=255)
                                                                                                          gaATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATa                                           >  1:242473/1‑150 (MQ=255)
                                                                                                           aaTCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCAACGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATAc                                          >  2:2029453/1‑150 (MQ=255)
                                                                                                             tcgatcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCa                                        <  2:1338154/150‑1 (MQ=255)
                                                                                                             tcgatcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCa                                        >  2:2931562/1‑150 (MQ=255)
                                                                                                               gatcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCaaa                                      >  1:2626360/1‑150 (MQ=255)
                                                                                                                atcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAac                                     >  2:4469559/1‑150 (MQ=255)
                                                                                                                 tcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAacg                                    <  1:2760139/150‑1 (MQ=255)
                                                                                                                 tcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAacg                                    >  2:2453825/1‑150 (MQ=255)
                                                                                                                 tcgGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAacg                                    <  2:109371/150‑1 (MQ=255)
                                                                                                                  cgGGTCATCGGTGGTGCCGACCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAacgg                                   >  1:2579570/1‑150 (MQ=255)
                                                                                                                       cATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACTGACGgc                              <  2:1916232/150‑1 (MQ=255)
                                                                                                                           ggtggtGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAg                          <  1:78210/150‑1 (MQ=255)
                                                                                                                           ggtggtGCCGACCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAg                          >  2:3932029/1‑150 (MQ=255)
                                                                                                                             tggtgCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATCTGCCAGTGATACCAAACGGACGGCGCAGtt                        <  1:3894671/150‑1 (MQ=255)
                                                                                                                               gtgCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCg                      <  1:3934083/150‑1 (MQ=255)
                                                                                                                               gtgCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCg                      >  2:1229722/1‑150 (MQ=255)
                                                                                                                               gtgCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCg                      >  1:2186061/1‑150 (MQ=255)
                                                                                                                                    gACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTg                 <  2:2171696/150‑1 (MQ=255)
                                                                                                                                     aCCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGc                >  1:3007045/1‑150 (MQ=255)
                                                                                                                                     aCCATTTTCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGc                >  1:1295354/1‑150 (MQ=255)
                                                                                                                                        aTTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTc             >  2:2741872/1‑150 (MQ=255)
                                                                                                                                         ttttCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTcc            <  2:1025743/150‑1 (MQ=255)
                                                                                                                                          tttCACGTTAATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTCCa           >  1:1207955/1‑150 (MQ=255)
                                                                                                                                             cACGTTCATCTGCTGCATGAAGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTCCAGTg        <  1:4643936/150‑1 (MQ=255)
                                                                                                                                                 ttCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTCCAGTGGtat    <  1:3458562/150‑1 (MQ=255)
                                                                                                                                                   cATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTCCAGTGGTATaa  >  1:374335/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGATGGTGAAAGAGCCGTCTTTGGCGATCTCTGCGTGATGCTCCAGAGAATCGATCGGGTCATCGGTGGTGCCGACCATTTTCACGTTCATCTGCTGCATGATGCCGCGCGCGGAGAAGTTATCCTGCGCCAGCAGTTCGTTACATTCGTTCCAGATTTCATCGGCAGTTGACGGAGAAAGCAATTTGCCAGTGATACCAAACGGACGGCGCAGTTCGAGGTGCGTCCAGTGGTATAA  >  CP009273/3237507‑3237803

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: