Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,821,696 A→T 68.0% L274Q (CTG→CAG)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,821,6960AT68.0% 405.3 / 282.3 287L274Q (CTG→CAG) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base A (60/32);  new base T (100/95);  total (160/127)
Fisher's exact test for biased strand distribution p-value = 3.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 287 total aligned reads displayed.

GAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATT  >  CP009273/1821550‑1821844
                                                                                                                                                  |                                                                                                                                                    
  ggCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAAc                                                                                                                                                 >  2:4245272/1‑150 (MQ=255)
     agcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACtgt                                                                                                                                              <  1:4199039/150‑1 (MQ=255)
         ccAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACt                                                                                                                                          <  2:843199/150‑1 (MQ=255)
         ccAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACt                                                                                                                                          <  1:3414391/150‑1 (MQ=255)
         ccAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGAATCAGACACGGAATCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACt                                                                                                                                          <  1:3689545/150‑1 (MQ=255)
          cAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTg                                                                                                                                         <  1:3201482/150‑1 (MQ=255)
          cAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTg                                                                                                                                         >  2:1296056/1‑150 (MQ=255)
           aaCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTgg                                                                                                                                        >  1:1612936/1‑150 (MQ=255)
           aaCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTgg                                                                                                                                        >  1:1949372/1‑150 (MQ=255)
           aaCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTgg                                                                                                                                        <  1:4239941/150‑1 (MQ=255)
              atcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcga                                                                                                                                     <  2:4084952/150‑1 (MQ=255)
                 atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCCTAAAACCATTGACCAGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgagcg                                                                                                                                  >  2:1341091/1‑150 (MQ=255)
                 atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgagcg                                                                                                                                  >  1:2994864/1‑150 (MQ=255)
                    gAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAAc                                                                                                                               >  2:3412715/1‑150 (MQ=255)
                      aCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGc                                                                                                                             >  1:214857/1‑150 (MQ=255)
                      aCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGc                                                                                                                             >  1:214805/1‑150 (MQ=255)
                      aCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGc                                                                                                                             >  1:4259643/1‑150 (MQ=255)
                        catcatCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCt                                                                                                                           >  2:1946483/1‑150 (MQ=255)
                              cgcgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTtggt                                                                                                                     <  1:2685786/150‑1 (MQ=255)
                               gcgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtggtg                                                                                                                    <  2:2642013/150‑1 (MQ=255)
                                cgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtggtgg                                                                                                                   <  2:3351525/150‑1 (MQ=255)
                                 ggCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                                  >  1:2662723/1‑150 (MQ=255)
                                  gCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCa                                                                                                                 <  1:2535113/150‑1 (MQ=255)
                                    tCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaa                                                                                                               >  1:1045102/1‑150 (MQ=255)
                                     cagcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaaa                                                                                                              <  1:4173271/150‑1 (MQ=255)
                                      agcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAAc                                                                                                             >  1:4460227/1‑150 (MQ=255)
                                        cagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCTCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAg                                                                                                           >  1:1280181/1‑150 (MQ=255)
                                            tGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAcc                                                                                                       >  1:3877392/1‑150 (MQ=255)
                                             ggCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCt                                                                                                      >  2:4327450/1‑150 (MQ=255)
                                              gCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTg                                                                                                     >  1:4334816/1‑150 (MQ=255)
                                                  ttAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCgg                                                                                                 >  1:229418/1‑150 (MQ=255)
                                                     aaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTa                                                                                              >  2:1555340/1‑150 (MQ=255)
                                                     aaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTa                                                                                              >  2:3813475/1‑150 (MQ=255)
                                                     aaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTa                                                                                              >  2:1601514/1‑150 (MQ=255)
                                                          aTAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                         >  2:2037029/1‑150 (MQ=255)
                                                          aTAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                         >  2:88330/1‑150 (MQ=255)
                                                          aTAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                         <  1:4250304/150‑1 (MQ=255)
                                                           tAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACg                                                                                        <  2:867735/150‑1 (MQ=255)
                                                           tAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCACAACAGCACCTGGCGGTTACTCACg                                                                                        >  2:548134/1‑150 (MQ=255)
                                                           tAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACg                                                                                        >  2:546320/1‑150 (MQ=255)
                                                            aGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACgg                                                                                       >  2:728147/1‑150 (MQ=255)
                                                            aGGTAGACAACGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACgg                                                                                       <  1:2039320/150‑1 (MQ=255)
                                                             ggTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGc                                                                                      <  2:2450608/150‑1 (MQ=255)
                                                             ggTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGc                                                                                      <  1:4158327/150‑1 (MQ=255)
                                                             ggTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGc                                                                                      >  1:2224984/1‑150 (MQ=255)
                                                                aGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAAt                                                                                   >  1:2802585/1‑150 (MQ=255)
                                                                   cacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCAc                                                                                <  2:2222980/150‑1 (MQ=255)
                                                                   cacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCAc                                                                                >  2:664005/1‑150 (MQ=255)
                                                                                acGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaaa                                                                   <  1:697388/150‑1 (MQ=255)
                                                                                 cGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaaaa                                                                  >  2:2021228/1‑150 (MQ=255)
                                                                                  ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                 >  1:2496170/1‑150 (MQ=255)
                                                                                    aaaCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGc                                                               >  2:4104390/1‑150 (MQ=255)
                                                                                      aCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                             >  2:2124598/1‑150 (MQ=255)
                                                                                         tGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTgg                                                          >  2:1576345/1‑150 (MQ=255)
                                                                                           aGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTc                                                        <  1:4496282/150‑1 (MQ=255)
                                                                                               gCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATa                                                    >  1:4515167/1‑150 (MQ=255)
                                                                                                 cGGAACTTCTATCGCCATAAATCCATTGACCCGCGCACGCAGGATTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCAGGCGGTTACTCACGTCAATCACGTCATTATGAAAAACGCCCTGGTCGATAAc                                                  >  2:2374476/1‑150 (MQ=255)
                                                                                                 cGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAAc                                                  <  2:4457087/150‑1 (MQ=255)
                                                                                                     aCTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTcc                                              <  1:4217124/150‑1 (MQ=255)
                                                                                                          tATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtt                                         >  1:862811/1‑150 (MQ=255)
                                                                                                          tATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtt                                         >  1:3565053/1‑150 (MQ=255)
                                                                                                          tATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtt                                         >  1:71292/1‑150 (MQ=255)
                                                                                                            tCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtttt                                       >  2:679252/1‑150 (MQ=255)
                                                                                                            tCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtttt                                       >  2:2635921/1‑150 (MQ=255)
                                                                                                            tCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtttt                                       <  1:602405/150‑1 (MQ=255)
                                                                                                            tCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtttt                                       <  1:4041451/150‑1 (MQ=255)
                                                                                                              gCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGc                                     <  2:3912120/150‑1 (MQ=255)
                                                                                                                 aTAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTgg                                  <  2:3315510/150‑1 (MQ=255)
                                                                                                                 aTAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTgg                                  <  2:1248265/150‑1 (MQ=255)
                                                                                                                 aTAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTgg                                  <  1:210628/150‑1 (MQ=255)
                                                                                                                  tAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTggg                                 >  1:428644/1‑150 (MQ=255)
                                                                                                                    aaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                               <  2:3463697/150‑1 (MQ=255)
                                                                                                                    aaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                               >  2:3524624/1‑150 (MQ=255)
                                                                                                                    aaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                               >  2:3525063/1‑150 (MQ=255)
                                                                                                                     aaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGa                              >  1:245287/1‑150 (MQ=255)
                                                                                                                      aCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaa                             <  1:3294007/150‑1 (MQ=255)
                                                                                                                       ccATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaa                            >  1:1776668/1‑150 (MQ=255)
                                                                                                                        cATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGATACTCACGGCAATCACGTCATTATGAAAAACGCCCAGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                           >  2:2357441/1‑150 (MQ=255)
                                                                                                                        cATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                           >  2:1383769/1‑150 (MQ=255)
                                                                                                                        cATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                           >  1:414411/1‑150 (MQ=255)
                                                                                                                        cATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                           >  1:2546429/1‑150 (MQ=255)
                                                                                                                         aTTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAAt                          >  1:58449/1‑150 (MQ=255)
                                                                                                                          ttGACCCGCGCACCCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATc                         <  2:953309/150‑1 (MQ=255)
                                                                                                                               ccGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTg                    <  2:3431874/150‑1 (MQ=255)
                                                                                                                               ccGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTg                    <  2:3431858/150‑1 (MQ=255)
                                                                                                                                  cgcacgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTg                 <  2:2078446/150‑1 (MQ=255)
                                                                                                                                    cacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTggg               <  2:1546063/150‑1 (MQ=255)
                                                                                                                                     acgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTgggg              <  2:4514716/150‑1 (MQ=255)
                                                                                                                                      cgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGa             <  2:4301395/150‑1 (MQ=255)
                                                                                                                                        caGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGAtt           >  1:3109341/1‑150 (MQ=255)
                                                                                                                                        caGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGAtt           >  2:2713694/1‑150 (MQ=255)
                                                                                                                                          ggTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCa         <  1:913557/150‑1 (MQ=255)
                                                                                                                                          ggTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCa         >  1:4249950/1‑150 (MQ=255)
                                                                                                                                           gTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGACAATCACCTGTTGGGGATTCAc        >  2:3697071/1‑150 (MQ=255)
                                                                                                                                           gTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAc        <  2:1113427/150‑1 (MQ=255)
                                                                                                                                           gTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAc        <  2:1746147/150‑1 (MQ=255)
                                                                                                                                              tGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTg     >  2:3012867/1‑150 (MQ=255)
                                                                                                                                              tGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTg     <  2:1885330/150‑1 (MQ=255)
                                                                                                                                              tGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTg     >  1:91566/1‑150 (MQ=255)
                                                                                                                                                 cAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGAtt  <  2:214857/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATT  >  CP009273/1821550‑1821844

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: