Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,664,370 G→T 6.7% A473A (GCG→GCT treF → cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,664,3700GT 6.7% 1372.3 / 50.6 358A473A (GCG→GCTtreFcytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base G (145/189);  new base T (10/14);  total (155/203)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 366 total aligned reads displayed.

TCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCG  >  CP009273/3664222‑3664518
                                                                                                                                                    |                                                                                                                                                    
tCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCg                                                                                                                                                      >  1:2936436/1‑149 (MQ=255)
  ggTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGAt                                                                                                                                                    >  1:1191665/1‑149 (MQ=255)
  ggTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGAt                                                                                                                                                    <  1:842140/149‑1 (MQ=255)
   gTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGAtt                                                                                                                                                   >  1:1862191/1‑149 (MQ=255)
    tATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACACTGGATGGCGATTc                                                                                                                                                  >  2:2467139/1‑149 (MQ=255)
     aTGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCa                                                                                                                                                 <  1:4106515/149‑1 (MQ=255)
       ggCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAgg                                                                                                                                               >  2:2255097/1‑149 (MQ=255)
       ggCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAgg                                                                                                                                               <  2:2395329/149‑1 (MQ=255)
         cGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGa                                                                                                                                             <  1:3490040/149‑1 (MQ=255)
         cGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAAGTGATGGCGATTAAGGGa                                                                                                                                             >  1:2645206/1‑149 (MQ=255)
               aTGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTaaa                                                                                                                                       <  2:145899/149‑1 (MQ=255)
               aTGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTaaa                                                                                                                                       >  1:3065179/1‑149 (MQ=255)
               aTGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTaaa                                                                                                                                       >  1:4477326/1‑149 (MQ=255)
                tGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTaaaa                                                                                                                                      >  2:190452/1‑149 (MQ=255)
                tGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTaaaa                                                                                                                                      >  2:2815941/1‑149 (MQ=255)
                 gAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAAt                                                                                                                                     <  1:2852474/149‑1 (MQ=255)
                   agtggCCGATCGTCTGTCAATCGCCGTGCGCAGTCGGTTTCTGACACCTGGCGGGATTCTGGCAAGCGAGTACTAAACGGGTTTACTGTTGGATAAACCCAACGGCTGGGCTCCGTTACAATGGATGGCGATTCAGGGATTTAAAATGt                                                                                                                                   <  2:1784147/146‑1 (MQ=255)
                   aCAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGt                                                                                                                                   >  1:1669140/1‑149 (MQ=255)
                     aGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTAc                                                                                                                                 <  2:1086732/149‑1 (MQ=255)
                      ggCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACg                                                                                                                                <  2:1511925/149‑1 (MQ=255)
                      ggCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACg                                                                                                                                >  1:189245/1‑149 (MQ=255)
                       gCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACgg                                                                                                                               >  1:1448125/1‑149 (MQ=255)
                        ccGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGc                                                                                                                              <  2:1217026/149‑1 (MQ=255)
                        ccGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGc                                                                                                                              >  2:4278308/1‑149 (MQ=255)
                          gATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGa                                                                                                                            >  2:5137665/1‑149 (MQ=255)
                          gATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGa                                                                                                                            >  1:2703826/1‑149 (MQ=255)
                           aTCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGAt                                                                                                                           <  1:741527/149‑1 (MQ=255)
                            tCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATg                                                                                                                          >  2:4641765/1‑149 (MQ=255)
                              gTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAc                                                                                                                        >  2:658653/1‑149 (MQ=255)
                               tCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAcc                                                                                                                       >  1:3201577/1‑149 (MQ=255)
                                 tGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCtt                                                                                                                     >  2:2176988/1‑149 (MQ=255)
                                   gCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAAACGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACCATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTGCt                                                                                                                   >  1:2028079/1‑149 (MQ=255)
                                     aaaCGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAGCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTgg                                                                                                                 <  2:4718294/149‑1 (MQ=255)
                                      aaCGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTggg                                                                                                                >  2:4862121/1‑149 (MQ=255)
                                      aaCGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTggg                                                                                                                <  1:4852836/149‑1 (MQ=255)
                                       aCGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGt                                                                                                               <  2:1685649/149‑1 (MQ=255)
                                        cGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGtg                                                                                                              <  2:1855026/149‑1 (MQ=255)
                                          ccGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGtgat                                                                                                            <  2:4319500/149‑1 (MQ=255)
                                            gTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGtgatga                                                                                                          >  1:3813458/1‑149 (MQ=255)
                                                gcAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATc                                                                                                      <  1:4813808/149‑1 (MQ=255)
                                                         tACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgc                                                                                             <  1:4051899/149‑1 (MQ=255)
                                                         tACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgc                                                                                             >  2:2448668/1‑149 (MQ=255)
                                                         tACTGACACCTGGAGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgc                                                                                             >  2:590405/1‑149 (MQ=255)
                                                            tGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctgg                                                                                          <  1:4674540/149‑1 (MQ=255)
                                                             gACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctggc                                                                                         >  2:4461265/1‑149 (MQ=255)
                                                             gACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctggc                                                                                         >  2:3836598/1‑149 (MQ=255)
                                                                 ccTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGaa                                                                                     <  1:4388838/149‑1 (MQ=255)
                                                                  cTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAg                                                                                    <  1:4502887/149‑1 (MQ=255)
                                                                    ggCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGAc                                                                                  >  2:4898085/1‑149 (MQ=255)
                                                                       gggATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGt                                                                               <  2:3063312/149‑1 (MQ=255)
                                                                           ttCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAAt                                                                           <  2:3514718/149‑1 (MQ=255)
                                                                            tCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATc                                                                          <  1:1501113/149‑1 (MQ=255)
                                                                             cTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCa                                                                         >  1:2101213/1‑149 (MQ=255)
                                                                             cTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCa                                                                         <  2:5034701/149‑1 (MQ=255)
                                                                               ggCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGt                                                                       <  2:439805/149‑1 (MQ=255)
                                                                                gCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTTCGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGtt                                                                      <  1:2411105/149‑1 (MQ=255)
                                                                                 cAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTc                                                                     <  2:4572643/149‑1 (MQ=255)
                                                                                  aaGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCt                                                                    <  2:4007232/149‑1 (MQ=255)
                                                                                  aaGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCt                                                                    >  2:304277/1‑149 (MQ=255)
                                                                                    gCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTAt                                                                  <  1:134880/149‑1 (MQ=255)
                                                                                    gCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTAt                                                                  >  1:224856/1‑149 (MQ=255)
                                                                                       aGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTg                                                               >  1:1230886/1‑149 (MQ=255)
                                                                                              aaCCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGc                                                        <  1:1022524/149‑1 (MQ=255)
                                                                                                 cGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGcaca                                                     <  1:4254009/149‑1 (MQ=255)
                                                                                                   gTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACaaa                                                   >  1:5122069/1‑149 (MQ=255)
                                                                                                    tGAACAGTGGTATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGTGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                  <  2:3152161/149‑1 (MQ=255)
                                                                                                    tGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                  <  2:1002692/149‑1 (MQ=255)
                                                                                                    tGAACAGGGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                  <  1:3108114/149‑1 (MQ=255)
                                                                                                       aCAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGa                                               <  1:1811585/149‑1 (MQ=255)
                                                                                                         aGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTCCTATCTGGAACAGCACAAACTGATc                                             <  1:935935/149‑1 (MQ=255)
                                                                                                           tGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGa                                           >  1:4230328/1‑149 (MQ=255)
                                                                                                           tGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGa                                           <  1:604407/149‑1 (MQ=255)
                                                                                                            gggATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGaa                                          <  1:4878255/149‑1 (MQ=255)
                                                                                                                tAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAAt                                      <  2:4454335/149‑1 (MQ=255)
                                                                                                                 aaaCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATa                                     >  1:3943372/1‑149 (MQ=255)
                                                                                                                  aaCCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAc                                    >  2:2618942/1‑149 (MQ=255)
                                                                                                                  aaCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAc                                    <  2:2254836/149‑1 (MQ=255)
                                                                                                                  aaCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAc                                    <  2:1571840/149‑1 (MQ=255)
                                                                                                                    cccAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCa                                  <  2:4001307/149‑1 (MQ=255)
                                                                                                                    cccAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCa                                  <  2:832174/149‑1 (MQ=255)
                                                                                                                    cccAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCa                                  <  2:1382537/149‑1 (MQ=255)
                                                                                                                     ccAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCat                                 >  2:3335681/1‑149 (MQ=255)
                                                                                                                      cAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCata                                <  1:3415298/149‑1 (MQ=255)
                                                                                                                       aaCGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCatat                               <  1:3559028/149‑1 (MQ=255)
                                                                                                                       aaCGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCatat                               >  1:3626694/1‑149 (MQ=255)
                                                                                                                        aCGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATAtt                              <  2:1804916/149‑1 (MQ=255)
                                                                                                                        aCGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATAtt                              >  1:4989113/1‑149 (MQ=255)
                                                                                                                               ggCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATg                       <  2:2101213/149‑1 (MQ=255)
                                                                                                                                 cACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGgt                     <  1:2543246/149‑1 (MQ=255)
                                                                                                                                   ccGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGgtgt                   <  2:3743379/149‑1 (MQ=255)
                                                                                                                                        aCAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTcccc              <  2:1987712/149‑1 (MQ=255)
                                                                                                                                        aCAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTcccc              <  2:28609/149‑1 (MQ=255)
                                                                                                                                          aaTGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCAAATTGCCGATGGTGTTCCCcgc            >  2:205115/1‑149 (MQ=255)
                                                                                                                                           atggatggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCcgcg           >  1:5223573/1‑149 (MQ=255)
                                                                                                                                            tggatggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGa          >  1:2394596/1‑149 (MQ=255)
                                                                                                                                             ggatggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGaa         <  2:3272791/149‑1 (MQ=255)
                                                                                                                                              gatggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATGCCATATTGCCGATGGTGTTCCCCGCGAAg        >  1:1098168/1‑149 (MQ=255)
                                                                                                                                                tggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggc      >  1:3282184/1‑149 (MQ=255)
                                                                                                                                                tggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggc      >  1:2466922/1‑149 (MQ=255)
                                                                                                                                                   cGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggcggc   <  1:3704047/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCG  >  CP009273/3664222‑3664518

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