Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,869,780 G→A 93.0% 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,7800GA93.0% 697.9 / 32.4 227intergenic (‑92/‑160)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (3/13);  new base A (92/119);  total (95/132)
Fisher's exact test for biased strand distribution p-value = 6.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTA  >  CP009273/2869636‑2869925
                                                                                                                                                |                                                                                                                                                 
    tGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGtt                                                                                                                                            >  1:1639982/1‑148 (MQ=255)
      cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTcc                                                                                                                                          <  1:2536412/148‑1 (MQ=255)
      cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTcc                                                                                                                                          >  2:2275049/1‑148 (MQ=255)
      cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTcc                                                                                                                                          <  1:882117/148‑1 (MQ=255)
         tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTa                                                                                                                                       <  2:1142001/148‑1 (MQ=255)
         tCCGCCTACAGTTGCCGCAGGTGAGTAAGTCTTATTTGATCCATAACCGTTACTTTGCAATACAGCTATTTTCTTGCCATCATATGTTTCGACTATAGGTAGCGTAAGAGAAAGAATGAAATTACCAATTAGAATAAGTATTTCCTTa                                                                                                                                       <  2:2794149/148‑1 (MQ=255)
          ccGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTaa                                                                                                                                      >  1:2659769/1‑148 (MQ=255)
           cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAAc                                                                                                                                     <  2:2022316/148‑1 (MQ=255)
              cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGa                                                                                                                                  <  1:358892/148‑1 (MQ=255)
              cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGa                                                                                                                                  <  1:3093904/148‑1 (MQ=255)
               tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGaa                                                                                                                                 >  1:1616606/1‑148 (MQ=255)
               tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGaa                                                                                                                                 >  1:1863942/1‑148 (MQ=255)
                 cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATa                                                                                                                               <  2:2872167/148‑1 (MQ=255)
                  aGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATaa                                                                                                                              >  1:1746104/1‑148 (MQ=255)
                    ttGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACg                                                                                                                            >  2:2225322/1‑148 (MQ=255)
                    ttGCCGCAGGTGAGTAAGTCGTATATGATCCATAACCGATCCTTTGCAATACCGCTATTTTCTTGGCCTCAGATGTTTCGACGATAGGGAGCGTAAGAGAACCAATGAAATTACCAATTAGAATAAGTAGTTCCATAACGGAATaaag                                                                                                                            >  2:1805725/1‑146 (MQ=255)
                     tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCTTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGTGATCGTAAGAGAACGAATGAAATGACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                           <  1:222655/148‑1 (MQ=255)
                     tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                           <  1:3050689/148‑1 (MQ=255)
                     tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                           <  1:222687/148‑1 (MQ=255)
                     tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                           <  1:1407643/148‑1 (MQ=255)
                                aGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGc                                                                                                                >  2:2950213/1‑148 (MQ=255)
                                 gTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCt                                                                                                               >  1:2789637/1‑148 (MQ=255)
                                 gTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCt                                                                                                               >  1:1713167/1‑148 (MQ=255)
                                  tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTa                                                                                                              <  2:1547842/148‑1 (MQ=255)
                                  tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTa                                                                                                              <  1:3176500/148‑1 (MQ=255)
                                   aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                                                                             <  2:1432345/148‑1 (MQ=255)
                                   aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                                                                             <  2:2574504/148‑1 (MQ=255)
                                    aGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAat                                                                                                            <  1:2549534/148‑1 (MQ=255)
                                      tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAatat                                                                                                          >  1:2614953/1‑148 (MQ=255)
                                      tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAatat                                                                                                          >  1:396183/1‑148 (MQ=255)
                                         tATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaa                                                                                                       >  1:1787215/1‑148 (MQ=255)
                                         tATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaa                                                                                                       >  1:2065797/1‑148 (MQ=255)
                                         tATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaa                                                                                                       >  2:2157121/1‑148 (MQ=255)
                                          aTTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaa                                                                                                      <  2:1179493/148‑1 (MQ=255)
                                          aTTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaa                                                                                                      <  2:2830349/148‑1 (MQ=255)
                                            ttGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                                    >  1:720490/1‑148 (MQ=255)
                                             tGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAg                                                                                                   <  1:2324502/148‑1 (MQ=255)
                                               aTCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTg                                                                                                 <  1:407657/148‑1 (MQ=255)
                                               aTCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTg                                                                                                 <  1:652090/148‑1 (MQ=255)
                                                 ccATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCGTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCt                                                                                               >  2:169824/1‑148 (MQ=255)
                                                     aaCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                                                                           >  2:1874291/1‑148 (MQ=255)
                                                       ccGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGc                                                                                         <  1:763844/148‑1 (MQ=255)
                                                        cGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                                                                        >  2:1688623/1‑148 (MQ=255)
                                                        cGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                                                                        <  1:2332986/148‑1 (MQ=255)
                                                        cGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                                                                        <  2:146191/148‑1 (MQ=255)
                                                         gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                                                                       <  1:135684/148‑1 (MQ=255)
                                                          ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAcc                                                                                      >  2:1467855/1‑148 (MQ=255)
                                                          ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTAAACGGAATAACGATTTGGCAACGCTAATATCAAAAAGAGCTTAAGGCAcc                                                                                      >  2:2044317/1‑148 (MQ=255)
                                                            ccTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                                                                    <  1:934648/148‑1 (MQ=255)
                                                            ccTTTGCAAGACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                                                                    <  1:2925648/148‑1 (MQ=255)
                                                                tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAttt                                                                                <  2:1220757/148‑1 (MQ=255)
                                                                  cAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                              <  2:1369079/148‑1 (MQ=255)
                                                                      aCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCg                                                                          <  1:2403588/148‑1 (MQ=255)
                                                                        cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGtt                                                                        <  1:2661437/148‑1 (MQ=255)
                                                                          cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGACAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTa                                                                      <  1:3128577/148‑1 (MQ=255)
                                                                                 cTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAt                                                               <  2:1292353/148‑1 (MQ=255)
                                                                                       aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa                                                         <  2:2520045/148‑1 (MQ=255)
                                                                                       aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa                                                         >  1:1919175/1‑148 (MQ=255)
                                                                                        tCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                        <  1:2700844/148‑1 (MQ=255)
                                                                                         cAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAtt                                                       >  2:1858813/1‑148 (MQ=255)
                                                                                         cAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAtt                                                       >  2:1525648/1‑148 (MQ=255)
                                                                                           gATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGt                                                     <  2:2985113/148‑1 (MQ=255)
                                                                                            aTGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                    >  1:2910173/1‑148 (MQ=255)
                                                                                            aTGTTTCGACTATAGGGAGCGTAACAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                    >  2:23726/1‑148 (MQ=255)
                                                                                              gTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtttt                                                  >  1:1166828/1‑148 (MQ=255)
                                                                                                  cGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAgcg                                              <  2:3251032/148‑1 (MQ=255)
                                                                                                   gACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAgcgc                                             <  2:3032361/148‑1 (MQ=255)
                                                                                                      tataGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACTGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGGTTAGGAAGATTTGAAATTGTTTTATCGCAga                                          >  2:2916246/1‑147 (MQ=255)
                                                                                                          gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                      <  2:2557965/148‑1 (MQ=255)
                                                                                                          gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                      <  1:2275049/148‑1 (MQ=255)
                                                                                                           ggAGCGTAAGAGAACGAATGAAATTACCAATTAGACTAATTAGTTCCTAAACGGAATAACGAATTGGCAAAGATAATATCAAAAAGTGCTTAAGTCACAGGATTTCGGGCTGTGAGGAAGAGTTGAAATAGTTTTAGAGCAGCGGCAg                                     >  2:1494707/1‑148 (MQ=255)
                                                                                                               cGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTc                                 <  2:1639836/148‑1 (MQ=255)
                                                                                                                gTAAGATAACGAATGAAATTAACAATTAGAATAAGTAGTTACTTAACGGAATAACGATTTGGCAAAGCTAATATAAATAAGTGCTTATTGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCa                                <  1:834324/148‑1 (MQ=255)
                                                                                                                gTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCa                                <  1:36571/148‑1 (MQ=255)
                                                                                                                  aaGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                              >  1:2986847/1‑148 (MQ=255)
                                                                                                                     agaACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTa                           <  2:47254/148‑1 (MQ=255)
                                                                                                                     agaACCAATGAAATTACCACTTAGAATAAGTACTTCCTTAACGGAATAACGATATGGCAAAGCTAATATCATAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTa                           <  2:1723565/148‑1 (MQ=255)
                                                                                                                       aaCGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                         <  2:2338644/148‑1 (MQ=255)
                                                                                                                       aaCGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                         <  2:478432/148‑1 (MQ=255)
                                                                                                                       aaCGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                         >  1:115302/1‑148 (MQ=255)
                                                                                                                         cGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGc                       <  1:3262267/148‑1 (MQ=255)
                                                                                                                           aaTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAAAGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCgg                     <  1:380008/148‑1 (MQ=255)
                                                                                                                           aaTGAAATTAACAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAATCTAATATCAAAAAGTGCTTAATGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCgg                     <  1:124861/148‑1 (MQ=255)
                                                                                                                            aTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGt                    <  2:1616606/148‑1 (MQ=255)
                                                                                                                            aTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGt                    <  1:277530/148‑1 (MQ=255)
                                                                                                                            aTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGt                    >  2:1888776/1‑148 (MQ=255)
                                                                                                                               aaaTTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaa                 <  1:2564304/148‑1 (MQ=255)
                                                                                                                                  ttACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaa              >  1:2091678/1‑148 (MQ=255)
                                                                                                                                   tACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaaa             <  2:2087036/148‑1 (MQ=255)
                                                                                                                                   tACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaaa             >  1:1373391/1‑148 (MQ=255)
                                                                                                                                    aCCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAt            <  2:3126707/148‑1 (MQ=255)
                                                                                                                                     ccAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAtt           <  2:2753232/148‑1 (MQ=255)
                                                                                                                                       aaTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTg         >  1:904809/1‑148 (MQ=255)
                                                                                                                                        aTTAGAATAAGTATTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc        <  2:776102/148‑1 (MQ=255)
                                                                                                                                        aTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc        >  1:1064093/1‑148 (MQ=255)
                                                                                                                                          tAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAt      <  1:2158179/148‑1 (MQ=255)
                                                                                                                                            gAATAAGTAGTTCCTTAACGGAATAACGATTTGTCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATgg    <  2:2783324/148‑1 (MQ=255)
                                                                                                                                            gAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATgg    >  2:2327511/1‑148 (MQ=255)
                                                                                                                                            gAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATgg    >  2:2160434/1‑148 (MQ=255)
                                                                                                                                              ttGAGTTGGTCCTTAACGTAATAAAGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGTGCGTTTAGGAAGATTTGAAATTGTTTTAGCTCAGCGGCAGTTTCATACTATGGCGGTAATAAAATTTGCATGGTa  <  2:2749739/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTA  >  CP009273/2869636‑2869925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: