Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,548,437 C→A 6.4% G303C (GGT→TGT)  iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,548,4370CA 6.4% 805.3 / 16.4 218G303C (GGT→TGT) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (114/89);  new base A (7/7);  total (121/97)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

GCCTGGCAGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGT  >  CP009273/4548291‑4548584
                                                                                                                                                  |                                                                                                                                                   
 ccTGGCAGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGa                                                                                                                                                   >  1:2134987/1‑148 (MQ=255)
   tGGCAGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGata                                                                                                                                                 <  1:2570399/148‑1 (MQ=255)
      cAGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATgg                                                                                                                                              <  2:1575745/148‑1 (MQ=255)
      cAGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATgg                                                                                                                                              <  1:2329199/148‑1 (MQ=255)
       aGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATggg                                                                                                                                             <  1:258635/148‑1 (MQ=255)
        gAATTTCGCCTGTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGt                                                                                                                                            <  2:775196/148‑1 (MQ=255)
          aTTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTa                                                                                                                                          <  1:2269766/148‑1 (MQ=255)
            ttCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTaaa                                                                                                                                        <  2:1712647/148‑1 (MQ=255)
            ttCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTaaa                                                                                                                                        >  1:585574/1‑148 (MQ=255)
            ttCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTaaa                                                                                                                                        >  1:2252343/1‑148 (MQ=255)
              cGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACAGATATGGGTTAAAtt                                                                                                                                      <  1:2077420/148‑1 (MQ=255)
                 cTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTccc                                                                                                                                   <  1:3145067/148‑1 (MQ=255)
                  tttCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCt                                                                                                                                  <  1:171442/148‑1 (MQ=255)
                  tttCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCt                                                                                                                                  >  2:1304155/1‑148 (MQ=255)
                       cGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTc                                                                                                                             <  1:2082623/148‑1 (MQ=255)
                             ggTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCgaa                                                                                                                       <  2:2309099/148‑1 (MQ=255)
                               ttAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCgaaga                                                                                                                     <  2:2376390/148‑1 (MQ=255)
                                tAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCgaagaa                                                                                                                    >  1:2550538/1‑148 (MQ=255)
                                   ggAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACgg                                                                                                                 <  1:2684566/148‑1 (MQ=255)
                                    gAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACggc                                                                                                                <  1:2283600/148‑1 (MQ=255)
                                       aCCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACggctgg                                                                                                             >  2:2624482/1‑148 (MQ=255)
                                       aCCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACggctgg                                                                                                             <  1:930339/148‑1 (MQ=255)
                                       aCCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACggctgg                                                                                                             >  1:2866413/1‑148 (MQ=255)
                                        ccGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACggctggc                                                                                                            >  1:2088241/1‑148 (MQ=255)
                                          ggCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTa                                                                                                          <  2:2057878/148‑1 (MQ=255)
                                            ctacgctacTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTAcc                                                                                                        >  1:965080/1‑148 (MQ=255)
                                             tacgctacTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCg                                                                                                       <  1:1131495/148‑1 (MQ=255)
                                              acgctacTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACAGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGt                                                                                                      >  2:728514/1‑148 (MQ=255)
                                              acgatacTGGTGTGCGGTCGCTGGGAATCGCTTATACTGAAATCATTGTCTTGGACCAGCACCTGCTGGGTTTCCAGCAACGTTTCAAACCCGGCAACACCGTTATGGGTTAAATTCCGTTCGTCATCTAAGAACGGCTGGCTACCGt                                                                                                      <  1:474212/148‑1 (MQ=255)
                                                   acTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCg                                                                                                 >  2:2080165/1‑148 (MQ=255)
                                                   acTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCg                                                                                                 >  2:1258235/1‑148 (MQ=255)
                                                   acTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACAGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCg                                                                                                 >  2:5517/1‑148 (MQ=255)
                                                     tGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTc                                                                                               >  2:2798633/1‑148 (MQ=255)
                                                       gTGAGCGGGCGCAGGGCATCTGTGATACTGAAATCATAGTCTTTGACCATCACCTGGACGGTTTCCAGCATCTTTTCAAAACCGGCAACACCTATATGGGTTAAATTCCCTTCGTCATCGAAGATCGGCTGGCTACCGTTGCCGTCgg                                                                                             <  2:718854/148‑1 (MQ=255)
                                                         gAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAg                                                                                           <  1:1874631/148‑1 (MQ=255)
                                                         gAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAg                                                                                           <  1:3014114/148‑1 (MQ=255)
                                                           gCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCt                                                                                         <  2:2017539/148‑1 (MQ=255)
                                                            cGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTg                                                                                        >  2:2765158/1‑148 (MQ=255)
                                                            cGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTg                                                                                        >  2:3223713/1‑148 (MQ=255)
                                                              ggCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAACCCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAg                                                                                      <  1:1204886/148‑1 (MQ=255)
                                                                cgcAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAggg                                                                                    >  1:403510/1‑148 (MQ=255)
                                                                cgcAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAggg                                                                                    <  2:1225778/148‑1 (MQ=255)
                                                                  cAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTg                                                                                  <  1:2661221/148‑1 (MQ=255)
                                                                     ggCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGAcg                                                                               >  2:781145/1‑148 (MQ=255)
                                                                     ggCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGAcg                                                                               >  2:1737909/1‑148 (MQ=255)
                                                                     ggCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGATTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTAGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGAcg                                                                               >  2:2722243/1‑148 (MQ=255)
                                                                        aTCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGt                                                                            >  2:2256286/1‑148 (MQ=255)
                                                                         tCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTg                                                                           >  2:1053382/1‑148 (MQ=255)
                                                                         tCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTg                                                                           >  1:170905/1‑148 (MQ=255)
                                                                         tCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTg                                                                           >  2:2952389/1‑148 (MQ=255)
                                                                          cGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGc                                                                          <  1:86623/148‑1 (MQ=255)
                                                                            cTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCa                                                                        >  2:1962839/1‑148 (MQ=255)
                                                                            cTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCa                                                                        >  1:1418520/1‑148 (MQ=255)
                                                                            cTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCa                                                                        >  2:1963004/1‑148 (MQ=255)
                                                                              gATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGc                                                                      >  2:1032113/1‑148 (MQ=255)
                                                                               aTACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCg                                                                     >  1:1905569/1‑148 (MQ=255)
                                                                               aTACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCg                                                                     >  1:81677/1‑148 (MQ=255)
                                                                                 aCTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGa                                                                   >  1:2931670/1‑148 (MQ=255)
                                                                                   tGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCAACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAAt                                                                 >  2:2169587/1‑148 (MQ=255)
                                                                                    gAAATCATTGTCTTTGACCAGCACCTGCACGGTTTCCAGCAAGGTTTCAAAACAGGAAACACAGATATGGTTTAAATTCCATTCGTCATGGTAGAACGGCTTGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTTCCAGCGGAATa                                                                <  2:1437689/148‑1 (MQ=255)
                                                                                    gAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATa                                                                >  1:252436/1‑148 (MQ=255)
                                                                                    gAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACAGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATa                                                                >  2:2701113/1‑148 (MQ=255)
                                                                                     aaaTCATAGTCTTTGCCCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATAc                                                               <  2:583012/148‑1 (MQ=255)
                                                                                       aTCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATAccc                                                             >  2:3269115/1‑148 (MQ=255)
                                                                                         cATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGc                                                           >  1:1434967/1‑148 (MQ=255)
                                                                                                    aCCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGgcgc                                                >  2:1903518/1‑148 (MQ=255)
                                                                                                      cAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCgg                                              >  1:3206003/1‑148 (MQ=255)
                                                                                                        gCACCTGCACGGTTTCCAGCACCGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGc                                            >  2:3001847/1‑148 (MQ=255)
                                                                                                        gCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGc                                            >  1:24205/1‑148 (MQ=255)
                                                                                                        gCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGc                                            >  2:3162757/1‑148 (MQ=255)
                                                                                                        gCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGc                                            >  2:2909843/1‑148 (MQ=255)
                                                                                                         cACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCa                                           <  2:3135119/148‑1 (MQ=255)
                                                                                                          aCCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAANGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCaa                                          <  1:1680170/148‑1 (MQ=255)
                                                                                                          aCCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCaa                                          <  1:2951733/148‑1 (MQ=255)
                                                                                                           ccTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAAt                                         <  1:3222430/148‑1 (MQ=255)
                                                                                                            cTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATa                                        <  2:1536589/148‑1 (MQ=255)
                                                                                                             tGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATAc                                       <  1:2202091/148‑1 (MQ=255)
                                                                                                              gCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATAcc                                      >  1:1980012/1‑148 (MQ=255)
                                                                                                              gCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATAcc                                      <  2:531935/148‑1 (MQ=255)
                                                                                                               cACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCt                                     <  1:1337766/148‑1 (MQ=255)
                                                                                                                tgggTTTCCAGCAACGTTTCAAAACCGGCAACACTGATATGTGTTAAATTCCCTTCGTCCTCTAAGAACGGCTGGCTACGGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCTGCGGAATACCCGCCTGAACGGCGCGGGCAATACCtt                                    <  1:728476/146‑1 (MQ=255)
                                                                                                                aCGGTTTCCAGCAACGTTTCAAAACCGGCAACACAGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCtt                                    <  1:728759/148‑1 (MQ=255)
                                                                                                                 cGGTTTCCAGCAACGTTTCAAAACCGGCAACACAGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTc                                   <  1:273657/148‑1 (MQ=255)
                                                                                                                          gCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGc                          >  1:2115229/1‑148 (MQ=255)
                                                                                                                            aaCGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGa                        <  1:265567/148‑1 (MQ=255)
                                                                                                                            aaCGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGa                        <  1:266040/148‑1 (MQ=255)
                                                                                                                                tttCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCgg                    >  2:2420600/1‑148 (MQ=255)
                                                                                                                                tttCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCgg                    <  2:34751/148‑1 (MQ=255)
                                                                                                                                tttCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCgg                    >  2:1671395/1‑148 (MQ=255)
                                                                                                                                tttCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCgg                    <  2:2661103/148‑1 (MQ=255)
                                                                                                                                    aaaaCCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCg                <  2:1178723/148‑1 (MQ=255)
                                                                                                                                     aaaCCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGt               >  1:1690226/1‑148 (MQ=255)
                                                                                                                                      aaCCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTc              >  1:2793721/1‑148 (MQ=255)
                                                                                                                                        ccGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCaa            <  1:1297136/148‑1 (MQ=255)
                                                                                                                                          ggCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAAtg          <  2:1090636/148‑1 (MQ=255)
                                                                                                                                                 ccGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTgg   >  1:2305891/1‑148 (MQ=255)
                                                                                                                                                 ccGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTgg   >  2:385675/1‑148 (MQ=255)
                                                                                                                                                 ccGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTgg   >  1:2345113/1‑148 (MQ=255)
                                                                                                                                                  cGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCACTACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGt  >  2:445401/1‑148 (MQ=255)
                                                                                                                                                  cGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGt  <  2:249317/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCTGGCAGAATTTCGCCTTTCCCGGTCAGGTTAAGGAAACCGGCTACGCTACTGGTGAGCGGGCGCAGGGCATCGCTGATACTGAAATCATAGTCTTTGACCAGCACCTGCACGGTTTCCAGCAACGTTTCAAAACCGGCAACACCGATATGGGTTAAATTCCCTTCGTCATCGAAGAACGGCTGGCTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGT  >  CP009273/4548291‑4548584

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: ‑