Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,829,319 C→A 100% I85I (ATC→ATA ydjZ → TVP38/TMEM64 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,829,3190CA98.4% 794.7 / 2.9 249I85I (ATC→ATAydjZTVP38/TMEM64 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (131/114);  total (132/117)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-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 249 total aligned reads displayed.

TCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTG  >  CP009273/1829174‑1829463
                                                                                                                                                 |                                                                                                                                                
tCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAc                                                                                                                                                >  2:1913053/1‑148 (MQ=255)
tCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGGTTTTGATAAc                                                                                                                                                >  1:853871/1‑148 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               >  2:928635/1‑148 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCCCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACAGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTGTTTGATAAccgt                                                                                                                                             >  1:3619970/1‑146 (MQ=255)
    tCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCttt                                                                                                                                            >  2:386070/1‑148 (MQ=255)
    tCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCttt                                                                                                                                            >  1:234546/1‑148 (MQ=255)
      gCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGc                                                                                                                                          <  2:878275/148‑1 (MQ=255)
             tGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGGTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCg                                                                                                                                   >  1:1555068/1‑148 (MQ=255)
              gcggcgTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGt                                                                                                                                  <  2:1157770/148‑1 (MQ=255)
                ggcgTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCg                                                                                                                                <  1:2616761/148‑1 (MQ=255)
                 gcgTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGc                                                                                                                               <  2:96101/148‑1 (MQ=255)
                     ttGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGtt                                                                                                                           <  1:1913053/148‑1 (MQ=255)
                            cgTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGcgcg                                                                                                                    <  1:2945937/148‑1 (MQ=255)
                                gACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCt                                                                                                                <  1:2039306/148‑1 (MQ=255)
                                gACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCt                                                                                                                >  1:2508219/1‑148 (MQ=255)
                                 aCCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTg                                                                                                               <  2:3555487/148‑1 (MQ=255)
                                   cAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGgcg                                                                                                                      <  1:3213200/139‑1 (MQ=255)
                                    aaCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTgggg                                                                                                            >  1:3529225/1‑148 (MQ=255)
                                    aaCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTgggg                                                                                                            <  1:309854/148‑1 (MQ=255)
                                     aCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggc                                                                                                           >  1:1386234/1‑148 (MQ=255)
                                       aGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcgg                                                                                                         <  2:462531/148‑1 (MQ=255)
                                        gggTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcggc                                                                                                        <  1:1024625/148‑1 (MQ=255)
                                        gggTATAGAACGCCTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcggc                                                                                                        <  2:1031174/148‑1 (MQ=255)
                                         ggTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCt                                                                                                       <  1:1068486/148‑1 (MQ=255)
                                            ataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTTCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTgg                                                                                                    >  1:353848/1‑147 (MQ=255)
                                               gAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGt                                                                                                 >  2:2831457/1‑148 (MQ=255)
                                               gAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGGTCTGGGGCGGCTTGCTGt                                                                                                 >  1:1395599/1‑148 (MQ=255)
                                                aaCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTc                                                                                                >  1:2938724/1‑148 (MQ=255)
                                                aaCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTc                                                                                                <  1:1137876/148‑1 (MQ=255)
                                                aaCGCTTTATTCAGTCTTACGGAGCACTGTCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGACATTGCTGCACCGCTACATGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGTGGCTGCTTGCTGTc                                                                                                >  2:3807496/1‑148 (MQ=255)
                                                       tATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGAcc                                                                                         <  1:386070/148‑1 (MQ=255)
                                                          tCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGt                                                                                      <  1:2260230/148‑1 (MQ=255)
                                                          tCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGt                                                                                      <  1:324185/148‑1 (MQ=255)
                                                            aGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTc                                                                                    <  1:639296/148‑1 (MQ=255)
                                                             gTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCg                                                                                   <  2:2693667/148‑1 (MQ=255)
                                                             gTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCg                                                                                   <  1:3311189/148‑1 (MQ=255)
                                                             gTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCg                                                                                   <  2:2966187/148‑1 (MQ=255)
                                                             gTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCg                                                                                   <  1:3499296/148‑1 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  <  1:2074893/148‑1 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  2:1413895/1‑148 (MQ=255)
                                                               cTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAt                                                                                 <  1:942242/148‑1 (MQ=255)
                                                                ttACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCGGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATg                                                                                >  1:1031384/1‑148 (MQ=255)
                                                                 tACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCGGTTTGGCGCGGTCTGGGGCGGCTTGCTGTCGTGGACCCGTTCGCTgg                                                                               >  1:976954/1‑148 (MQ=255)
                                                                  aCGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGc                                                                              <  1:2856140/148‑1 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            >  1:2214710/1‑148 (MQ=255)
                                                                            ggcggcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggcgcggcg                                                                    <  2:1672003/148‑1 (MQ=255)
                                                                            ggcggcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggcgcggcg                                                                    <  1:1344088/148‑1 (MQ=255)
                                                                             gcggcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATTGCCggcgcggcgc                                                                   >  2:1493347/1‑148 (MQ=255)
                                                                             gcggcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCGGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCggcgcggcgc                                                                   >  2:2121926/1‑148 (MQ=255)
                                                                                gcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCtgt                                                                >  2:2127258/1‑148 (MQ=255)
                                                                                 cAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGAGGACCAGTTCGATGGCCGGCGCGGCGCtgtg                                                               <  2:3192488/148‑1 (MQ=255)
                                                                                  aGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGc                                                              <  2:688102/148‑1 (MQ=255)
                                                                                  aGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGc                                                              <  2:925000/148‑1 (MQ=255)
                                                                                       tctcGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCttttt                                                         >  1:2799499/1‑148 (MQ=255)
                                                                                       tctcGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCttttt                                                         <  1:2842034/148‑1 (MQ=255)
                                                                                        ctcGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtttttt                                                        >  2:2859337/1‑148 (MQ=255)
                                                                                         tcGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTa                                                       <  1:2477576/148‑1 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      >  2:487583/1‑148 (MQ=255)
                                                                                          cGTTCTTGTTGATGATATTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATTCGTCGCTGTTTGGAGCGTTCTTGGGCGGCTTGCTTTCGTGGACCAGTTCGATGGCCGGCGCGGCCCTGTGCTTTTTTAt                                                      >  2:2274006/1‑148 (MQ=255)
                                                                                            ttCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCg                                                    >  2:1338301/1‑148 (MQ=255)
                                                                                              cTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGcc                                                  <  2:2004247/148‑1 (MQ=255)
                                                                                                 gttgATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCaga                                               <  1:1722049/148‑1 (MQ=255)
                                                                                                   tgatgaTTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGt                                             >  2:2014321/1‑148 (MQ=255)
                                                                                                   tgatgaTTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGt                                             >  2:754023/1‑148 (MQ=255)
                                                                                                   tgatgaTTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATAGCCAGAGt                                             >  1:3435061/1‑148 (MQ=255)
                                                                                                     atgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGa                                           >  2:956897/1‑148 (MQ=255)
                                                                                                     atgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGa                                           >  2:3370290/1‑148 (MQ=255)
                                                                                                      tgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGAt                                          <  1:1236693/148‑1 (MQ=255)
                                                                                                      tgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGAt                                          <  1:3140569/148‑1 (MQ=255)
                                                                                                          tttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGc                                      >  2:3060363/1‑148 (MQ=255)
                                                                                                          tttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGc                                      >  1:2137244/1‑148 (MQ=255)
                                                                                                           ttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGcc                                     >  2:1408878/1‑148 (MQ=255)
                                                                                                               aGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGa                                 <  1:1288818/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               >  2:1816764/1‑148 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAAAGCGTCGCTGTTTTGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTAGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               >  1:3747757/1‑148 (MQ=255)
                                                                                                                  ccATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTGTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgt                              <  1:2011004/148‑1 (MQ=255)
                                                                                                                    aTTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtgg                            <  2:3111527/148‑1 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           >  2:2502040/1‑148 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  2:2534788/148‑1 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           >  1:2325080/1‑148 (MQ=255)
                                                                                                                     ttGCTGCACCGCTAACTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTAGTGGACCAGTTCCATGGCCGGCGCGGCGCTGTGCTTTTTTATCTCCAGACTGATGGGCCGCGCAttggt                           >  2:3608816/1‑148 (MQ=255)
                                                                                                                      tgctgcACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGTCCAGTTCGATGGCCGGCGCGGCGCGGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggtg                          >  1:423791/1‑148 (MQ=255)
                                                                                                                      tgctgcACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGGTTTTTTATAGCCAGAGTGATGGGCCGCGACgtggtg                          >  1:1804186/1‑148 (MQ=255)
                                                                                                                      tgctgcACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCGGTGCTTTTTTATAGCCCGAGTGATGGGCCGCGAAgtggtg                          >  1:2005567/1‑148 (MQ=255)
                                                                                                                        ctgcACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGACGTGGTGGa                        >  2:2262546/1‑148 (MQ=255)
                                                                                                                         tgcACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaa                       <  2:2938724/148‑1 (MQ=255)
                                                                                                                           cACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaaaa                     >  2:1372006/1‑148 (MQ=255)
                                                                                                                           cACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaaaa                     >  2:3272818/1‑148 (MQ=255)
                                                                                                                            aCCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaaaaa                    <  2:2913264/148‑1 (MQ=255)
                                                                                                                            aCCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaaaaa                    <  1:832397/148‑1 (MQ=255)
                                                                                                                                 tACCTGCGTTTTTGATAACCTTTGCCACTGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAc               <  1:2328321/148‑1 (MQ=255)
                                                                                                                                  aCCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAcc              >  1:613310/1‑148 (MQ=255)
                                                                                                                                    cTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCgg            >  2:2458282/1‑148 (MQ=255)
                                                                                                                                        gtgtttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCaaa        <  2:2415847/145‑1 (MQ=255)
                                                                                                                                          ttttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAAc      >  2:612394/1‑148 (MQ=255)
                                                                                                                                          ttttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAAc      >  2:1236789/1‑148 (MQ=255)
                                                                                                                                          ttttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAAc      >  2:3149787/1‑148 (MQ=255)
                                                                                                                                          ttttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAACGGCAAAAc      >  1:434908/1‑148 (MQ=255)
                                                                                                                                             tGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGGTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGTAAAATTAACCGGCAAAACCGt   >  1:2574008/1‑148 (MQ=255)
                                                                                                                                              gATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTg  <  1:2878967/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTG  >  CP009273/1829174‑1829463

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: