Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,941,341 G→A 100% L343L (CTC→CTT amiC ← N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,941,3410GA98.0% 653.0 / 2.7 205L343L (CTC→CTTamiCN‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (96/103);  total (100/105)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 205 total aligned reads displayed.

GAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTTGT  >  CP009273/2941197‑2941482
                                                                                                                                                |                                                                                                                                             
 aaaGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCt                                                                                                                                           <  1:1173559/148‑1 (MQ=255)
    gTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTAt                                                                                                                                        <  1:3133949/148‑1 (MQ=255)
     tCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAAAGCCGTTTCGACCAGAATGGAGGGCATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTAACAAGCTTAtt                                                                                                                                       >  2:70887/1‑148 (MQ=255)
     tCGCCGTTTTCAGTTTACGCTCGTCCTCCACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTAtt                                                                                                                                       <  2:2050866/148‑1 (MQ=255)
       gCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCa                                                                                                                                     <  2:2971156/148‑1 (MQ=255)
        ccGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAg                                                                                                                                    <  2:1986898/148‑1 (MQ=255)
          gTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTa                                                                                                                                  >  2:19980/1‑148 (MQ=255)
           ttttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTAAGTAc                                                                                                                                 >  1:3334750/1‑148 (MQ=255)
            tttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTAcc                                                                                                                                >  2:1494663/1‑148 (MQ=255)
             ttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCg                                                                                                                               <  1:2547183/148‑1 (MQ=255)
             ttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCg                                                                                                                               <  1:2339025/148‑1 (MQ=255)
                aGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCtt                                                                                                                            <  2:779694/148‑1 (MQ=255)
                aGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCtt                                                                                                                            <  1:2612401/148‑1 (MQ=255)
                  tttACGCTCTTCCTCAACTTTTCTGATAAACGCCGTTTCGACCATAATGGAGGGAATATCTGGTGCCTTTAGTACTGCAAACCCGGCCTGTTCAACTTGATTTTTATTCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                          <  2:83555/148‑1 (MQ=255)
                  tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                          >  1:2847169/1‑148 (MQ=255)
                       gCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCaaa                                                                                                                     <  2:190646/148‑1 (MQ=255)
                       gCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCaaa                                                                                                                     >  1:537106/1‑148 (MQ=255)
                       gCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCaaa                                                                                                                     >  1:2304857/1‑148 (MQ=255)
                        ctctTCCTCAACGTTACTGATACACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAAc                                                                                                                    <  2:2715478/148‑1 (MQ=255)
                         tctTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACt                                                                                                                   <  2:3311701/148‑1 (MQ=255)
                           ttCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTc                                                                                                                 <  1:242562/148‑1 (MQ=255)
                            tCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCa                                                                                                                >  1:1924285/1‑148 (MQ=255)
                                cAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCt                                                                                                            <  2:408755/148‑1 (MQ=255)
                                 aaCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTg                                                                                                           <  2:2142644/148‑1 (MQ=255)
                                  aCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGt                                                                                                          >  2:3085508/1‑148 (MQ=255)
                                  aCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGt                                                                                                          >  2:127739/1‑148 (MQ=255)
                                  aCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGt                                                                                                          >  1:1032074/1‑148 (MQ=255)
                                  aCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGt                                                                                                          >  2:1338064/1‑148 (MQ=255)
                                    gTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCg                                                                                                        >  1:2749938/1‑148 (MQ=255)
                                     ttACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCTGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCgg                                                                                                       <  1:3072938/148‑1 (MQ=255)
                                     ttACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCgg                                                                                                       >  1:972412/1‑148 (MQ=255)
                                     ttACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCgg                                                                                                       >  1:375919/1‑148 (MQ=255)
                                      tACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGc                                                                                                      >  1:3202286/1‑148 (MQ=255)
                                         tGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAAt                                                                                                   >  1:1443403/1‑148 (MQ=255)
                                           aTAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATgg                                                                                                 <  2:2978119/148‑1 (MQ=255)
                                           aTAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATgg                                                                                                 <  2:1600842/148‑1 (MQ=255)
                                             aaaCGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTc                                                                                               >  2:73949/1‑148 (MQ=255)
                                              aaCGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCa                                                                                              >  2:552373/1‑148 (MQ=255)
                                                cGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGc                                                                                            >  1:222650/1‑148 (MQ=255)
                                                  ccGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGa                                                                                          <  1:650970/148‑1 (MQ=255)
                                                   cGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGAc                                                                                         <  1:1752388/148‑1 (MQ=255)
                                                   cGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGAc                                                                                         <  1:1752760/148‑1 (MQ=255)
                                                    gTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACt                                                                                        >  2:1008582/1‑148 (MQ=255)
                                                     tttCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTg                                                                                       <  2:1352501/148‑1 (MQ=255)
                                                           ccAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCat                                                                                 >  1:2693300/1‑148 (MQ=255)
                                                            cAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCata                                                                                >  1:1407132/1‑148 (MQ=255)
                                                                 gggAGGTAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCCGCGACTGTACCATATCGaa                                                                           <  1:2083073/147‑1 (MQ=255)
                                                                   gAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACa                                                                         >  1:1634649/1‑148 (MQ=255)
                                                                     gggAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATg                                                                       >  1:2882385/1‑148 (MQ=255)
                                                                     gggAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATg                                                                       <  1:3103919/148‑1 (MQ=255)
                                                                      ggAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATgg                                                                      >  1:731401/1‑148 (MQ=255)
                                                                        aaTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGCCAATGGTCAGCGACTGTACCATATCGAACATGgtg                                                                    >  1:51256/1‑148 (MQ=255)
                                                                          tatCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGtgtg                                                                  <  1:280114/148‑1 (MQ=255)
                                                                           atCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                 <  2:2389884/148‑1 (MQ=255)
                                                                           atCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                 >  2:891056/1‑148 (MQ=255)
                                                                              tGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCg                                                              <  2:2760726/148‑1 (MQ=255)
                                                                                gTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGAc                                                            <  2:1175097/148‑1 (MQ=255)
                                                                                gTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGAc                                                            <  2:683340/148‑1 (MQ=255)
                                                                                     tttAGTACGGCAACCCCGGCTTGTTCAACTTTATTTTTTTGCAGATTGTTGATTTTACAAAGATTATTCAGAACCGCTTTACCAAACTTCAGGCTGTCGGAAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGc                                                       <  1:3266375/148‑1 (MQ=255)
                                                                                     tttAGTACGGCAAACCCGGCCTTTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGc                                                       <  1:3512437/148‑1 (MQ=255)
                                                                                     tttAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCATGCTTATTCTGTACCGCTTTCCCAAACTTCCGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGc                                                       <  1:1567079/148‑1 (MQ=255)
                                                                                          tACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCa                                                  >  1:1160296/1‑148 (MQ=255)
                                                                                            cGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCAcc                                                >  2:581704/1‑148 (MQ=255)
                                                                                             ggCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCg                                               <  2:2354771/148‑1 (MQ=255)
                                                                                             ggCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCg                                               <  1:2569764/148‑1 (MQ=255)
                                                                                              gCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGc                                              <  2:3018444/148‑1 (MQ=255)
                                                                                              gCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGc                                              >  2:2142226/1‑148 (MQ=255)
                                                                                              gCAAACCCGGCCTGTTCAACTTGAATTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGACAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGc                                              <  1:1263223/148‑1 (MQ=255)
                                                                                               cAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCt                                             <  2:1098224/148‑1 (MQ=255)
                                                                                               cAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACAAAGCTTATTCAGTACCGCTATACCAAACTTCATGCTGTCGGCAATGGTCAGCGACTGTACCAAATCGAACAAGGTGTGGTCGACATAGCGGTCACCGCt                                             >  1:864354/1‑148 (MQ=255)
                                                                                                aaaCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtt                                            <  2:3233595/148‑1 (MQ=255)
                                                                                                 aaCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCttt                                           <  1:3260741/148‑1 (MQ=255)
                                                                                                  aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtttt                                          >  1:848779/1‑148 (MQ=255)
                                                                                                  aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtttt                                          >  2:3494277/1‑148 (MQ=255)
                                                                                                  aCCCGGCCTGTTCAACTTGATTTTTAAGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGTTAAGCGACTGTAACAAATCGAAAATGGTGTGGTCGACATATCGGTCACCGCttta                                          >  2:3200627/1‑147 (MQ=255)
                                                                                                     cGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCCGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCt                                       <  2:2178973/148‑1 (MQ=255)
                                                                                                     cGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCt                                       >  1:812506/1‑148 (MQ=255)
                                                                                                             tCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccac                               >  2:3508188/1‑148 (MQ=255)
                                                                                                               aaCTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccacca                             >  1:1952650/1‑148 (MQ=255)
                                                                                                                   tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                         >  2:1805242/1‑148 (MQ=255)
                                                                                                                     aTTTTTATGGAGGTTGTTGTTTTTACCTAGGTTATTCAGTTCCGCTTTTCCAATGTTGAGGGTGTCGGCATTTTTCAGCGACTGTACCATATCGAACATTGTTTGGTCGACATAGCGGTCACCGATTTTGCTCACGCCACCAATCTAg                       <  2:377699/148‑1 (MQ=255)
                                                                                                                      tttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGt                      >  2:2374448/1‑148 (MQ=255)
                                                                                                                         ttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCg                   <  2:2550712/148‑1 (MQ=255)
                                                                                                                           aTGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAg                 >  2:2456553/1‑148 (MQ=255)
                                                                                                                              cAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCg              <  2:972412/148‑1 (MQ=255)
                                                                                                                              cAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCg              <  2:2653911/148‑1 (MQ=255)
                                                                                                                                 cTTGTTTATTTTACCAAGCTTATTCAGTACCGCTTTCCCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTAACGCCACCAATCAAGTCCGAGGCGTTc           <  1:321764/148‑1 (MQ=255)
                                                                                                                                 cTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTc           <  1:127797/148‑1 (MQ=255)
                                                                                                                                    gttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTgg        <  2:804246/148‑1 (MQ=255)
                                                                                                                                     ttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAACCTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTggg       <  1:1884550/148‑1 (MQ=255)
                                                                                                                                      tgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGt      <  2:537106/148‑1 (MQ=255)
                                                                                                                                      tgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGt      <  2:2594117/148‑1 (MQ=255)
                                                                                                                                      tgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGt      <  2:2304857/148‑1 (MQ=255)
                                                                                                                                       gATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGtt     <  1:146569/148‑1 (MQ=255)
                                                                                                                                       gATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGtt     <  2:1147857/148‑1 (MQ=255)
                                                                                                                                        aTTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGttt    >  1:798619/1‑148 (MQ=255)
                                                                                                                                         ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg   >  1:2911875/1‑148 (MQ=255)
                                                                                                                                         ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg   <  1:2594052/148‑1 (MQ=255)
                                                                                                                                         ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg   <  1:835335/148‑1 (MQ=255)
                                                                                                                                         ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg   <  2:640829/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTTGT  >  CP009273/2941197‑2941482

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: