Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,967,572 A→G 13.9% F601S (TTT→TCT)  aas ← fused 2‑acylglycerophospho‑ethanolamine acyl transferase/acyl‑acyl carrier protein synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,967,5720AG13.9% 486.8 / 67.2 180F601S (TTT→TCT) aasfused 2‑acylglycerophospho‑ethanolamine acyl transferase/acyl‑acyl carrier protein synthetase
Reads supporting (aligned to +/- strand):  ref base A (86/69);  new base G (15/10);  total (101/79)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 186 total aligned reads displayed.

ATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCT  >  CP009273/2967427‑2967719
                                                                                                                                                 |                                                                                                                                                   
aTCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCaaaa                                                                                                                                                   <  2:2050299/148‑1 (MQ=255)
    ctctTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCTCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGca                                                                                                                                               >  1:1684903/1‑148 (MQ=255)
     tctTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGcac                                                                                                                                              >  1:2637923/1‑148 (MQ=255)
      ctTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGcaca                                                                                                                                             <  2:2833864/148‑1 (MQ=255)
        tAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAat                                                                                                                                           >  1:690870/1‑148 (MQ=255)
        tAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAat                                                                                                                                           <  1:2847701/148‑1 (MQ=255)
        tAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAat                                                                                                                                           >  1:1799557/1‑148 (MQ=255)
        tAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAat                                                                                                                                           <  1:229720/148‑1 (MQ=255)
        tAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAat                                                                                                                                           <  1:228744/148‑1 (MQ=255)
         aaTCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAata                                                                                                                                          <  2:1289177/148‑1 (MQ=255)
         aaTCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAata                                                                                                                                          >  2:3086200/1‑148 (MQ=255)
            cGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCg                                                                                                                                       >  2:152237/1‑148 (MQ=255)
             gCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGc                                                                                                                                      <  1:3090738/148‑1 (MQ=255)
                 tGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAgt                                                                                                                                  >  1:908225/1‑148 (MQ=255)
                 tGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAgt                                                                                                                                  >  2:692880/1‑148 (MQ=255)
                  ggCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAgtg                                                                                                                                 <  2:339755/148‑1 (MQ=255)
                  ggCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAgtg                                                                                                                                 <  2:3208801/148‑1 (MQ=255)
                   gCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAgtgt                                                                                                                                >  1:2075015/1‑148 (MQ=255)
                      tGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCAt                                                                                                                             <  1:3219111/148‑1 (MQ=255)
                       ggACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATa                                                                                                                            <  1:2673961/148‑1 (MQ=255)
                       ggACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATa                                                                                                                            >  2:2557044/1‑148 (MQ=255)
                        gACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCACTTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATAc                                                                                                                           >  2:2945688/1‑148 (MQ=255)
                        gACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATAc                                                                                                                           >  2:2945637/1‑148 (MQ=255)
                        gACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATAc                                                                                                                           >  2:2945639/1‑148 (MQ=255)
                         aCTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATAcc                                                                                                                          <  2:2606793/148‑1 (MQ=255)
                         aCTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATAcc                                                                                                                          >  1:518004/1‑148 (MQ=255)
                            ttATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGc                                                                                                                       >  2:310910/1‑148 (MQ=255)
                            ttATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGc                                                                                                                       >  2:1352799/1‑148 (MQ=255)
                            ttATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGc                                                                                                                       >  2:1352569/1‑148 (MQ=255)
                                cTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCcgcg                                                                                                                   >  1:1252400/1‑148 (MQ=255)
                                cTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCcgcg                                                                                                                   >  2:3143491/1‑148 (MQ=255)
                                cTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCcgcg                                                                                                                   >  2:3142547/1‑148 (MQ=255)
                                cTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCcgcg                                                                                                                   >  1:2901347/1‑148 (MQ=255)
                                  ggCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCt                                                                                                                 >  1:1013223/1‑148 (MQ=255)
                                  ggCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCt                                                                                                                 >  1:1011975/1‑148 (MQ=255)
                                   gCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTc                                                                                                                <  2:247697/148‑1 (MQ=255)
                                       aaCACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCAtt                                                                                                            >  1:123505/1‑148 (MQ=255)
                                        acacCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCAttt                                                                                                           <  2:430621/148‑1 (MQ=255)
                                        acacCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCAttt                                                                                                           <  1:1495266/148‑1 (MQ=255)
                                         cacCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTc                                                                                                          >  2:2336796/1‑148 (MQ=255)
                                           ccAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCGAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTcgc                                                                                                        <  2:3192655/148‑1 (MQ=255)
                                           ccAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTcgc                                                                                                        >  2:2212641/1‑148 (MQ=255)
                                             aaGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTcgccg                                                                                                      >  1:1995691/1‑148 (MQ=255)
                                              aGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTcgccgc                                                                                                     >  2:170656/1‑148 (MQ=255)
                                              aGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTcgccgc                                                                                                     >  2:170055/1‑148 (MQ=255)
                                              aGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTcgccgc                                                                                                     >  2:2975887/1‑148 (MQ=255)
                                               gTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCcgcg                                                                                                    >  2:1951371/1‑148 (MQ=255)
                                                tGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGa                                                                                                   >  1:671490/1‑148 (MQ=255)
                                                  ccAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAAc                                                                                                 >  1:755505/1‑148 (MQ=255)
                                                   cAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACa                                                                                                >  2:176417/1‑148 (MQ=255)
                                                    aGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACAt                                                                                               >  2:207763/1‑148 (MQ=255)
                                                      ttgttCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATtc                                                                                             <  1:860453/148‑1 (MQ=255)
                                                      ttgttCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATtc                                                                                             <  2:1799557/148‑1 (MQ=255)
                                                        gttCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATtctc                                                                                           <  1:2974884/148‑1 (MQ=255)
                                                          tCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTcgg                                                                                         <  1:2085649/148‑1 (MQ=255)
                                                                 tttCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTggg                                                                                  >  2:110156/1‑148 (MQ=255)
                                                                  ttCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGc                                                                                 >  2:2070360/1‑148 (MQ=255)
                                                                         gACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCa                                                                          <  1:467123/148‑1 (MQ=255)
                                                                          acacCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAg                                                                         >  1:691536/1‑148 (MQ=255)
                                                                              cATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTaca                                                                     <  2:2409369/148‑1 (MQ=255)
                                                                               aTTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTacac                                                                    <  1:3123844/148‑1 (MQ=255)
                                                                                    gCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTAGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCTTTTCGCCTCGAACATTCTCGGCGGTGGGCACTTCCTGTAAACCTGGc                                                               <  2:1579686/148‑1 (MQ=255)
                                                                                       tGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTtc                                                            >  2:807225/1‑148 (MQ=255)
                                                                                        gCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCTCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCCCTTCCAGTACACCTGGCTtct                                                           >  2:2491378/1‑148 (MQ=255)
                                                                                           aTTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCa                                                        >  2:528569/1‑148 (MQ=255)
                                                                                           aTTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCa                                                        >  1:1309721/1‑148 (MQ=255)
                                                                                                  cAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAg                                                 <  1:529736/148‑1 (MQ=255)
                                                                                                    aaGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGAt                                               >  2:468663/1‑148 (MQ=255)
                                                                                                    aaGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGAt                                               <  1:344326/148‑1 (MQ=255)
                                                                                                    aaGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGAt                                               >  2:1653017/1‑148 (MQ=255)
                                                                                                      gCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAg                                             <  2:2292611/148‑1 (MQ=255)
                                                                                                      gCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAg                                             <  1:2194711/148‑1 (MQ=255)
                                                                                                       cGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGc                                            >  1:1381717/1‑148 (MQ=255)
                                                                                                       cGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGc                                            <  2:2108814/148‑1 (MQ=255)
                                                                                                         ttttGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCg                                          <  1:1158038/148‑1 (MQ=255)
                                                                                                           ttGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGtt                                        <  1:2072308/148‑1 (MQ=255)
                                                                                                               gcgGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCAtt                                    <  2:1441830/148‑1 (MQ=255)
                                                                                                                cgGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTa                                   <  2:981763/148‑1 (MQ=255)
                                                                                                                cgGCACTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTa                                   <  1:3232150/148‑1 (MQ=255)
                                                                                                                  gCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATg                                 >  2:95925/1‑148 (MQ=255)
                                                                                                                  gCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATg                                 >  1:1373632/1‑148 (MQ=255)
                                                                                                                  gCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATg                                 <  2:518004/148‑1 (MQ=255)
                                                                                                                  gCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATg                                 >  2:1859399/1‑148 (MQ=255)
                                                                                                                   cccTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGt                                >  1:399910/1‑148 (MQ=255)
                                                                                                                       gAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCgg                            <  1:1352569/148‑1 (MQ=255)
                                                                                                                        aaTCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGc                           <  1:3063324/148‑1 (MQ=255)
                                                                                                                           cTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCt                        <  1:3182773/148‑1 (MQ=255)
                                                                                                                              cacaAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTc                     >  1:897439/1‑148 (MQ=255)
                                                                                                                              cacaAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTc                     >  1:898355/1‑148 (MQ=255)
                                                                                                                               acaAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCa                    <  1:788128/148‑1 (MQ=255)
                                                                                                                                 aaaGCCCTGCTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGt                  >  2:1297343/1‑148 (MQ=255)
                                                                                                                                   aGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTg                >  2:1675100/1‑148 (MQ=255)
                                                                                                                                   aGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTg                >  2:740955/1‑148 (MQ=255)
                                                                                                                                   aGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTg                >  2:2675793/1‑148 (MQ=255)
                                                                                                                                     cccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCa              >  1:957800/1‑148 (MQ=255)
                                                                                                                                     cccTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCa              >  1:956211/1‑148 (MQ=255)
                                                                                                                                          cTCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGt         >  2:91512/1‑148 (MQ=255)
                                                                                                                                           tCGTCAGAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTc        <  2:2990449/148‑1 (MQ=255)
                                                                                                                                                aaaaCGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGccc   >  2:2786361/1‑148 (MQ=255)
                                                                                                                                                 aaaCGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCct  <  2:1011975/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGACACCATTTCGCCTGCAATTTTGGCAAAGCGTTTTGCGCGGCCCTGAATCTGCACAAAGCCCTGCTCGTCAAAACGCACAATATCGCCAGTGTCATACCAGCCGCGCTCCATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCT  >  CP009273/2967427‑2967719

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: