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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,821,6960AT97.4% 597.2 / 4.6 189L274Q (CTG→CAG) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base T (95/89);  total (99/90)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01
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 189 total aligned reads displayed.

AGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACC  >  CP009273/1821551‑1821839
                                                                                                                                                 |                                                                                                                                                
aGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGt                                                                                                                                                >  1:1154670/1‑148 (MQ=255)
aGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGt                                                                                                                                                >  1:1156436/1‑148 (MQ=255)
    agcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACt                                                                                                                                            <  2:3186472/148‑1 (MQ=255)
      caCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACtgt                                                                                                                                          >  1:2075010/1‑148 (MQ=255)
           aCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTg                                                                                                                                     >  2:101634/1‑148 (MQ=255)
                atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgag                                                                                                                                >  2:2448416/1‑148 (MQ=255)
                 tGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGgcgagc                                                                                                                               >  2:3554352/1‑148 (MQ=255)
                    aaCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGaa                                                                                                                            <  2:3431129/148‑1 (MQ=255)
                     aCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAAc                                                                                                                           >  2:898639/1‑148 (MQ=255)
                      ccATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACg                                                                                                                          >  2:415302/1‑148 (MQ=255)
                        atcatcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGcc                                                                                                                        <  1:962963/148‑1 (MQ=255)
                        atcatcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGcc                                                                                                                        >  2:2471016/1‑148 (MQ=255)
                        atcatcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGcc                                                                                                                        <  1:1620751/148‑1 (MQ=255)
                         tcatcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGACCGCCt                                                                                                                       >  1:3445182/1‑148 (MQ=255)
                         tcatcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCt                                                                                                                       >  2:1858596/1‑148 (MQ=255)
                          catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTg                                                                                                                      >  2:3344079/1‑148 (MQ=255)
                          catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTg                                                                                                                      >  2:2347992/1‑148 (MQ=255)
                          catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTg                                                                                                                      <  2:1751378/148‑1 (MQ=255)
                           atcGCGGCTCAGCATTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACTGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGt                                                                                                                     <  1:2183230/148‑1 (MQ=255)
                            tcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                    >  2:2690972/1‑148 (MQ=255)
                            tcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                    >  1:1960377/1‑148 (MQ=255)
                             cgcgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtg                                                                                                                   >  1:3160043/1‑148 (MQ=255)
                             cgcgGCTCAGCAGTTGGCTGTTAAACAGAGAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtg                                                                                                                   <  2:459039/148‑1 (MQ=255)
                              gcgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtgg                                                                                                                  >  2:1916057/1‑148 (MQ=255)
                               cgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtggt                                                                                                                 >  1:1050197/1‑148 (MQ=255)
                                  cTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                              >  2:2783096/1‑148 (MQ=255)
                                    cagcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaa                                                                                                            >  2:2415218/1‑148 (MQ=255)
                                          ttGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGc                                                                                                     <  1:2281541/148‑1 (MQ=255)
                                            ggCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAc                                                                                                   <  2:359132/148‑1 (MQ=255)
                                              cTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCt                                                                                                 >  2:1204675/1‑148 (MQ=255)
                                                 ttAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGc                                                                                              >  1:1978338/1‑148 (MQ=255)
                                                  tAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCg                                                                                             >  1:503414/1‑148 (MQ=255)
                                                    aaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGt                                                                                           <  1:2285264/148‑1 (MQ=255)
                                                    aaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGt                                                                                           >  2:287706/1‑148 (MQ=255)
                                                      cAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTa                                                                                         >  2:1345643/1‑148 (MQ=255)
                                                        gATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACt                                                                                       >  2:3416479/1‑148 (MQ=255)
                                                           aGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                    <  2:1155351/148‑1 (MQ=255)
                                                            ggTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACg                                                                                   >  2:1873073/1‑148 (MQ=255)
                                                               aGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCa                                                                                <  2:2153570/148‑1 (MQ=255)
                                                                gACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCaa                                                                               >  2:312355/1‑148 (MQ=255)
                                                                 acacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAAt                                                                              <  1:1712783/148‑1 (MQ=255)
                                                                  cacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                             >  1:17454/1‑148 (MQ=255)
                                                                  cacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                             >  1:3321333/1‑148 (MQ=255)
                                                                   acCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCa                                                                            >  1:2983483/1‑148 (MQ=255)
                                                                     cGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACg                                                                          <  2:1051837/148‑1 (MQ=255)
                                                                     cGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACg                                                                          <  2:1747546/148‑1 (MQ=255)
                                                                        aTCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                       <  2:1435547/148‑1 (MQ=255)
                                                                        aTCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                       <  1:926570/148‑1 (MQ=255)
                                                                        aTCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                       <  1:2793562/148‑1 (MQ=255)
                                                                          cAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGGCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCAtt                                                                     >  1:2816370/1‑148 (MQ=255)
                                                                          cAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCAtt                                                                     >  2:2507273/1‑148 (MQ=255)
                                                                           aGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTa                                                                    <  1:2643093/148‑1 (MQ=255)
                                                                            gACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTAt                                                                   >  2:553122/1‑148 (MQ=255)
                                                                             acacGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATg                                                                  >  1:3548485/1‑148 (MQ=255)
                                                                               acGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaa                                                                >  1:2906775/1‑148 (MQ=255)
                                                                               acGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTCTGACTGGCGAGCGAAAGCCTGTTGTTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaa                                                                >  2:299384/1‑148 (MQ=255)
                                                                                cGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaa                                                               <  1:1496237/148‑1 (MQ=255)
                                                                                cGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaa                                                               >  1:3309467/1‑148 (MQ=255)
                                                                                  gAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaaaa                                                             <  1:1135488/148‑1 (MQ=255)
                                                                                    aaCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACg                                                           <  1:828147/148‑1 (MQ=255)
                                                                                       cTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                        <  2:2091133/148‑1 (MQ=255)
                                                                                       cTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                        >  1:58360/1‑148 (MQ=255)
                                                                                         gAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCCTTATGAAAAACGCCCTg                                                      >  2:2024354/1‑148 (MQ=255)
                                                                                         gAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTg                                                      >  1:3328685/1‑148 (MQ=255)
                                                                                              gCCGTAAGTTCTATCGCCATAAAACAATTGACCCGCGCACGCTGTTTTGCCTGTAACTGTGACTGGCGAGCTAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTACTCACGGAAATCACGTAATTATGAAAAACGCCCTGGTCGa                                                 <  1:2118706/148‑1 (MQ=255)
                                                                                              gCCGGAACTTCTATCGCCAAAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCGGGGGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGa                                                 <  1:63077/148‑1 (MQ=255)
                                                                                              gCCGGAACTTCTATCCCCATAAAACCATGGACCCGCGCACGCAGGTTTGCCTGTAACTGTGAGTGGCGAGCGAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGa                                                 <  1:2118669/148‑1 (MQ=255)
                                                                                               ccGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGAt                                                >  1:30680/1‑148 (MQ=255)
                                                                                                cGGAACTTCTATCGCCATAAAACCATTGACCCGCGCCCGCAGTTTTGCCTTTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTACTCACGGAAATCACGTCATTATTAAAAACGCCCTGGTCGATa                                               <  1:3362666/148‑1 (MQ=255)
                                                                                                    aCTTCTATCGCAATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGt                                           <  1:2836830/148‑1 (MQ=255)
                                                                                                      ttCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTcc                                         >  1:2217027/1‑148 (MQ=255)
                                                                                                       tCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCg                                        <  2:65936/148‑1 (MQ=255)
                                                                                                         tATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCggg                                      <  2:3225175/148‑1 (MQ=255)
                                                                                                          aTCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGt                                     >  2:1040070/1‑148 (MQ=255)
                                                                                                               cATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCCCGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGc                                <  2:2075010/148‑1 (MQ=255)
                                                                                                                aTAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCt                               >  1:3362282/1‑148 (MQ=255)
                                                                                                                 tAAAACCTTTGACCCGCGCACGATGTTTTGCCTGTAACTTTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCCCCTGGCGGTTACTCACGGAAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTg                              <  2:629243/148‑1 (MQ=255)
                                                                                                                     aCCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCg                          >  1:2318874/1‑148 (MQ=255)
                                                                                                                        aTTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaa                       <  1:2348037/148‑1 (MQ=255)
                                                                                                                         ttGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                      >  1:1001456/1‑148 (MQ=255)
                                                                                                                         ttGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                      >  1:1001439/1‑148 (MQ=255)
                                                                                                                         ttGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                      >  1:2865617/1‑148 (MQ=255)
                                                                                                                         tgcccccGCGCACGCAGTTTTGCCTGTACCTGTGAGGGGCGAGCGACCGCCTGTTGGTGGAAAAACCGCACGGGGCGTTTACCCCCGGAAATCACTTCATTATCAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                      <  1:1548959/144‑1 (MQ=255)
                                                                                                                                gcgcACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTTTGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTg               <  1:2934272/148‑1 (MQ=255)
                                                                                                                                gcgcACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTg               <  2:2949596/148‑1 (MQ=255)
                                                                                                                                gcgcACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTg               <  2:381889/148‑1 (MQ=255)
                                                                                                                                 cgcacgctGGTTTGCCTGTACCTGTGACTGGCGAGCGACCGCCTGTGGGTGGAAAAACCCCACCTGCCGGTTACCCACGGCAATCACGTCATAATCAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt              <  1:2506539/148‑1 (MQ=255)
                                                                                                                                 cgcacgctGGTTTGCCTGTAACGGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCGGGCGGTTACTCCCGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt              <  2:340120/148‑1 (MQ=255)
                                                                                                                                 cgcacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt              <  2:3309467/148‑1 (MQ=255)
                                                                                                                                 cgcacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt              <  2:692842/148‑1 (MQ=255)
                                                                                                                                 cgcacgcaGGTTTGCCTGAAACTGTGACGGGCGAGCGAACGCCTTTTGGGGCAAAACCAGCACCGGGCGGTTACTCACGGAAATCACGTCATTATTAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt              <  1:3359588/148‑1 (MQ=255)
                                                                                                                                  gcacgctGGTTTGCCTTTAACTGTGACGGGCGAGCGACCGCCTGTTGGGGGAAAAACAGCACCTGGGGTTTACTCACGGCAATCACTAAATTATTAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGtt             <  1:755361/148‑1 (MQ=255)
                                                                                                                                  gcacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGtt             >  2:1660119/1‑148 (MQ=255)
                                                                                                                                  gcacgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGtt             <  1:3235654/148‑1 (MQ=255)
                                                                                                                                     cgcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTggg          >  2:3172462/1‑148 (MQ=255)
                                                                                                                                      gcaGGTTTGCCTTTAAGTGTGACTGGCGAGCGAACGCCTGTGGGTGGAAAAACAGCACCGGGGGGTTACTCACGGAAATCACGTCTTTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTgggg         <  1:2448416/148‑1 (MQ=255)
                                                                                                                                       caGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGa        >  2:3435800/1‑148 (MQ=255)
                                                                                                                                         ggTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGAtt      >  1:1810962/1‑148 (MQ=255)
                                                                                                                                             tGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAcc  >  2:1303392/1‑148 (MQ=255)
                                                                                                                                             tGCCTGTAACTGTGACTGGCGAGCGAACGCCTGATGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAcc  >  2:1271192/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACC  >  CP009273/1821551‑1821839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: