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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,821,6960AT11.7% 743.3 / 74.7 244L274Q (CTG→CAG) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base A (119/95);  new base T (10/19);  total (129/115)
Fisher's exact test for biased strand distribution p-value = 4.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 244 total aligned reads displayed.

TGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGAT  >  CP009273/1821549‑1821843
                                                                                                                                                   |                                                                                                                                                   
 gAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTa                                                                                                                                                  <  1:563227/150‑1 (MQ=255)
  aGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTaa                                                                                                                                                 <  1:3994876/150‑1 (MQ=255)
    gcagcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACt                                                                                                                                               <  2:3082934/150‑1 (MQ=255)
      agcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACtgt                                                                                                                                             <  1:3660225/150‑1 (MQ=255)
        caCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGa                                                                                                                                           >  2:1155690/1‑150 (MQ=255)
          ccAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACt                                                                                                                                         <  2:2765328/150‑1 (MQ=255)
            aaCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTgg                                                                                                                                       <  1:2339949/150‑1 (MQ=255)
             aCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgc                                                                                                                                      <  1:1382453/150‑1 (MQ=255)
             aCATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgc                                                                                                                                      >  2:2273596/1‑150 (MQ=255)
               atcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgcga                                                                                                                                    >  1:3298720/1‑150 (MQ=255)
               atcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgcga                                                                                                                                    <  1:4499033/150‑1 (MQ=255)
               atcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgcga                                                                                                                                    >  1:3298757/1‑150 (MQ=255)
               atcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgcga                                                                                                                                    <  1:476848/150‑1 (MQ=255)
                  atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAGACCATTGACCCGCGCACGCAGGTTTGCCAGTCACTGTGACTGgcgagcg                                                                                                                                 >  2:2679503/1‑150 (MQ=255)
                         catcatCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCt                                                                                                                          >  1:3523667/1‑150 (MQ=255)
                            catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                       >  2:533535/1‑150 (MQ=255)
                             atcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTtg                                                                                                                      >  2:60919/1‑150 (MQ=255)
                                 cgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTtggtgg                                                                                                                  >  1:1057024/1‑150 (MQ=255)
                                  ggCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                                 >  1:3751891/1‑150 (MQ=255)
                                  ggCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                                 >  1:1468981/1‑150 (MQ=255)
                                  ggCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGc                                                                                                                 >  2:2928676/1‑150 (MQ=255)
                                   gCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCa                                                                                                                <  1:2653697/150‑1 (MQ=255)
                                    cTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaa                                                                                                               <  2:1539948/150‑1 (MQ=255)
                                    cTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaa                                                                                                               >  1:3141738/1‑150 (MQ=255)
                                      cagcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaaa                                                                                                             >  2:743094/1‑150 (MQ=255)
                                       agcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAAc                                                                                                            <  1:3244504/150‑1 (MQ=255)
                                           gTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCa                                                                                                        >  1:1041195/1‑150 (MQ=255)
                                            ttGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAc                                                                                                       <  1:2762415/150‑1 (MQ=255)
                                             tGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAcc                                                                                                      >  2:4363675/1‑150 (MQ=255)
                                             tGGCAGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCAcc                                                                                                      <  2:3468098/150‑1 (MQ=255)
                                              ggCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCt                                                                                                     >  1:842416/1‑150 (MQ=255)
                                                  gTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGGGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCg                                                                                                 >  2:1613041/1‑150 (MQ=255)
                                                     aaaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGtt                                                                                              >  1:4421622/1‑150 (MQ=255)
                                                     aaaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGtt                                                                                              <  2:1695221/150‑1 (MQ=255)
                                                      aaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTa                                                                                             <  2:786385/150‑1 (MQ=255)
                                                       aCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTAc                                                                                            <  2:2332110/150‑1 (MQ=255)
                                                        cAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACt                                                                                           <  2:3783793/150‑1 (MQ=255)
                                                        cAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACt                                                                                           <  1:2504549/150‑1 (MQ=255)
                                                         aGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTc                                                                                          <  1:3255245/150‑1 (MQ=255)
                                                           aTAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                        >  2:4447280/1‑150 (MQ=255)
                                                           aTAGGTAGACACCGTATCAGACACGGAAACCTGAGGTGCCGGACCTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTGGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                        >  1:3054619/1‑150 (MQ=255)
                                                            tAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACg                                                                                       >  2:100892/1‑150 (MQ=255)
                                                                   acacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCa                                                                                >  1:2833282/1‑150 (MQ=255)
                                                                   acacCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCa                                                                                >  1:3641954/1‑150 (MQ=255)
                                                                      ccGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGt                                                                             <  2:4146269/150‑1 (MQ=255)
                                                                       cGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTc                                                                            >  2:2475911/1‑150 (MQ=255)
                                                                        gTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                           >  1:1334823/1‑150 (MQ=255)
                                                                        gTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                           <  2:3875995/150‑1 (MQ=255)
                                                                        gTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                           <  2:2512844/150‑1 (MQ=255)
                                                                        gTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                           >  1:2916195/1‑150 (MQ=255)
                                                                         tATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCAt                                                                          <  1:3636514/150‑1 (MQ=255)
                                                                           tCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTa                                                                        <  1:301765/150‑1 (MQ=255)
                                                                              gACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGa                                                                     >  1:656294/1‑150 (MQ=255)
                                                                                cacGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaa                                                                   <  1:1936108/150‑1 (MQ=255)
                                                                                cacGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaa                                                                   <  1:4210959/150‑1 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                <  2:483113/150‑1 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                <  2:3697710/150‑1 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                >  2:3489564/1‑150 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                >  2:4425298/1‑150 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                >  2:2149607/1‑150 (MQ=255)
                                                                                   ggAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAAc                                                                >  1:2710031/1‑150 (MQ=255)
                                                                                    gAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACg                                                               <  1:1829629/150‑1 (MQ=255)
                                                                                     aaaCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGc                                                              >  1:2842245/1‑150 (MQ=255)
                                                                                      aaCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGcc                                                             <  2:3216119/150‑1 (MQ=255)
                                                                                       aCCTGAGGTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCCCGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAACACGccc                                                            <  2:4009100/150‑1 (MQ=255)
                                                                                           gAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGt                                                        <  2:2140330/150‑1 (MQ=255)
                                                                                           gAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGt                                                        <  2:1399812/150‑1 (MQ=255)
                                                                                           gAGTTGCCGGAACTTATATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCCGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCGGGCGGTTACTCACGGCAATCACGTCTTTATGAAAAACGCCCTGGt                                                        <  1:2181562/150‑1 (MQ=255)
                                                                                            aGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTc                                                       >  1:733545/1‑150 (MQ=255)
                                                                                                 ccGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATaa                                                  <  2:1926190/150‑1 (MQ=255)
                                                                                                  cGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAAc                                                 >  2:218845/1‑150 (MQ=255)
                                                                                                   gtaaCTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGAAAAACAGCACCTGGCGGTTAATCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                                <  1:1240921/148‑1 (MQ=255)
                                                                                                   gtaaCTTATATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATTAAAAACGCCCTGGTCGATAACg                                                <  1:3526416/148‑1 (MQ=255)
                                                                                                   ggAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                                >  2:459640/1‑150 (MQ=255)
                                                                                                   ggAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                                >  1:4159513/1‑150 (MQ=255)
                                                                                                     aaCTTCTATCGCCATAAAACCATTGACCCGCGCACGCTGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTGGGTGGCAAAACAGCACCTGGCGGTTACTCCCGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTc                                              <  1:3944419/150‑1 (MQ=255)
                                                                                                        ttCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCgg                                           >  2:3385577/1‑150 (MQ=255)
                                                                                                              cGCCATAAAACCATTGACCCGCGCACGCTGTTTTGCCATTAACTGTGAGTGGCGAGCGAACGCCTTTTGGTGGAAAAACAGCACCTGGGGTTTACTAACGGAAATCACTTATTTATTAAAAACGCCCTGGTCGATAACGTCCGGGTTTTg                                     <  1:4323666/150‑1 (MQ=255)
                                                                                                                ccATAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCt                                   >  2:2983382/1‑150 (MQ=255)
                                                                                                                   tAAAACCATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTggg                                <  2:3547660/150‑1 (MQ=255)
                                                                                                                    aaaaCCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGc                               >  2:4098652/1‑150 (MQ=255)
                                                                                                                    aaaaCCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGc                               >  1:4244527/1‑150 (MQ=255)
                                                                                                                       aCCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaa                            >  1:3876366/1‑150 (MQ=255)
                                                                                                                        ccATTGACCCGCGCACGCAGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaa                           <  2:2893331/150‑1 (MQ=255)
                                                                                                                        ccATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaa                           >  2:1499868/1‑150 (MQ=255)
                                                                                                                         cATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGaaaa                          >  1:4275067/1‑150 (MQ=255)
                                                                                                                          aTTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAAt                         >  1:2378277/1‑150 (MQ=255)
                                                                                                                          aTTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAAAGCCTGTTGGTGGCAAAACAGCACCTGTCGGTTACTCACTGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAAt                         >  2:1183270/1‑150 (MQ=255)
                                                                                                                                 cgcgCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt                  >  2:2542566/1‑150 (MQ=255)
                                                                                                                                   cgcacgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGGTg                >  2:2861050/1‑150 (MQ=255)
                                                                                                                                    gcacgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTgg               >  1:1134709/1‑150 (MQ=255)
                                                                                                                                     cacgcaGGTTTGCCAGTAACTGTGACTGGCGCGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTggg              <  1:1302193/150‑1 (MQ=255)
                                                                                                                                     cacgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTggg              >  1:2601875/1‑150 (MQ=255)
                                                                                                                                      acgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTgggg             >  1:2986976/1‑150 (MQ=255)
                                                                                                                                        gcaGGTTTGCCTGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGAt           >  1:3205203/1‑150 (MQ=255)
                                                                                                                                          aGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTc         >  1:3113285/1‑150 (MQ=255)
                                                                                                                                           ggTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCa        <  1:759270/150‑1 (MQ=255)
                                                                                                                                            gTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAc       <  2:1657652/150‑1 (MQ=255)
                                                                                                                                              ttGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCt     <  2:4035112/150‑1 (MQ=255)
                                                                                                                                                 ccAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGAt  <  1:2163504/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGAGGCAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGAT  >  CP009273/1821549‑1821843

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: