Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,820,888 C→A 100% C98F (TGC→TTC)  astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,820,8880CA98.3% 559.3 / 1.1 175C98F (TGC→TTC) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (78/94);  total (81/94)
Fisher's exact test for biased strand distribution p-value = 9.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 178 total aligned reads displayed.

TAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCG  >  CP009273/1820743‑1821032
                                                                                                                                                 |                                                                                                                                                
tAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGaaa                                                                                                                                                >  1:3363636/1‑148 (MQ=255)
tAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGaaa                                                                                                                                                <  2:2702598/148‑1 (MQ=255)
   aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATaa                                                                                                                                             <  2:3158889/148‑1 (MQ=255)
     cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACg                                                                                                                                           >  1:132725/1‑148 (MQ=255)
     cAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACg                                                                                                                                           <  1:557205/148‑1 (MQ=255)
      aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGt                                                                                                                                          >  1:1793029/1‑148 (MQ=255)
       aGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGtt                                                                                                                                         <  1:1723133/148‑1 (MQ=255)
       aGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGtt                                                                                                                                         <  1:1724449/148‑1 (MQ=255)
        ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                        <  2:3117894/148‑1 (MQ=255)
        ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTATCCGCTTTCAGCAAATAGCTGCAAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGttt                                                                                                                                        <  2:2063426/148‑1 (MQ=255)
         gTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTc                                                                                                                                       >  2:1119114/1‑148 (MQ=255)
               gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                                 <  1:335581/148‑1 (MQ=255)
               gcgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGc                                                                                                                                 >  2:2258764/1‑148 (MQ=255)
                cgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCt                                                                                                                                >  1:2997197/1‑148 (MQ=255)
                 gcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCtt                                                                                                                               <  2:1903404/148‑1 (MQ=255)
                  cacaGAGTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTc                                                                                                                              <  2:1872622/148‑1 (MQ=255)
                   acaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCa                                                                                                                             >  1:991005/1‑148 (MQ=255)
                       aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAgcgc                                                                                                                         >  2:566641/1‑148 (MQ=255)
                         tCTTGGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCag                                                                                                                       <  2:3376468/148‑1 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      <  1:1638314/148‑1 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCTGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      <  1:1574140/148‑1 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      >  1:1907050/1‑148 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      <  2:2877736/148‑1 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      >  1:1906531/1‑148 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      >  2:2339651/1‑148 (MQ=255)
                          cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGGTTCCCTTGCTTCAGCGCagg                                                                                                                      <  2:2878519/148‑1 (MQ=255)
                          cTTGGCCCTGGACATCAAAATCTTCCCGCCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCagg                                                                                                                      <  1:3263866/148‑1 (MQ=255)
                            ttGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCaggag                                                                                                                    >  1:883688/1‑148 (MQ=255)
                            ttGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCaggag                                                                                                                    >  2:1690375/1‑148 (MQ=255)
                             tGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCaggagg                                                                                                                   <  2:3229010/148‑1 (MQ=255)
                              gCCCTGGTCATAAAAATCTTCCCGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAggagga                                                                                                                  <  2:3178808/148‑1 (MQ=255)
                              gCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAggagga                                                                                                                  <  1:415995/148‑1 (MQ=255)
                              gCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAggagga                                                                                                                  <  2:1987408/148‑1 (MQ=255)
                              gCCCGGGTCATAACAATCTTCCCGGCACTGTTCCAGTTCGCGGGCCCGACAGGTTTCTCGGCGTTCCGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCAGATCGCTATGGAAAGAACGTTTCCCTTGCTTCAGCGCAggagga                                                                                                                  <  2:2976862/148‑1 (MQ=255)
                                 cTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGGTTc                                                                                                               <  2:477394/148‑1 (MQ=255)
                                 cTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTc                                                                                                               <  2:895462/148‑1 (MQ=255)
                                 cTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTc                                                                                                               >  2:2524393/1‑148 (MQ=255)
                                   ggTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTccc                                                                                                             <  1:2561906/148‑1 (MQ=255)
                                    gTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCg                                                                                                            <  1:1007554/148‑1 (MQ=255)
                                           aaaTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCa                                                                                                     >  2:2851152/1‑148 (MQ=255)
                                             aTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAcc                                                                                                   >  2:2679145/1‑148 (MQ=255)
                                             aTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAcc                                                                                                   >  1:2740950/1‑148 (MQ=255)
                                             aTCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAcc                                                                                                   >  1:2741050/1‑148 (MQ=255)
                                              tCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCa                                                                                                  >  2:872182/1‑148 (MQ=255)
                                               cTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAg                                                                                                 >  2:1282408/1‑148 (MQ=255)
                                               cTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAg                                                                                                 >  1:2225814/1‑148 (MQ=255)
                                                    aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCCGCAAcc                                                                                            >  2:2622408/1‑148 (MQ=255)
                                                    aGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAcc                                                                                            >  2:2286281/1‑148 (MQ=255)
                                                     ggCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCg                                                                                           >  1:2702939/1‑148 (MQ=255)
                                                     ggCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCg                                                                                           >  2:1530956/1‑148 (MQ=255)
                                                             tCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTa                                                                                   >  2:1406461/1‑148 (MQ=255)
                                                                   tCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGa                                                                             <  1:1798809/148‑1 (MQ=255)
                                                                    cgcgcgCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAt                                                                            <  2:3154375/148‑1 (MQ=255)
                                                                      cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                          <  2:40509/148‑1 (MQ=255)
                                                                      cgcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtct                                                                          <  2:2432665/148‑1 (MQ=255)
                                                                             gACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTg                                                                   >  2:3239646/1‑148 (MQ=255)
                                                                              aCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTga                                                                  <  2:2698765/148‑1 (MQ=255)
                                                                                  gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGata                                                              >  1:2474380/1‑148 (MQ=255)
                                                                                   tttCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGatat                                                             >  2:88866/1‑148 (MQ=255)
                                                                                   tttCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGatat                                                             >  2:1667616/1‑148 (MQ=255)
                                                                                    ttCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATc                                                            <  1:3342755/148‑1 (MQ=255)
                                                                                     tctcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATcg                                                           >  2:847386/1‑148 (MQ=255)
                                                                                       tcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATcgcg                                                         <  2:213958/148‑1 (MQ=255)
                                                                                       tcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATcgcg                                                         >  1:1155653/1‑148 (MQ=255)
                                                                                        ccGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCgcgc                                                        >  1:2183531/1‑148 (MQ=255)
                                                                                         cGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGcc                                                       >  2:631345/1‑148 (MQ=255)
                                                                                         cGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGcc                                                       >  1:668847/1‑148 (MQ=255)
                                                                                            tttCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAg                                                    >  2:1501711/1‑148 (MQ=255)
                                                                                              tCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGAAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCa                                                  >  2:1432039/1‑148 (MQ=255)
                                                                                                       aGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCa                                         >  1:2393078/1‑148 (MQ=255)
                                                                                                        gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAg                                        >  2:1495826/1‑148 (MQ=255)
                                                                                                        gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAg                                        <  1:3138041/148‑1 (MQ=255)
                                                                                                         cTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGc                                       <  2:1793029/148‑1 (MQ=255)
                                                                                                            ccAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtc                                    <  2:990966/148‑1 (MQ=255)
                                                                                                             cAACGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCAtct                                   <  1:3381260/148‑1 (MQ=255)
                                                                                                                cGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCca                                <  2:2297520/148‑1 (MQ=255)
                                                                                                                cGACCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCca                                >  2:3287195/1‑148 (MQ=255)
                                                                                                                  aCCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCACCAAGCAACGCGTCCAGCATCTCcaca                              >  1:825197/1‑148 (MQ=255)
                                                                                                                  aCCGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCACCAAGCAACGCGTCCAGCATCTCcaca                              >  1:2607212/1‑148 (MQ=255)
                                                                                                                    cGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAgg                            <  2:1155653/148‑1 (MQ=255)
                                                                                                                    cGCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAgg                            <  2:3410625/148‑1 (MQ=255)
                                                                                                                    cGCCAAACATTCGATTCATATAGCTATGGAATTAACGTTTCCCTTGCTTCAGCGAATGAGGATTCCCGAGGATCACCAGCAACAGACAACGTAACGGGATCTCGACTTGAGATTTCGTGGAATGCAAAGCTTAAATCATCTCCACAgg                            <  2:122030/148‑1 (MQ=255)
                                                                                                                     gCCAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGt                           <  1:1007647/148‑1 (MQ=255)
                                                                                                                      ccAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTg                          <  1:1758207/148‑1 (MQ=255)
                                                                                                                      ccAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTg                          <  1:2129359/148‑1 (MQ=255)
                                                                                                                       cAAACATTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGc                         <  1:22714/148‑1 (MQ=255)
                                                                                                                            aTTCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtct                    >  2:3261305/1‑148 (MQ=255)
                                                                                                                             ttCGATTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtctc                   >  2:1221566/1‑148 (MQ=255)
                                                                                                                                 aTTCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTa               <  2:188690/148‑1 (MQ=255)
                                                                                                                                  ttCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAc              >  1:709219/1‑148 (MQ=255)
                                                                                                                                   tCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAcc             <  2:2579334/148‑1 (MQ=255)
                                                                                                                                   tCATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAcc             >  2:796981/1‑148 (MQ=255)
                                                                                                                                    cATATCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCg            <  1:2258764/148‑1 (MQ=255)
                                                                                                                                       atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt         <  1:2716864/148‑1 (MQ=255)
                                                                                                                                       atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt         <  1:3245835/148‑1 (MQ=255)
                                                                                                                                       atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt         <  1:2524393/148‑1 (MQ=255)
                                                                                                                                       atCGCTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCgtgt         <  1:3245863/148‑1 (MQ=255)
                                                                                                                                           cTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTc     >  1:2317040/1‑148 (MQ=255)
                                                                                                                                           cTATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTc     >  2:1497475/1‑148 (MQ=255)
                                                                                                                                            tATGGAAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTcc    <  1:2321778/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCG  >  CP009273/1820743‑1821032

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: ‑