Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,445 T→A 6.0% Q272L (CAG→CTG)  astB ← N‑succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,4450TA 6.0% 323.6 / 7.8 100Q272L (CAG→CTG) astBN‑succinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (48/46);  new base A (4/2);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 104 total aligned reads displayed.

CAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCTT  >  NC_000913/1827297‑1827594
                                                                                                                                                    |                                                                                                                                                     
cagcagCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAAGAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATTAAAGCATTGACCCGCGCACGCAGGTTTGCCAGTAACtg                                                                                                                                                      <  2:781245/150‑1 (MQ=255)
 agcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACtgt                                                                                                                                                     <  1:48523/150‑1 (MQ=255)
  gcaCCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACtgtg                                                                                                                                                    >  2:711545/1‑150 (MQ=255)
           tcatGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTTgcgag                                                                                                                                           >  2:152125/1‑150 (MQ=255)
           taatGGATCCATCATCGCGGCTCATCATTGGGGTGTTAAACAGATACGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACACTCGAAAGCAGGTTTGCCAGTAACTGTGACTGgcgag                                                                                                                                           <  2:1126005/148‑1 (MQ=255)
             atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTCACAGTTACTGgcaaacc                                                                                                                                         <  2:35265/150‑6 (MQ=255)
             atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTagcgagcg                                                                                                                                         <  2:1582195/150‑1 (MQ=255)
             atGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGgcgagcg                                                                                                                                         <  2:1033181/150‑1 (MQ=255)
               ggAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGaa                                                                                                                                       >  2:559366/1‑150 (MQ=255)
                  aCCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGc                                                                                                                                    <  1:905315/150‑1 (MQ=255)
                   ccATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGcc                                                                                                                                   <  1:1317105/150‑1 (MQ=255)
                     atcatcGCGGCTCAGCAGTTGGCAGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTg                                                                                                                                 <  2:815788/150‑1 (MQ=255)
                       catcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGtt                                                                                                                               >  1:861996/1‑150 (MQ=255)
                         tcGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTtgg                                                                                                                             <  1:258701/150‑1 (MQ=255)
                          cgcgGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTtggt                                                                                                                            >  2:1151757/1‑150 (MQ=255)
                           gcgGCTCAGCAGTTGGCTGTTAAACAGATAGGGAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTtggtg                                                                                                                           <  1:876415/150‑1 (MQ=255)
                                 cagcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaaa                                                                                                                     >  2:76231/1‑150 (MQ=255)
                                 cagcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCACCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaaa                                                                                                                     >  1:438433/1‑150 (MQ=255)
                                 cagcacTTGGCCGTTAAACACCTAGGTAGACACCGTATCACACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCaaaa                                                                                                                     <  2:751825/150‑1 (MQ=255)
                                   gcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGTAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACa                                                                                                                   >  2:491790/1‑150 (MQ=255)
                                   gcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACa                                                                                                                   <  1:633808/150‑1 (MQ=255)
                                   gcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACa                                                                                                                   <  1:575018/150‑1 (MQ=255)
                                   gcagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACa                                                                                                                   <  1:888685/150‑1 (MQ=255)
                                    cagTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAg                                                                                                                  >  2:1015210/1‑150 (MQ=255)
                                     agTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGc                                                                                                                 <  1:1229869/150‑1 (MQ=255)
                                     agTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGc                                                                                                                 <  1:1151757/150‑1 (MQ=255)
                                            tGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGc                                                                                                          >  1:427308/1‑150 (MQ=255)
                                             gTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCg                                                                                                         >  1:187674/1‑150 (MQ=255)
                                             gTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCg                                                                                                         >  1:1113366/1‑150 (MQ=255)
                                              ttAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCgg                                                                                                        >  2:1468634/1‑150 (MQ=255)
                                               tAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGt                                                                                                       >  2:1234179/1‑150 (MQ=255)
                                                aaaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGtt                                                                                                      >  1:358151/1‑150 (MQ=255)
                                                aaaCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGtt                                                                                                      >  2:1275142/1‑150 (MQ=255)
                                                  aCAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCACCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTAc                                                                                                    <  2:438433/150‑1 (MQ=255)
                                                   cAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACt                                                                                                   >  2:1210200/1‑150 (MQ=255)
                                                   cAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACt                                                                                                   >  1:993163/1‑150 (MQ=255)
                                                      aTAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCAc                                                                                                <  2:1339367/150‑1 (MQ=255)
                                                        aGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCCCCTGGCGGTTACTCACgg                                                                                              >  2:1075894/1‑150 (MQ=255)
                                                         ggTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGc                                                                                             >  2:325803/1‑150 (MQ=255)
                                                           tAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCaa                                                                                           >  2:817942/1‑150 (MQ=255)
                                                             gACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                                         >  1:408264/1‑150 (MQ=255)
                                                                acCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACTGCAATCACg                                                                                      <  1:513017/150‑1 (MQ=255)
                                                                acCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACg                                                                                      <  1:559366/150‑1 (MQ=255)
                                                                   gTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCa                                                                                   <  1:1063983/150‑1 (MQ=255)
                                                                      tCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGAGCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTa                                                                                <  2:212820/150‑1 (MQ=255)
                                                                      tCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTa                                                                                <  1:1188739/150‑1 (MQ=255)
                                                                       cAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGtt                                                                                                      <  1:303445/127‑1 (MQ=255)
                                                                       cAGACACGGAAACCTGACTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTAt                                                                               <  1:903686/150‑1 (MQ=255)
                                                                        aGACACGGAAACCAGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATg                                                                              >  2:485704/1‑150 (MQ=255)
                                                                         gACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGa                                                                             >  2:1644362/1‑150 (MQ=255)
                                                                           cacGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaa                                                                           >  1:67928/1‑150 (MQ=255)
                                                                            acGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaaa                                                                          <  1:1234179/150‑1 (MQ=255)
                                                                            acGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaaa                                                                          <  2:408264/150‑1 (MQ=255)
                                                                             cGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGaaaaa                                                                         >  2:82964/1‑150 (MQ=255)
                                                                                  aCCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                                    <  2:427308/150‑1 (MQ=255)
                                                                                   ccTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCt                                                                   <  2:1304858/150‑1 (MQ=255)
                                                                                      gAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGt                                                                <  1:1210200/150‑1 (MQ=255)
                                                                                       aGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTc                                                               <  1:711545/150‑1 (MQ=255)
                                                                                         ttGCCTGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGa                                                             <  1:568702/150‑1 (MQ=255)
                                                                                           gcaggAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATa                                                           >  2:578665/4‑150 (MQ=255)
                                                                                           gCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGCCATTATGAAAAACGCCCTGGTCGATa                                                           >  1:613575/1‑150 (MQ=255)
                                                                                              agAACTTCTATAGCCATAAAACAATTGACCCGAGAACGCAGCTTTGCCAGTAACTGTGACTGACGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                                        >  2:235469/2‑150 (MQ=255)
                                                                                                aaCTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTc                                                      >  1:890211/1‑150 (MQ=255)
                                                                                                    tCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCCCGTCATTATGAAAAACGCCCTGGTCGATAACGTCCggg                                                  >  1:1428523/1‑150 (MQ=255)
                                                                                                    tCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCggg                                                  >  1:1562586/1‑150 (MQ=255)
                                                                                                      tATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtt                                                >  1:277026/1‑150 (MQ=255)
                                                                                                       aTCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGttt                                               <  1:902520/150‑1 (MQ=255)
                                                                                                        tCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGtttt                                              <  1:604468/150‑1 (MQ=255)
                                                                                                          gCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                                         <  2:1621013/105‑1 (MQ=255)
                                                                                                          gCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATc                                                                                         >  1:1621013/1‑105 (MQ=255)
                                                                                                          gCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGc                                            >  1:1479892/1‑150 (MQ=255)
                                                                                                           ccATAAAACCATTGACCCGCGCACGCAGGTTTGCGAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCt                                           <  2:358151/150‑1 (MQ=255)
                                                                                                             aTAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTgg                                         <  1:1306594/150‑1 (MQ=255)
                                                                                                                aaaCCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                                    >  2:472996/1‑120 (MQ=255)
                                                                                                                aaaCCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGccc                                                                    <  1:472996/120‑1 (MQ=255)
                                                                                                                       tGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCAGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCa                               >  2:517098/1‑150 (MQ=255)
                                                                                                                         aCCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCAcc                             >  1:1263215/1‑150 (MQ=255)
                                                                                                                         aCCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCAcc                             <  1:1015210/150‑1 (MQ=255)
                                                                                                                         aCCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCAcc                             <  2:316472/150‑1 (MQ=255)
                                                                                                                          cccGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCt                            >  2:134951/1‑150 (MQ=255)
                                                                                                                            cgcgCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGt                          <  2:245108/150‑1 (MQ=255)
                                                                                                                                cacgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTggg                      <  1:1644362/150‑1 (MQ=255)
                                                                                                                                  cgcaGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGa                    >  2:1009248/1‑150 (MQ=255)
                                                                                                                                     aGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTc                 <  1:1275142/150‑1 (MQ=255)
                                                                                                                                      ggTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCa                >  1:13594/1‑150 (MQ=255)
                                                                                                                                       gctcgCCAGTCACAGTTACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAc               >  2:812657/5‑150 (MQ=255)
                                                                                                                                       gctcgCCAGTCACAGTTACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAc               >  2:784558/5‑150 (MQ=255)
                                                                                                                                       gTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCCCCTGTTGGGGATTCAc               >  2:1505394/1‑150 (MQ=255)
                                                                                                                                        ctcgCCAGTCACAGTTACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                                        <  2:780071/105‑1 (MQ=255)
                                                                                                                                        ctcgCCAGTCACAGTTACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACg                                                        >  1:780071/4‑108 (MQ=255)
                                                                                                                                        ctcgCCAGTCACAGTTACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCAcc              >  1:1414400/4‑150 (MQ=255)
                                                                                                                                           gCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGa           >  1:1311638/1‑150 (MQ=255)
                                                                                                                                           gCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGa           >  2:1649268/1‑150 (MQ=255)
                                                                                                                                             cAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGAtt         >  2:70642/1‑150 (MQ=255)
                                                                                                                                             cAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGAtt         <  1:1149102/150‑1 (MQ=255)
                                                                                                                                              aGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTc        >  1:818115/1‑150 (MQ=255)
                                                                                                                                                tAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAg      >  2:232229/1‑150 (MQ=255)
                                                                                                                                                 aaCTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGc     <  1:6780/150‑1 (MQ=255)
                                                                                                                                                   cTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCt   <  2:67928/150‑1 (MQ=255)
                                                                                                                                                    tgtgACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCtt  >  2:1099781/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCACCAACATCATGGAACCATCATCGCGGCTCAGCAGTTGGCTGTTAAACAGATAGGTAGACACCGTATCAGACACGGAAACCTGAGTTGCCGGAACTTCTATCGCCATAAAACCATTGACCCGCGCACGCAGGTTTGCCAGTAACTGTGACTGGCGAGCGAACGCCTGTTGGTGGCAAAACAGCACCTGGCGGTTACTCACGGCAATCACGTCATTATGAAAAACGCCCTGGTCGATAACGTCCGGGTTTTGCTGGGCGAAAATCACCTGTTGGGGATTCACCTGATTCAGCCTT  >  NC_000913/1827297‑1827594

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