Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,528,263 A→T 7.0% V107E (GTA→GAA)  nanM ← N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,528,2630AT 7.0% 958.4 / 47.9 298V107E (GTA→GAA) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (121/155);  new base T (7/14);  total (129/169)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 303 total aligned reads displayed.

TTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTG  >  CP009273/4528115‑4528411
                                                                                                                                                    |                                                                                                                                                    
ttCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTa                                                                                                                                                      >  1:3486671/1‑149 (MQ=255)
    aaaTAGCCATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTc                                                                                                                                                  <  1:1117491/149‑1 (MQ=255)
         gCCATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTaa                                                                                                                                             <  1:598673/149‑1 (MQ=255)
         gCCATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTaa                                                                                                                                             >  1:3457751/1‑149 (MQ=255)
           cATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAAt                                                                                                                                           >  2:4575100/1‑149 (MQ=255)
            aTTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTTCGTCATTAAATa                                                                                                                                          <  2:3466793/149‑1 (MQ=255)
            aTTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATa                                                                                                                                          <  1:4285924/149‑1 (MQ=255)
             ttGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATAc                                                                                                                                         <  2:1184818/149‑1 (MQ=255)
             ttGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATAc                                                                                                                                         >  1:3403445/1‑149 (MQ=255)
               gAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCt                                                                                                                                       >  2:4048604/1‑149 (MQ=255)
               gAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCt                                                                                                                                       <  2:1263608/149‑1 (MQ=255)
               gAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCt                                                                                                                                       <  2:682021/149‑1 (MQ=255)
                  gATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTTCGTCATTAAATACCTGAg                                                                                                                                    <  1:2937677/149‑1 (MQ=255)
                  gATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAg                                                                                                                                    >  2:1219164/1‑149 (MQ=255)
                  gATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAg                                                                                                                                    <  1:409589/149‑1 (MQ=255)
                   atatTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGt                                                                                                                                   <  2:1653030/149‑1 (MQ=255)
                    tatTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTc                                                                                                                                  <  1:4335171/149‑1 (MQ=255)
                        cTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCAACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGc                                                                                                                              <  1:1533108/149‑1 (MQ=255)
                          ggTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGccc                                                                                                                            >  2:3315889/1‑149 (MQ=255)
                           gTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCt                                                                                                                           <  2:233662/149‑1 (MQ=255)
                              aaCACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGc                                                                                                                        >  1:756575/1‑149 (MQ=255)
                                caccacCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTg                                                                                                                      <  1:2542318/149‑1 (MQ=255)
                                caccacCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTg                                                                                                                      <  1:1484645/149‑1 (MQ=255)
                                caccacCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTg                                                                                                                      >  1:1969124/1‑149 (MQ=255)
                                  ccaccaGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGtt                                                                                                                    <  2:2222345/149‑1 (MQ=255)
                                  ccaccaGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGtt                                                                                                                    <  2:1389862/149‑1 (MQ=255)
                                    accaGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGtttt                                                                                                                  <  1:682706/149‑1 (MQ=255)
                                       aGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGc                                                                                                               >  1:935066/1‑149 (MQ=255)
                                       aGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGc                                                                                                               <  2:3673809/149‑1 (MQ=255)
                                          aaCATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCaa                                                                                                            <  2:4709333/149‑1 (MQ=255)
                                           aCATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAAt                                                                                                           <  2:1817682/149‑1 (MQ=255)
                                             aTAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTTCGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgc                                                                                                         <  1:3582942/149‑1 (MQ=255)
                                             aTAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgc                                                                                                         >  2:724439/1‑149 (MQ=255)
                                              tAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgcc                                                                                                        >  2:4382544/1‑149 (MQ=255)
                                                aGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTGTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgccgc                                                                                                      >  1:2262014/1‑149 (MQ=255)
                                                 gCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgccgcc                                                                                                     >  1:277576/1‑149 (MQ=255)
                                                    ttGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCaaa                                                                                                  >  1:3247354/1‑149 (MQ=255)
                                                      gCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAaca                                                                                                >  2:3940282/1‑149 (MQ=255)
                                                          ttGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATa                                                                                            >  1:3254610/1‑149 (MQ=255)
                                                            gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACa                                                                                          <  2:143609/149‑1 (MQ=255)
                                                                  aaaaGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGc                                                                                    <  1:39582/149‑1 (MQ=255)
                                                                   aaaGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGcc                                                                                   >  1:3106388/1‑149 (MQ=255)
                                                                   aaaGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGcc                                                                                   <  2:1935428/149‑1 (MQ=255)
                                                                    aaGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCa                                                                                  <  1:957208/149‑1 (MQ=255)
                                                                     aGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCAt                                                                                 <  1:792226/149‑1 (MQ=255)
                                                                      gTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATc                                                                                <  1:3675422/149‑1 (MQ=255)
                                                                        cacaTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCaa                                                                              <  2:1110970/149‑1 (MQ=255)
                                                                        cacaTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCaa                                                                              <  1:3366287/149‑1 (MQ=255)
                                                                        cacaTGGCCCGCCATGCCCATCGGCGCGCGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCaa                                                                              >  2:580251/1‑149 (MQ=255)
                                                                           aTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATaa                                                                           >  2:3550851/1‑149 (MQ=255)
                                                                            tGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATaaa                                                                          <  2:3617289/149‑1 (MQ=255)
                                                                            tGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATaaa                                                                          <  1:2723481/149‑1 (MQ=255)
                                                                            tGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATaaa                                                                          >  2:2598823/1‑149 (MQ=255)
                                                                             ggCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAAt                                                                         >  2:4262355/1‑149 (MQ=255)
                                                                               cccGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTTCGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGc                                                                       >  1:3717188/1‑149 (MQ=255)
                                                                                 cGCCATGCCCATCGGCGCGTGCGACATCAATTTCACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCag                                                                     <  1:1410318/149‑1 (MQ=255)
                                                                                  gCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCaga                                                                    >  1:455171/1‑149 (MQ=255)
                                                                                   ccATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTTCGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCagag                                                                   >  2:1778088/1‑149 (MQ=255)
                                                                                   ccATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCagag                                                                   <  1:4382639/149‑1 (MQ=255)
                                                                                    cATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTTCGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAgg                                                                  <  1:522525/149‑1 (MQ=255)
                                                                                     aTGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                 <  2:3365223/149‑1 (MQ=255)
                                                                                        cccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGc                                                              <  1:4048604/149‑1 (MQ=255)
                                                                                         ccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCt                                                             >  2:3251989/1‑149 (MQ=255)
                                                                                         ccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCt                                                             <  1:4695359/149‑1 (MQ=255)
                                                                                         ccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCt                                                             <  1:1669989/149‑1 (MQ=255)
                                                                                         ccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGAAGAGGTTGCt                                                             <  1:4156791/149‑1 (MQ=255)
                                                                                           aTCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                           <  1:1562085/149‑1 (MQ=255)
                                                                                           aTCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                           <  2:3518123/149‑1 (MQ=255)
                                                                                              ggCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGAtc                                                        <  1:3005225/149‑1 (MQ=255)
                                                                                               gcgcGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGAtct                                                       <  1:1245174/149‑1 (MQ=255)
                                                                                                   gTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTg                                                   <  1:2845847/149‑1 (MQ=255)
                                                                                                    tGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTgg                                                  <  2:4127175/149‑1 (MQ=255)
                                                                                                         cATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGc                                             >  1:2003300/1‑149 (MQ=255)
                                                                                                         cATCAATATAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTTCTTGATCTCTTGGTCCga                                             >  2:2840612/1‑148 (MQ=255)
                                                                                                           tCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCa                                           <  1:332920/149‑1 (MQ=255)
                                                                                                             aaTTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTCCGTCATTAAATACCTGAGTCAAGCCCTCGCTGTGTTTGCCAATGCCGCCCAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAgg                                         >  2:2519378/1‑149 (MQ=255)
                                                                                                              aTTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAggg                                        <  1:989805/149‑1 (MQ=255)
                                                                                                              aTTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTATTGCCAATGCCGCCAAACACATACAGATTGCCATCAACAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAggg                                        <  2:4351580/149‑1 (MQ=255)
                                                                                                                  aaCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATg                                    >  1:4754427/1‑149 (MQ=255)
                                                                                                                   aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGc                                   >  2:3034275/1‑149 (MQ=255)
                                                                                                                    cccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCa                                  <  2:2467110/149‑1 (MQ=255)
                                                                                                                     ccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAg                                 >  1:3881912/1‑149 (MQ=255)
                                                                                                                     ccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAg                                 >  1:2533774/1‑149 (MQ=255)
                                                                                                                     ccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAg                                 >  1:2171863/1‑149 (MQ=255)
                                                                                                                      cAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGc                                <  2:4782691/149‑1 (MQ=255)
                                                                                                                      cAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGc                                >  2:1975012/1‑149 (MQ=255)
                                                                                                                       aaCTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCt                               <  2:3890873/149‑1 (MQ=255)
                                                                                                                          tATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAAc                            <  2:1555839/149‑1 (MQ=255)
                                                                                                                                ttttGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCATCTAACGCTGTc                      >  2:2754773/1‑149 (MQ=255)
                                                                                                                                   tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAt                   <  1:1158550/149‑1 (MQ=255)
                                                                                                                                     gggTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAttt                 <  2:28012/149‑1 (MQ=255)
                                                                                                                                       gTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAttttt               <  1:2462735/149‑1 (MQ=255)
                                                                                                                                         tGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTa             <  1:2410047/149‑1 (MQ=255)
                                                                                                                                          gTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTAt            <  1:187986/149‑1 (MQ=255)
                                                                                                                                           tATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGAAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc           >  1:1475223/1‑149 (MQ=255)
                                                                                                                                            aTTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCt          >  1:535204/1‑149 (MQ=255)
                                                                                                                                             tttGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCtt         >  1:298216/1‑149 (MQ=255)
                                                                                                                                              ttGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCttt        <  2:2545128/149‑1 (MQ=255)
                                                                                                                                                   tACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCt   <  2:812332/149‑1 (MQ=255)
                                                                                                                                                    aCGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTCTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCtg  >  1:4106439/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTG  >  CP009273/4528115‑4528411

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: