Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,907 C→A 93.1% E298* (GAA→TAA)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,9070CA93.1% 574.3 / 31.1 203E298* (GAA→TAA) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (9/5);  new base A (107/82);  total (116/87)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 207 total aligned reads displayed.

TGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACC  >  CP009273/2808762‑2809051
                                                                                                                                                 |                                                                                                                                                
 gctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTCACGCGTTTTGGACGAATCGGCGCGTACagttagtt                                                                                                                                               >  1:510883/1‑148 (MQ=255)
   tgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGtttt                                                                                                                                             >  1:966418/1‑148 (MQ=255)
     ttCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCa                                                                                                                                           <  2:1609068/148‑1 (MQ=255)
      tCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCaa                                                                                                                                          >  2:623290/1‑148 (MQ=255)
      tCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCaa                                                                                                                                          >  2:1257342/1‑148 (MQ=255)
      tCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGGTAGTTTTCaa                                                                                                                                          >  2:2906125/1‑148 (MQ=255)
        aTCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATc                                                                                                                                        >  1:2867486/1‑148 (MQ=255)
         tCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTACCGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCt                                                                                                                                       >  1:748224/1‑148 (MQ=255)
         tCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCt                                                                                                                                       >  1:3251665/1‑148 (MQ=255)
          cTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTg                                                                                                                                      <  1:2927150/148‑1 (MQ=255)
            acacCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCa                                                                                                                                    >  1:1331842/1‑148 (MQ=255)
             cacCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAGGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCaa                                                                                                                                   >  1:2259636/1‑148 (MQ=255)
              acCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAac                                                                                                                                  <  2:230884/148‑1 (MQ=255)
                cAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAacac                                                                                                                                >  2:2818234/1‑148 (MQ=255)
                cAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAacac                                                                                                                                >  2:3189763/1‑148 (MQ=255)
                cAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAacac                                                                                                                                >  1:846398/1‑148 (MQ=255)
                cAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAacac                                                                                                                                >  1:1035429/1‑148 (MQ=255)
                 aaTCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACg                                                                                                                               <  1:208065/148‑1 (MQ=255)
                 aaTCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACg                                                                                                                               <  1:1492880/148‑1 (MQ=255)
                 aaTCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACg                                                                                                                               >  2:2983236/1‑148 (MQ=255)
                         agaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGtt                                                                                                                       >  2:889210/1‑148 (MQ=255)
                          gaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTg                                                                                                                      <  2:3105941/148‑1 (MQ=255)
                           aaCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGa                                                                                                                     >  2:2719324/1‑148 (MQ=255)
                            aCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAt                                                                                                                    >  1:1876268/1‑148 (MQ=255)
                               ggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTg                                                                                                                 >  2:2452038/1‑148 (MQ=255)
                                gTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGc                                                                                                                <  2:1026046/148‑1 (MQ=255)
                                 ttGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCa                                                                                                               >  2:1620240/1‑148 (MQ=255)
                                 ttGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCa                                                                                                               >  2:3260540/1‑148 (MQ=255)
                                   gATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGc                                                                                                             >  2:2859274/1‑148 (MQ=255)
                                   gATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGc                                                                                                             <  2:520529/148‑1 (MQ=255)
                                     tGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCct                                                                                                           >  2:2308655/1‑148 (MQ=255)
                                      gTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCctc                                                                                                          >  2:1313055/1‑148 (MQ=255)
                                      gTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCctc                                                                                                          <  1:635368/148‑1 (MQ=255)
                                       tCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCctct                                                                                                         <  1:1163723/148‑1 (MQ=255)
                                        ccAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCtt                                                                                                        >  2:2661818/1‑148 (MQ=255)
                                         cAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTa                                                                                                       <  1:1405625/148‑1 (MQ=255)
                                            cGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTAGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGAAGCCTCTTACCg                                                                                                    >  1:1783553/1‑148 (MQ=255)
                                              aaCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTc                                                                                                  >  2:1866881/1‑148 (MQ=255)
                                                cGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCtt                                                                                                <  2:2425664/148‑1 (MQ=255)
                                                   aCTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtct                                                                                             >  2:1123028/1‑148 (MQ=255)
                                                     ttCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCa                                                                                           >  2:286430/1‑148 (MQ=255)
                                                     ttCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCa                                                                                           >  1:1012612/1‑148 (MQ=255)
                                                      tCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCaa                                                                                          <  2:591071/148‑1 (MQ=255)
                                                       cAATATATTCAATGCCGCCACGTAACAGCGCATCAGAGGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAAt                                                                                         <  2:2259636/148‑1 (MQ=255)
                                                       cAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAAt                                                                                         <  2:19439/148‑1 (MQ=255)
                                                       cAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAAt                                                                                         <  2:2222156/148‑1 (MQ=255)
                                                        aaTATCTTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATa                                                                                        <  2:350153/148‑1 (MQ=255)
                                                         atatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAc                                                                                       >  2:260200/1‑148 (MQ=255)
                                                         atatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAc                                                                                       <  1:709501/148‑1 (MQ=255)
                                                          tatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                      >  1:2680329/1‑148 (MQ=255)
                                                          tatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                      >  2:192138/1‑148 (MQ=255)
                                                           atatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCa                                                                                     >  1:2356066/1‑148 (MQ=255)
                                                           atatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCa                                                                                     <  2:2748390/148‑1 (MQ=255)
                                                           atatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCa                                                                                     >  2:695868/1‑148 (MQ=255)
                                                            tatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCaa                                                                                    <  2:1674859/148‑1 (MQ=255)
                                                             atTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAAt                                                                                   >  2:1115847/1‑148 (MQ=255)
                                                             atTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAAt                                                                                   <  1:2435736/148‑1 (MQ=255)
                                                               tCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCt                                                                                 >  1:1621842/1‑148 (MQ=255)
                                                               tCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCt                                                                                 >  1:3039829/1‑148 (MQ=255)
                                                                    gccgccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcgcg                                                                            >  1:941194/1‑148 (MQ=255)
                                                                        ccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTAct                                                                        >  2:435675/1‑148 (MQ=255)
                                                                         cACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTActc                                                                       >  2:1379528/1‑148 (MQ=255)
                                                                          aCGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTActct                                                                      >  1:3223870/1‑148 (MQ=255)
                                                                           cGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTCACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                     >  1:171110/1‑148 (MQ=255)
                                                                           cGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                     <  2:1035429/148‑1 (MQ=255)
                                                                           cGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                     >  2:2055078/1‑148 (MQ=255)
                                                                           cGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                     <  1:2066835/148‑1 (MQ=255)
                                                                            gTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTc                                                                    <  1:2661818/148‑1 (MQ=255)
                                                                                cAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGAt                                                                <  1:2982450/148‑1 (MQ=255)
                                                                                 aGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACAAATCTTCGCGTACTCTTCCGATg                                                               >  2:1614645/1‑148 (MQ=255)
                                                                                       tCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGttt                                                         <  1:1284918/148‑1 (MQ=255)
                                                                                           aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGAt                                                     >  1:359434/1‑148 (MQ=255)
                                                                                           aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGAt                                                     >  2:1745614/1‑148 (MQ=255)
                                                                                             ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                   >  1:2082564/1‑148 (MQ=255)
                                                                                                 aCTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                               <  2:1325545/148‑1 (MQ=255)
                                                                                                  cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGt                                              <  2:736259/148‑1 (MQ=255)
                                                                                                  cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGt                                              <  2:2877899/148‑1 (MQ=255)
                                                                                                              ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCt                                  >  2:1574241/1‑148 (MQ=255)
                                                                                                                tAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATAAGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCATTATCAATAACCAACTTCGCGTACACTTCCCATTGCGTATGGCTTGACTGTTTTACGCCCGCGAc                                >  1:3064106/1‑148 (MQ=255)
                                                                                                                 aaCGCGTTTTGGACGAATCGGCGCGAACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAAACATCGCGTACTCTTCCCAAGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACg                               >  1:1105426/1‑148 (MQ=255)
                                                                                                                   cgcgTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTa                             >  2:852479/1‑148 (MQ=255)
                                                                                                                     cgTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCNTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACt                           <  1:1680448/148‑1 (MQ=255)
                                                                                                                     cgTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACt                           >  1:2191874/1‑148 (MQ=255)
                                                                                                                           ggACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                     <  2:2202756/148‑1 (MQ=255)
                                                                                                                             aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                   >  1:243911/1‑148 (MQ=255)
                                                                                                                             aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCGTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                   <  2:917626/148‑1 (MQ=255)
                                                                                                                                aaTCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCg                <  2:2669339/148‑1 (MQ=255)
                                                                                                                                aaTCGGCGAGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCg                <  2:2342918/148‑1 (MQ=255)
                                                                                                                                   cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg             <  1:2646154/148‑1 (MQ=255)
                                                                                                                                   cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg             <  1:2645900/148‑1 (MQ=255)
                                                                                                                                   cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg             >  2:1299000/1‑148 (MQ=255)
                                                                                                                                    ggCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGa            <  2:1297239/148‑1 (MQ=255)
                                                                                                                                     gcgcGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGaa           <  1:3189763/148‑1 (MQ=255)
                                                                                                                                      cgcgTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg          <  1:2442313/148‑1 (MQ=255)
                                                                                                                                         gTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTa       <  1:174402/148‑1 (MQ=255)
                                                                                                                                            cAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa    <  1:1620240/148‑1 (MQ=255)
                                                                                                                                            cAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa    <  1:3260540/148‑1 (MQ=255)
                                                                                                                                            cAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa    <  2:1331855/148‑1 (MQ=255)
                                                                                                                                              gTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATaac  >  2:775060/1‑146 (MQ=255)
                                                                                                                                              gTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAcc  >  1:2735078/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACC  >  CP009273/2808762‑2809051

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: