Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 991,204 T→A 9.7% S81C (AGC→TGC)  ssuD ← alkanesulfonate monooxygenase, FMNH(2)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273991,2040TA 9.7% 514.5 / 35.3 171S81C (AGC→TGC) ssuDalkanesulfonate monooxygenase, FMNH(2)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (72/82);  new base A (8/9);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

CGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATTTGTTGCAG  >  CP009273/991057‑991351
                                                                                                                                                   |                                                                                                                                                   
cggaGCGCTCGATATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGAACAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCTGCCTGTCGGGCGGCAACGGTAGGTGAGGTTACGCt                                                                                                                                                     <  1:1887153/145‑1 (MQ=255)
cGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCt                                                                                                                                                     <  2:1590120/148‑1 (MQ=255)
  tAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTgg                                                                                                                                                   >  1:1436470/1‑148 (MQ=255)
  tAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTgg                                                                                                                                                   >  2:2118561/1‑148 (MQ=255)
   aGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTggg                                                                                                                                                  <  2:248519/148‑1 (MQ=255)
       ctcgctATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGc                                                                                                                                              >  2:2629003/1‑148 (MQ=255)
         cgctATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAg                                                                                                                                            >  2:2855261/1‑148 (MQ=255)
           ctATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAggg                                                                                                                                          >  2:2333953/1‑148 (MQ=255)
             aTGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCg                                                                                                                                        >  2:2710752/1‑148 (MQ=255)
              tGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGa                                                                                                                                       <  1:440630/148‑1 (MQ=255)
              tGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGa                                                                                                                                       <  1:1163109/148‑1 (MQ=255)
               gATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGAc                                                                                                                                      >  1:824759/1‑148 (MQ=255)
               gATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGAc                                                                                                                                      >  1:3020112/1‑148 (MQ=255)
               gATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGAc                                                                                                                                      >  1:824857/1‑148 (MQ=255)
                aTCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACa                                                                                                                                     <  2:2983464/148‑1 (MQ=255)
                   aaGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGa                                                                                                                                  <  1:634235/148‑1 (MQ=255)
                     ggAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGaaa                                                                                                                                >  2:2108907/1‑148 (MQ=255)
                         cacTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTc                                                                                                                            >  2:2183247/1‑148 (MQ=255)
                         cacTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTc                                                                                                                            <  2:243345/148‑1 (MQ=255)
                         cacTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTc                                                                                                                            <  2:3220962/148‑1 (MQ=255)
                            tCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAACCTTCAGc                                                                                                                         >  2:3293537/1‑148 (MQ=255)
                            tCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGc                                                                                                                         >  1:1424553/1‑148 (MQ=255)
                             ccGTCGCCTGCCAGCTCTGGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGcc                                                                                                                        >  1:1604735/1‑148 (MQ=255)
                               gTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGc                                                                                                                      <  1:3042918/148‑1 (MQ=255)
                               gTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGc                                                                                                                      <  1:2944610/148‑1 (MQ=255)
                               gTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGc                                                                                                                      <  1:1912106/148‑1 (MQ=255)
                                 cGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTg                                                                                                                    >  1:2556481/1‑148 (MQ=255)
                                 cGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTg                                                                                                                    >  1:434946/1‑148 (MQ=255)
                                      gCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCa                                                                                                               <  1:1325586/148‑1 (MQ=255)
                                       ccAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCAc                                                                                                              >  2:1968367/1‑148 (MQ=255)
                                         aGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCg                                                                                                            >  1:2659348/1‑148 (MQ=255)
                                          gCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCgg                                                                                                           <  2:10485/148‑1 (MQ=255)
                                          gCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCgg                                                                                                           >  2:1179221/1‑148 (MQ=255)
                                               tgtgGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGatca                                                                                                      <  1:2836243/148‑1 (MQ=255)
                                               tgtgGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGatca                                                                                                      <  1:2837726/148‑1 (MQ=255)
                                                  ggATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCg                                                                                                   <  2:428192/148‑1 (MQ=255)
                                                   gATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGa                                                                                                  <  2:497512/148‑1 (MQ=255)
                                                     tCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATg                                                                                                <  1:1273709/148‑1 (MQ=255)
                                                     tCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATg                                                                                                <  1:2088921/148‑1 (MQ=255)
                                                               tgACCAGGTTAAACCACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCGGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAg                                                                                      <  1:320255/148‑1 (MQ=255)
                                                                 aCCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGcc                                                                                    <  2:1744252/148‑1 (MQ=255)
                                                                   cAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCAc                                                                                  >  1:2500501/1‑148 (MQ=255)
                                                                       ttAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCAt                                                                              <  1:2528025/148‑1 (MQ=255)
                                                                       ttAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCAt                                                                              >  2:2030722/1‑148 (MQ=255)
                                                                        tAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATc                                                                             <  2:2663390/148‑1 (MQ=255)
                                                                         aaaCAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCt                                                                            >  2:1588388/1‑148 (MQ=255)
                                                                         aaaCAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCt                                                                            <  1:2233674/148‑1 (MQ=255)
                                                                          aacaacGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCtt                                                                           >  2:1125223/1‑148 (MQ=255)
                                                                           acaacGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTc                                                                          >  1:1389537/1‑148 (MQ=255)
                                                                            caacGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCg                                                                         >  1:1049717/1‑148 (MQ=255)
                                                                            caacGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCg                                                                         >  2:2377960/1‑148 (MQ=255)
                                                                             aacGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGc                                                                        <  2:2542532/148‑1 (MQ=255)
                                                                              acGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCa                                                                       <  1:903771/148‑1 (MQ=255)
                                                                                gCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAgg                                                                     <  1:2315200/148‑1 (MQ=255)
                                                                                   cGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAgc                                                                  <  2:2691937/148‑1 (MQ=255)
                                                                                       cATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGccc                                                              <  1:1186670/148‑1 (MQ=255)
                                                                                       cATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGccc                                                              >  1:1957430/1‑148 (MQ=255)
                                                                                         tttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCg                                                            >  1:1319645/1‑148 (MQ=255)
                                                                                         tttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCg                                                            >  1:3293869/1‑148 (MQ=255)
                                                                                         tttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCg                                                            >  1:1565818/1‑148 (MQ=255)
                                                                                         tttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCg                                                            >  1:2899155/1‑148 (MQ=255)
                                                                                          ttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGt                                                           >  2:2428196/1‑148 (MQ=255)
                                                                                          ttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGt                                                           >  1:2453498/1‑148 (MQ=255)
                                                                                          ttGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGt                                                           >  2:2764253/1‑148 (MQ=255)
                                                                                           tGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGtt                                                          <  2:1957430/148‑1 (MQ=255)
                                                                                              gagaCGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGa                                                       <  2:124539/148‑1 (MQ=255)
                                                                                              gagaCGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGa                                                       <  2:305497/148‑1 (MQ=255)
                                                                                               agaCGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGaa                                                      >  2:32409/1‑148 (MQ=255)
                                                                                               agaCGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGaa                                                      >  1:2764333/1‑148 (MQ=255)
                                                                                               agaCGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGGAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGaa                                                      <  2:2533984/148‑1 (MQ=255)
                                                                                                     tCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGca                                                >  2:743907/1‑148 (MQ=255)
                                                                                                         gCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACAcc                                            >  2:2092688/1‑148 (MQ=255)
                                                                                                         gCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACAcc                                            >  2:507188/1‑148 (MQ=255)
                                                                                                           gTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCgg                                          <  2:1693743/148‑1 (MQ=255)
                                                                                                             ggcggcCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGta                                        >  2:2822520/1‑148 (MQ=255)
                                                                                                             ggcggcCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGta                                        >  1:993448/1‑148 (MQ=255)
                                                                                                               cggcCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGtata                                      <  2:86589/148‑1 (MQ=255)
                                                                                                                   cTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCa                                  <  2:1604735/148‑1 (MQ=255)
                                                                                                                   cTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCa                                  <  1:347859/148‑1 (MQ=255)
                                                                                                                     ggcgggcgGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAg                                <  2:1424553/148‑1 (MQ=255)
                                                                                                                     ggcgggcgGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAg                                <  1:2277010/148‑1 (MQ=255)
                                                                                                                      gcgggcggCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGa                               >  1:731295/1‑148 (MQ=255)
                                                                                                                             gcAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATccg                        <  1:1056799/148‑1 (MQ=255)
                                                                                                                             gcAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATccg                        <  1:2710752/148‑1 (MQ=255)
                                                                                                                              cAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATccgc                       >  1:446550/1‑148 (MQ=255)
                                                                                                                              cAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATccgc                       >  1:446529/1‑148 (MQ=255)
                                                                                                                               aaCGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATccgcc                      >  1:222843/1‑148 (MQ=255)
                                                                                                                                aCGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCCTTGGAATTAGCACACCGGTATAGCCAAGACGATccgccg                     <  1:3146010/148‑1 (MQ=255)
                                                                                                                                 cGGTAGGTGAGGTTACGCAGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCcgccgc                    <  2:2556481/148‑1 (MQ=255)
                                                                                                                                    tAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTg                 <  1:963876/148‑1 (MQ=255)
                                                                                                                                          cggTTACGCAGGCCCGCAGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAAtt          <  1:636606/147‑1 (MQ=255)
                                                                                                                                          aGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAAt           <  2:728209/148‑1 (MQ=255)
                                                                                                                                            gTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAAttt         <  1:689070/148‑1 (MQ=255)
                                                                                                                                            gTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAAttt         <  1:2798303/148‑1 (MQ=255)
                                                                                                                                            gTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAAttt         <  1:2350720/148‑1 (MQ=255)
                                                                                                                                             ttACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATttg        >  1:913501/1‑148 (MQ=255)
                                                                                                                                              tACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCTTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATttgt       <  2:226484/148‑1 (MQ=255)
                                                                                                                                              tACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATttgt       <  1:587592/148‑1 (MQ=255)
                                                                                                                                                 gCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATTTGTTGc    <  1:2118561/148‑1 (MQ=255)
                                                                                                                                                   tGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATTTGTTGCAg  <  1:1588432/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGGTGAGGTTACGCTGGGACGCAGGGCGACAAGAAACTTCAGCCGCTGCGTCACCGGGATCATCGATGCGGCAACCAGCCACGCATCTTCGCAGGAGCGCCCCGTTGGAATTAGCACACCGGTATAGCCAAGACGATCCGCCGCTTGCGCAATTTGTTGCAG  >  CP009273/991057‑991351

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: