Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,869,780 G→A 87.4% intergenic (‑92/‑160) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,869,7800GA87.4% 903.6 / 107.7 318intergenic (‑92/‑160)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base A (148/130);  total (166/152)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 318 total aligned reads displayed.

TGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTAT  >  CP009273/2869634‑2869926
                                                                                                                                                  |                                                                                                                                                  
tgtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGTTGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAag                                                                                                                                                  >  2:976452/1‑149 (MQ=255)
tgtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAag                                                                                                                                                  >  2:3080888/1‑149 (MQ=255)
 gtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                                                                                                                                                 <  2:2891944/149‑1 (MQ=255)
 gtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagt                                                                                                                                                 <  1:254875/149‑1 (MQ=255)
  tgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagta                                                                                                                                                >  1:3504482/1‑149 (MQ=255)
    gATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagtagt                                                                                                                                              >  1:133555/1‑149 (MQ=255)
     aTGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGtt                                                                                                                                             <  1:4017811/149‑1 (MQ=255)
      tGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTc                                                                                                                                            >  2:2051724/1‑149 (MQ=255)
       gCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTcc                                                                                                                                           <  2:4648424/149‑1 (MQ=255)
       gCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTcc                                                                                                                                           <  1:4187095/149‑1 (MQ=255)
        cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCt                                                                                                                                          <  2:4495461/149‑1 (MQ=255)
         tttCCTCCTCCATTTGCCGCAGGTGAGTAAGTCGTATTTGATACATAACCGGTCCTTTGCAATACAGCTATTTTCTTGCAATCAGTTGATTCGACTATAGGGATCTTATGAGAACGAATGAAATTACCAATTAGAATAAGTTGTTCCtt                                                                                                                                         <  2:4412908/149‑1 (MQ=255)
         tttCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCtt                                                                                                                                         <  2:3405324/149‑1 (MQ=255)
           tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTaa                                                                                                                                       >  2:3887520/1‑149 (MQ=255)
             cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACg                                                                                                                                     <  1:4661169/149‑1 (MQ=255)
              gCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACgg                                                                                                                                    >  1:3762995/1‑149 (MQ=255)
                cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGaa                                                                                                                                  >  2:1514212/1‑149 (MQ=255)
                       tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGAt                                                                                                                           <  2:4609473/149‑1 (MQ=255)
                        gCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGAtt                                                                                                                          <  2:3944523/149‑1 (MQ=255)
                            cAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGc                                                                                                                      <  1:1014254/149‑1 (MQ=255)
                            cAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGc                                                                                                                      <  2:604706/149‑1 (MQ=255)
                             aGGTGAGTAAGTCGTATTTGATCCATACCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                                                                                     <  2:43117/149‑1 (MQ=255)
                             aGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                                                                                     <  2:1144277/149‑1 (MQ=255)
                              ggTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                                                                    <  1:1665873/149‑1 (MQ=255)
                               gTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCaaa                                                                                                                   >  2:128766/1‑149 (MQ=255)
                               gTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCaaa                                                                                                                   <  1:3012262/149‑1 (MQ=255)
                                 gAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGc                                                                                                                 <  1:500281/149‑1 (MQ=255)
                                  aGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCt                                                                                                                >  1:1369730/1‑149 (MQ=255)
                                  aGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGAAATAACGATTTGGCAAAGCt                                                                                                                >  1:3155217/1‑149 (MQ=255)
                                   gTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTa                                                                                                               <  2:2348879/149‑1 (MQ=255)
                                    tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                                                                              <  1:1710717/149‑1 (MQ=255)
                                     aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAat                                                                                                             <  1:1381774/149‑1 (MQ=255)
                                     aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAat                                                                                                             >  2:1954808/1‑149 (MQ=255)
                                            aTTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaa                                                                                                      >  1:3632107/1‑149 (MQ=255)
                                             tttGATCCATAACCGTTCCTTTGCAAGACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                                     <  1:4772095/149‑1 (MQ=255)
                                                gATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTg                                                                                                  >  1:814586/1‑149 (MQ=255)
                                                 aTCCATAACCGGTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGc                                                                                                 <  1:2002159/149‑1 (MQ=255)
                                                  tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCt                                                                                                >  1:1483322/1‑149 (MQ=255)
                                                   ccATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCtt                                                                                               >  2:3516832/1‑149 (MQ=255)
                                                   ccATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCtt                                                                                               <  1:123676/149‑1 (MQ=255)
                                                   ccATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAACGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCtt                                                                                               >  2:1181133/1‑149 (MQ=255)
                                                    cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTa                                                                                              >  2:1639984/1‑149 (MQ=255)
                                                    cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTa                                                                                              >  2:125530/1‑149 (MQ=255)
                                                     aTAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTGAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTaa                                                                                             <  1:4286157/149‑1 (MQ=255)
                                                      tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                                                                            >  2:2470273/1‑149 (MQ=255)
                                                       aaCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                                                                           >  2:931660/1‑149 (MQ=255)
                                                        aCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCCAAGCTAATATCAAAAAGTGCTTAAGGc                                                                                          >  2:4157007/1‑149 (MQ=255)
                                                        aCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGc                                                                                          <  1:1351147/149‑1 (MQ=255)
                                                           gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAcc                                                                                       >  1:3445055/1‑149 (MQ=255)
                                                           gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAcc                                                                                       <  2:2985281/149‑1 (MQ=255)
                                                               cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAt                                                                                   >  2:503607/1‑149 (MQ=255)
                                                                tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAtt                                                                                  >  2:4779793/1‑149 (MQ=255)
                                                                tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAtt                                                                                  >  2:1825767/1‑149 (MQ=255)
                                                                  tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTc                                                                                <  2:1911263/149‑1 (MQ=255)
                                                                    cAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCgg                                                                              <  1:1880553/149‑1 (MQ=255)
                                                                    cAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCgg                                                                              <  2:934735/149‑1 (MQ=255)
                                                                         ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGtt                                                                         <  1:711927/149‑1 (MQ=255)
                                                                         ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGtt                                                                         <  1:3901520/149‑1 (MQ=255)
                                                                         ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGtt                                                                         <  2:2496954/149‑1 (MQ=255)
                                                                            cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAg                                                                      >  2:702030/1‑149 (MQ=255)
                                                                              aTTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGa                                                                    >  2:37214/1‑149 (MQ=255)
                                                                                  tCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAt                                                                >  1:4020563/1‑149 (MQ=255)
                                                                                      gCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGa                                                            <  2:631396/149‑1 (MQ=255)
                                                                                        cATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa                                                          >  2:3365582/1‑149 (MQ=255)
                                                                                         aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                         >  2:3285863/1‑149 (MQ=255)
                                                                                         aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                         <  1:2168587/149‑1 (MQ=255)
                                                                                          tCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAtt                                                        <  2:2718753/149‑1 (MQ=255)
                                                                                            aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGt                                                      >  1:3430264/1‑149 (MQ=255)
                                                                                              aTGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGttt                                                    <  1:2839256/149‑1 (MQ=255)
                                                                                              aTGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGttt                                                    >  2:1416214/1‑149 (MQ=255)
                                                                                               tGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtttt                                                   >  2:3749090/1‑149 (MQ=255)
                                                                                                gTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTa                                                  <  1:1957266/149‑1 (MQ=255)
                                                                                                      aCTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAg                                            >  2:1007552/1‑149 (MQ=255)
                                                                                                       cTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGc                                           >  1:4316281/1‑149 (MQ=255)
                                                                                                        tataGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                          >  1:3835900/1‑149 (MQ=255)
                                                                                                        tataGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                          >  1:1447485/1‑149 (MQ=255)
                                                                                                          taGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGc                                        <  2:857867/149‑1 (MQ=255)
                                                                                                             ggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGt                                     <  1:4118205/149‑1 (MQ=255)
                                                                                                              gAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                                    >  1:4341934/1‑149 (MQ=255)
                                                                                                                  gTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCAt                                <  1:3873191/149‑1 (MQ=255)
                                                                                                                  gTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCAt                                >  2:11036/1‑149 (MQ=255)
                                                                                                                  gTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCAt                                <  1:4502024/149‑1 (MQ=255)
                                                                                                                   tAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                               <  1:931660/149‑1 (MQ=255)
                                                                                                                   tAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                               <  1:1575391/149‑1 (MQ=255)
                                                                                                                    aaGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATAc                              <  2:1334237/149‑1 (MQ=255)
                                                                                                                     agagAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACt                             >  1:2784339/1‑149 (MQ=255)
                                                                                                                        gaACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                          <  1:2741114/149‑1 (MQ=255)
                                                                                                                             aaTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGt                     >  2:2262630/1‑149 (MQ=255)
                                                                                                                              aTGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGtc                    >  2:3487339/1‑148 (MQ=255)
                                                                                                                                gAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAACAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaa                  <  1:2590290/149‑1 (MQ=255)
                                                                                                                                gAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaa                  >  2:2831799/1‑149 (MQ=255)
                                                                                                                                 aaaTTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaa                 >  2:3819548/1‑149 (MQ=255)
                                                                                                                                  aaTTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaa                >  1:2062089/1‑149 (MQ=255)
                                                                                                                                       ccAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAttt           <  2:3643929/149‑1 (MQ=255)
                                                                                                                                        cAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTg          <  2:554385/149‑1 (MQ=255)
                                                                                                                                         aaTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc         <  2:1793786/149‑1 (MQ=255)
                                                                                                                                           ttAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAt       >  2:3172790/1‑149 (MQ=255)
                                                                                                                                             aGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATgg     >  2:261445/1‑149 (MQ=255)
                                                                                                                                              gAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGt    >  2:3186261/1‑149 (MQ=255)
                                                                                                                                               aaTGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTa   <  1:4355137/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTAT  >  CP009273/2869634‑2869926

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: