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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,821,6960AT89.3% 774.3 / 66.9 250L274Q (CTG→CAG) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base A (14/13);  new base T (120/103);  total (134/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTG  >  CP009273/1821549‑1821841
                                                                                                                                                   |                                                                                                                                                 
tGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGt                                                                                                                                                 <  1:3037330/150‑1 (MQ=255)
   ggCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAAc                                                                                                                                              >  1:163057/1‑150 (MQ=255)
     cagcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACtg                                                                                                                                            >  2:1001834/1‑150 (MQ=255)
        caCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGa                                                                                                                                         <  1:1236518/150‑1 (MQ=255)
         aCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGAc                                                                                                                                        >  2:3204823/1‑150 (MQ=255)
               atcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcga                                                                                                                                  >  1:791433/1‑150 (MQ=255)
                tcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgag                                                                                                                                 >  1:2340673/1‑150 (MQ=255)
                 catGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgagc                                                                                                                                <  2:211729/150‑1 (MQ=255)
                 catGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgagc                                                                                                                                <  2:1298940/150‑1 (MQ=255)
                    ggAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGaa                                                                                                                             <  1:2998321/150‑1 (MQ=255)
                       aCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGACCGc                                                                                                                          >  2:1435771/1‑150 (MQ=255)
                        ccATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGTTAGGTAGACACCGTATCAGACACGGAATCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCTTTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGcc                                                                                                                         <  1:2502769/150‑1 (MQ=255)
                        ccATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGcc                                                                                                                         >  1:1387837/1‑150 (MQ=255)
                         catcatCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCt                                                                                                                        <  1:2501270/150‑1 (MQ=255)
                         catcatCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCt                                                                                                                        >  1:3401962/1‑150 (MQ=255)
                            catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                     <  1:2512586/150‑1 (MQ=255)
                            catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                     <  2:1600954/150‑1 (MQ=255)
                            catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                     <  1:460726/150‑1 (MQ=255)
                            catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                     >  2:813494/1‑150 (MQ=255)
                             atcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtg                                                                                                                    >  2:1858001/1‑150 (MQ=255)
                               cgcgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtggt                                                                                                                  <  2:2542648/150‑1 (MQ=255)
                                  ggCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                               >  1:2896044/1‑150 (MQ=255)
                                  ggCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                               >  1:1543752/1‑150 (MQ=255)
                                     tCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaa                                                                                                            >  2:531346/1‑150 (MQ=255)
                                          agTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGc                                                                                                       >  2:4394788/1‑150 (MQ=255)
                                            ttGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAc                                                                                                     >  1:276999/1‑150 (MQ=255)
                                            ttGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAc                                                                                                     >  1:2809013/1‑150 (MQ=255)
                                             tGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCACAACAGCAcc                                                                                                    >  2:2562695/1‑150 (MQ=255)
                                               gCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTg                                                                                                  <  2:2291573/150‑1 (MQ=255)
                                               gCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTg                                                                                                  <  1:2839164/150‑1 (MQ=255)
                                                cTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTgg                                                                                                 <  2:4021699/150‑1 (MQ=255)
                                                  gTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCg                                                                                               >  2:330339/1‑150 (MQ=255)
                                                       aCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTAc                                                                                          >  2:3034695/1‑150 (MQ=255)
                                                       aCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTAc                                                                                          <  1:1116309/150‑1 (MQ=255)
                                                        cAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACt                                                                                         <  2:2928179/150‑1 (MQ=255)
                                                          gATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCa                                                                                       <  2:3160212/150‑1 (MQ=255)
                                                           aTAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                      >  2:176446/1‑150 (MQ=255)
                                                            tAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACg                                                                                     >  1:2256849/1‑150 (MQ=255)
                                                             aGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACgg                                                                                    >  1:215686/1‑150 (MQ=255)
                                                             aGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACgg                                                                                    >  2:4064081/1‑150 (MQ=255)
                                                              ggTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGc                                                                                   >  1:3858258/1‑150 (MQ=255)
                                                               gTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTGGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCa                                                                                  <  2:2272704/150‑1 (MQ=255)
                                                                 aGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAAt                                                                                <  1:1070872/150‑1 (MQ=255)
                                                                  gACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                               <  2:1690676/150‑1 (MQ=255)
                                                                  gACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                               <  2:807569/150‑1 (MQ=255)
                                                                   acacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCa                                                                              <  2:4610395/150‑1 (MQ=255)
                                                                    cacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCTTAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCAc                                                                             <  2:2695106/150‑1 (MQ=255)
                                                                    cacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCAc                                                                             <  2:2694923/150‑1 (MQ=255)
                                                                     acCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACg                                                                            >  2:3886390/1‑150 (MQ=255)
                                                                       cGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTc                                                                          <  2:3351430/150‑1 (MQ=255)
                                                                       cGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTc                                                                          <  2:1509785/150‑1 (MQ=255)
                                                                          aTCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCAtt                                                                       <  1:1656298/150‑1 (MQ=255)
                                                                             aGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATg                                                                    <  1:1865196/150‑1 (MQ=255)
                                                                                cacGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaa                                                                 >  2:4644916/1‑150 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                              >  1:1359407/1‑150 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                              <  2:1153130/150‑1 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                              <  1:125831/150‑1 (MQ=255)
                                                                                    gAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACg                                                             <  1:3847194/150‑1 (MQ=255)
                                                                                    gAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACg                                                             >  2:951397/1‑150 (MQ=255)
                                                                                     aaaCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGc                                                            >  2:3087217/1‑150 (MQ=255)
                                                                                       aCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                          <  2:11084/150‑1 (MQ=255)
                                                                                       aCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                          <  1:3561187/150‑1 (MQ=255)
                                                                                            aGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTc                                                     <  2:214304/150‑1 (MQ=255)
                                                                                              ttGCTGTAACTTCTATCGCAATAAAACCTTTGACCCGCGCACGAAGTTTTGTCTGTAACTGTGACTGGCGAGCTAACGCCTTTTGGTGGCAAAACAGCACGTGGCGGTTACTCACGGCAATCACGTCTTTTTTAAAAACGCCCTGGTCGa                                                   <  1:866762/150‑1 (MQ=255)
                                                                                              ttGCCGTAAATTCTATCGCAATAAAACCATTGTCCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGAACCTGGCGGTTACTCACGGATATCTCGTCATTATGAAAAACGCCCTGGTCGa                                                   <  1:731178/150‑1 (MQ=255)
                                                                                              ttGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGa                                                   <  1:3137001/150‑1 (MQ=255)
                                                                                              ttGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGa                                                   >  1:4023891/1‑150 (MQ=255)
                                                                                              ttGCCGGAACTTCTATCGCCATAAAACCATAGACCCGCGCACTCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGa                                                   >  1:4290476/1‑150 (MQ=255)
                                                                                                gCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATa                                                 <  1:3431023/150‑1 (MQ=255)
                                                                                                gCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATa                                                 <  1:2752816/150‑1 (MQ=255)
                                                                                                  cGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAAc                                               >  1:2365764/1‑150 (MQ=255)
                                                                                                   ggAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                              >  2:2348983/1‑150 (MQ=255)
                                                                                                    gAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGt                                             >  2:1035864/1‑150 (MQ=255)
                                                                                                      aCTTCTATCGCCATAAAACCATTGACCCGCGCACTCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTcc                                           >  1:1926659/1‑150 (MQ=255)
                                                                                                      aCTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTcc                                           >  2:4541094/1‑150 (MQ=255)
                                                                                                      aCTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTcc                                           <  2:4258537/150‑1 (MQ=255)
                                                                                                       cTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCg                                          <  2:2430259/150‑1 (MQ=255)
                                                                                                        ttCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCgg                                         <  2:3496188/150‑1 (MQ=255)
                                                                                                           tATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtt                                      >  1:2153772/1‑150 (MQ=255)
                                                                                                           tATCGCCATAAAACAAAAGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGATCGAACGCATGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtt                                      >  2:1891740/1‑150 (MQ=255)
                                                                                                             tCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtttt                                    <  2:919866/150‑1 (MQ=255)
                                                                                                              cGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTg                                   >  1:3269820/1‑150 (MQ=255)
                                                                                                               gCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGc                                  >  1:2560551/1‑150 (MQ=255)
                                                                                                                  aTAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTgg                               >  2:3970168/1‑150 (MQ=255)
                                                                                                                    aaaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGc                             >  2:3426422/1‑150 (MQ=255)
                                                                                                                     aaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                            <  2:280326/150‑1 (MQ=255)
                                                                                                                     aaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                            <  2:3799479/150‑1 (MQ=255)
                                                                                                                     aaaCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                            >  1:121093/1‑150 (MQ=255)
                                                                                                                         cATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                        <  2:163057/150‑1 (MQ=255)
                                                                                                                            tGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCa                     >  2:3514927/1‑150 (MQ=255)
                                                                                                                            tGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCa                     >  1:713498/1‑150 (MQ=255)
                                                                                                                            tGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCa                     >  2:2073989/1‑150 (MQ=255)
                                                                                                                            tGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCa                     <  1:4467734/150‑1 (MQ=255)
                                                                                                                                   cgcacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTg              >  2:1531580/1‑150 (MQ=255)
                                                                                                                                    gcacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTgg             >  2:2030164/1‑150 (MQ=255)
                                                                                                                                          aGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTc       >  1:4303200/1‑150 (MQ=255)
                                                                                                                                          aGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTc       >  2:3606961/1‑150 (MQ=255)
                                                                                                                                            gTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAc     >  1:3922129/1‑150 (MQ=255)
                                                                                                                                             tttGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAcc    <  1:2378614/150‑1 (MQ=255)
                                                                                                                                               tGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTg  >  2:2846631/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTG  >  CP009273/1821549‑1821841

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