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,3190CA96.9% 593.2 / 7.7 196I85I (ATC→ATAydjZTVP38/TMEM64 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (88/102);  total (89/107)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-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 196 total aligned reads displayed.

TCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTG  >  CP009273/1829174‑1829463
                                                                                                                                                 |                                                                                                                                                
tCAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAc                                                                                                                                                <  2:692348/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  2:1834932/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  2:1938257/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  1:1928276/148‑1 (MQ=255)
 cAATCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAcc                                                                                                                                               <  2:2541982/148‑1 (MQ=255)
  aaTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCt                                                                                                                                              <  2:2164183/148‑1 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCCACAGGGTATAGAACGCTTTATTCAGTCTTAAGGCGCACTGGCGGCAGTTGTCTCGTTCTTGGTGATGATTTTGCAGGCCAATGCTGCACCGCAACCTGCGTTTTTGATAACCtt                                                                                                                                             >  1:1099191/1‑148 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCtt                                                                                                                                             <  2:1393393/148‑1 (MQ=255)
   aTCGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCtt                                                                                                                                             >  2:1272003/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:1799142/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:134007/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:132943/1‑148 (MQ=255)
     cGCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGGTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTg                                                                                                                                           >  2:1799905/1‑148 (MQ=255)
      gCAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGc                                                                                                                                          <  2:2906646/148‑1 (MQ=255)
       cAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCCGTTGTCTCGTTCTTGTTGATGTTTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTGTGATAACCTTTGcc                                                                                                                                         >  2:2727717/1‑148 (MQ=255)
       cAGCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAAACTTTGcc                                                                                                                                         >  1:3199897/1‑148 (MQ=255)
         gCGTTGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCaa                                                                                                                                       <  1:91042/148‑1 (MQ=255)
            ttGCGGCGTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGc                                                                                                                                    >  2:2117996/1‑148 (MQ=255)
              gcggcgTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGt                                                                                                                                  <  2:486264/148‑1 (MQ=255)
                   gTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTg                                                                                                                             <  1:84754/148‑1 (MQ=255)
                   gTTTGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTg                                                                                                                             <  1:603526/148‑1 (MQ=255)
                    tttGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGt                                                                                                                            >  2:1494871/1‑148 (MQ=255)
                     ttGCCGCCGTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGtt                                                                                                                           <  1:2564284/148‑1 (MQ=255)
                           ccgTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGgcgc                                                                                                                     >  2:2306247/1‑148 (MQ=255)
                           ccgTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGGTCTTGTTGATGATTTTGCAGACCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGgggc                                                                                                                     >  1:2690003/1‑148 (MQ=255)
                             gTGGACCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGt                                                                                                                   <  2:360103/148‑1 (MQ=255)
                                 aCCAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTg                                                                                                               >  2:845954/1‑148 (MQ=255)
                                 aCCAACAGGGTATAGAACGCTTTATTCAGTCATACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTg                                                                                                               <  1:249083/148‑1 (MQ=255)
                                  ccAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTgg                                                                                                              <  1:2833880/148‑1 (MQ=255)
                                   cAACAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTggg                                                                                                             >  1:1481973/1‑148 (MQ=255)
                                     aCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggc                                                                                                           >  2:3524642/1‑148 (MQ=255)
                                     aCAGGGTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGAGTTCTGGggc                                                                                                           >  2:2523255/1‑148 (MQ=255)
                                        gggTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcggc                                                                                                        >  1:2154496/1‑148 (MQ=255)
                                        gggTATAGAACGCTGTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGggcggc                                                                                                        >  2:1427481/1‑148 (MQ=255)
                                         ggTATAGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCt                                                                                                       <  2:972308/148‑1 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     >  1:1571307/1‑148 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     <  2:2042199/148‑1 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     <  2:540649/148‑1 (MQ=255)
                                           tataGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTg                                                                                                     >  2:125464/1‑148 (MQ=255)
                                             taGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCt                                                                                                   <  2:598454/148‑1 (MQ=255)
                                             taGAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCt                                                                                                   <  2:1469391/148‑1 (MQ=255)
                                               gAACGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGt                                                                                                 >  2:2263599/1‑148 (MQ=255)
                                                 aCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCg                                                                                               <  2:3258530/148‑1 (MQ=255)
                                                 aCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCg                                                                                               >  2:1915519/1‑148 (MQ=255)
                                                 aCGCTTTATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCg                                                                                               >  2:1915526/1‑148 (MQ=255)
                                                       tATTCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGAcc                                                                                         <  2:3185283/148‑1 (MQ=255)
                                                          tCAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGt                                                                                      <  1:2235019/148‑1 (MQ=255)
                                                           cAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGtt                                                                                     >  1:3108677/1‑148 (MQ=255)
                                                           cAGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGtt                                                                                     <  2:2125722/148‑1 (MQ=255)
                                                            aGTCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCGTGCTGTCGTGGACCAGGTc                                                                                    >  1:2988111/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  2:2402793/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  1:1776171/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  >  1:3537772/1‑148 (MQ=255)
                                                              tCTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGa                                                                                  <  2:3551077/148‑1 (MQ=255)
                                                               cTTACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGGTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAt                                                                                 >  1:1315672/1‑148 (MQ=255)
                                                                 tACGGCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATgg                                                                               >  1:2596282/1‑148 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            <  1:1897564/148‑1 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            <  2:2424395/148‑1 (MQ=255)
                                                                    ggCGCACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcg                                                                            <  2:2424701/148‑1 (MQ=255)
                                                                      cgcACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggc                                                                          >  1:2220975/1‑148 (MQ=255)
                                                                      cgcACTGGCGGCAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggc                                                                          >  1:864866/1‑148 (MQ=255)
                                                                            ggcggcAGTTGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGAGTCGCTGGTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCcggcgcggcg                                                                    >  2:2501784/1‑148 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            >  1:3485000/1‑148 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            >  2:2717429/1‑148 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            <  1:2014938/148‑1 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCtt                                                            <  1:1088628/148‑1 (MQ=255)
                                                                                    ttGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGAATGCCGGCGCGGCGCAGTGCtt                                                            >  1:1964491/1‑148 (MQ=255)
                                                                                     tGTCTCGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCttt                                                           <  2:2567180/148‑1 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTTGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      >  1:1246614/1‑148 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      >  1:255435/1‑148 (MQ=255)
                                                                                          cGTTCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTAt                                                      <  1:2994265/148‑1 (MQ=255)
                                                                                             tCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGc                                                   >  2:3511352/1‑148 (MQ=255)
                                                                                             tCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGc                                                   >  1:3169861/1‑148 (MQ=255)
                                                                                             tCTTGTTGATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGc                                                   >  2:2440492/1‑148 (MQ=255)
                                                                                                  ttgATGATTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCagag                                              >  1:2004034/1‑148 (MQ=255)
                                                                                                    gatgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTg                                            <  2:334540/148‑1 (MQ=255)
                                                                                                    gatgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTg                                            <  2:727333/148‑1 (MQ=255)
                                                                                                      tgatTTTGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGAt                                          >  2:220396/1‑148 (MQ=255)
                                                                                                          tttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGc                                      <  1:2144019/148‑1 (MQ=255)
                                                                                                          tttGCAGGCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCGGGTTGGGGCGTGCTGGGGCGGCGTGGGGTGGTGGTCCCGGTCGATGGGCGGCGCGGCGCGGTGGTTTTTTTTAGCCCGAGGGGTGGGc                                      >  1:6235/1‑148 (MQ=255)
                                                                                                                ggCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGaa                                <  1:2402774/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               <  1:1799905/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               <  1:1799142/148‑1 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               >  2:1624418/1‑148 (MQ=255)
                                                                                                                 gCCATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAg                               <  2:3503745/148‑1 (MQ=255)
                                                                                                                  ccATTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCGGTGCTTTTTTATCGCCAGAGGGATGGGCCGCGCAgt                              >  1:1840177/1‑148 (MQ=255)
                                                                                                                    aTTGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtgg                            >  2:2682979/1‑148 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  1:1069098/148‑1 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  1:331490/148‑1 (MQ=255)
                                                                                                                     ttGCTGCACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAgtggt                           <  1:332674/148‑1 (MQ=255)
                                                                                                                           cACCGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGaaaa                     <  2:1739307/148‑1 (MQ=255)
                                                                                                                             ccGCTACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAAt                   >  2:2056314/1‑148 (MQ=255)
                                                                                                                                cTACCTGCGTTTTTGATCACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTaa                <  2:3251787/148‑1 (MQ=255)
                                                                                                                                 tACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAc               <  1:906044/148‑1 (MQ=255)
                                                                                                                                 tACCTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAAc               <  1:1157928/148‑1 (MQ=255)
                                                                                                                                    cTGCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCCGAGTGATGGGCCGCGCAGTGGTGGAAAAATTAAcccg            >  2:1206594/1‑146 (MQ=255)
                                                                                                                                      gCGTTTTTGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATAGCCACAGTGATGGGCCGCGCAGGGGTGGTAAAAATAACAGGCa          >  1:2710720/1‑148 (MQ=255)
                                                                                                                                         tttttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCaaaa       <  2:2154496/148‑1 (MQ=255)
                                                                                                                                            ttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCg    <  2:39666/148‑1 (MQ=255)
                                                                                                                                            ttGATAACCTTTGCCAATGCGTCGCTGTTTGGCGCGTTCTGGGGCGGCTTGCTGTCGTGGACCAGTTCGATGGCCGGCGCGGCGCTGTGCTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCg    <  2:3169915/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: