Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,282,928 T→A 7.3% I318N (ATC→AAC) ‡ agaS → putative galactosamine‑6‑phosphate deaminase/isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,282,9280TA 7.3% 428.9 / 17.1 137I318N (ATC→AAC) ‡agaSputative galactosamine‑6‑phosphate deaminase/isomerase
Reads supporting (aligned to +/- strand):  ref base T (58/68);  new base A (6/4);  total (64/73)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 146 total aligned reads displayed.

AATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGTG  >  NC_000913/3282781‑3283077
                                                                                                                                                   |                                                                                                                                                     
aaTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcg                                                                                                                                                     >  2:151580/1‑150 (MQ=255)
  tCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgtc                                                                                                                                                   <  2:262780/150‑1 (MQ=255)
   cGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgtcg                                                                                                                                                  <  1:2287491/150‑1 (MQ=255)
   cGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgtcg                                                                                                                                                  <  1:1222647/150‑1 (MQ=255)
   cGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcgtcg                                                                                                                                                  >  1:714801/1‑150 (MQ=255)
    gCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                                 <  2:406599/150‑1 (MQ=255)
      tGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTg                                                                                                                                               <  2:1015646/150‑1 (MQ=255)
        gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                             <  1:1309921/150‑2 (MQ=255)
            aTGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                                                                                         <  2:2193964/150‑1 (MQ=255)
             tGACGTAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTcc                                                                                                                                        <  1:485904/150‑1 (MQ=255)
             tGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTcc                                                                                                                                        <  2:2324319/150‑1 (MQ=255)
               aCGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCcac                                                                                                                                      >  1:256529/1‑150 (MQ=255)
                cGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCcaca                                                                                                                                     <  2:1732089/150‑1 (MQ=255)
                  aaaCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACata                                                                                                                                   <  1:1668074/150‑1 (MQ=255)
                   aaCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACatat                                                                                                                                  <  1:2241717/150‑1 (MQ=255)
                    aCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatc                                                                                                                                 <  1:1367745/150‑1 (MQ=255)
                         ggtggtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCGATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCt                                                                                                                            >  1:1294217/1‑150 (MQ=255)
                         ggtggtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGCCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCt                                                                                                                            >  1:323989/1‑150 (MQ=255)
                          ggggtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAAGTTCGCCGTGACAACCAGTCAATGCTTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTg                                                                                                                           <  2:428275/148‑1 (MQ=255)
                             gtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGATGCCGGTCCACATATCATCCTGCCa                                                                                                                        <  2:690386/150‑1 (MQ=255)
                                 tATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgt                                                                                                                    <  2:1006431/150‑1 (MQ=255)
                                 tATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgt                                                                                                                    >  1:1917720/1‑150 (MQ=255)
                                 tATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgt                                                                                                                    >  1:1496529/1‑150 (MQ=255)
                                  aTTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtc                                                                                                                   <  1:240154/150‑1 (MQ=255)
                                   tttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtca                                                                                                                  >  1:1785089/1‑150 (MQ=255)
                                    ttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcac                                                                                                                 <  2:1462774/150‑1 (MQ=255)
                                     tGTCTCCAGCCACCCTTACACGCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCGACcgtcacg                                                                                                                <  2:154049/150‑1 (MQ=255)
                                      gTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCGGGCTGAACTTCGCCGTGACAACCAGGCAATGCGGGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                                                                                               <  1:632399/150‑1 (MQ=255)
                                           cAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACtt                                                                                                          <  2:1806155/150‑1 (MQ=255)
                                                aCCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGTCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCg                                                                                                     <  2:1122528/150‑1 (MQ=255)
                                                       cacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGa                                                                                              >  1:561174/1‑150 (MQ=255)
                                                         cccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGc                                                                                            >  1:246508/1‑150 (MQ=255)
                                                           cGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAg                                                                                          <  1:447070/150‑1 (MQ=255)
                                                            gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAgg                                                                                         >  2:1954534/1‑150 (MQ=255)
                                                               aGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAt                                                                                      >  2:1467430/1‑150 (MQ=255)
                                                                gTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAtt                                                                                     >  2:1798491/1‑150 (MQ=255)
                                                                 tATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAttt                                                                                    >  2:1615135/1‑150 (MQ=255)
                                                                  aTGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAtttt                                                                                   <  1:1187424/150‑1 (MQ=255)
                                                                   tgatcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGACAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTg                                                                                  >  1:1979438/1‑150 (MQ=255)
                                                                   tgatcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTg                                                                                  <  2:1933885/150‑1 (MQ=255)
                                                                     atcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCt                                                                                >  2:2415070/1‑150 (MQ=255)
                                                                     aacttcatctTCTGGCTGAACTTAGCCGTGACAACAAGGCAATGCTTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCt                                                                                >  2:2093326/3‑150 (MQ=255)
                                                                      tcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGCCGTTGAGCAGGCATTTTGCtt                                                                               >  2:142841/1‑150 (MQ=255)
                                                                       cttgatctGCTGGCGGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAGCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCAGTTTATCGACGTTGAGCAGGCATTTTGCTTc                                                                              <  1:530078/150‑1 (MQ=255)
                                                                        ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTcc                                                                             <  2:1412619/150‑1 (MQ=255)
                                                                        ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTcc                                                                             <  1:2176892/150‑1 (MQ=255)
                                                                        ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTcc                                                                             <  1:2073801/150‑1 (MQ=255)
                                                                         tgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCt                                                                            <  1:2066780/150‑1 (MQ=255)
                                                                          gatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGc                                                                                 >  1:851072/1‑144 (MQ=255)
                                                                          gatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTg                                                                           <  1:2449101/150‑1 (MQ=255)
                                                                             ctgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCTTCTCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACTTTGAGCAGGCATTTTGCTTCCTGATg                                                                        <  2:1399436/150‑1 (MQ=255)
                                                                             ctgctgGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATg                                                                        <  1:1608190/150‑1 (MQ=255)
                                                                               gctggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTa                                                                      >  1:1149979/1‑150 (MQ=255)
                                                                                 tggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACg                                                                    <  2:927634/150‑1 (MQ=255)
                                                                                 tggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACg                                                                    <  1:151580/150‑1 (MQ=255)
                                                                                   gctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGcc                                                                  <  2:555698/150‑1 (MQ=255)
                                                                                    ctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGccc                                                                 >  1:1119056/1‑150 (MQ=255)
                                                                                       aaCTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTCCAACCGTCACGTCACTTTAACGACGTTGAGAAGGCATTTTGCTTCCTGATGTACGCCCAGa                                                              >  2:641311/1‑150 (MQ=255)
                                                                                        aCTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCt                                                                            <  1:1746405/135‑1 (MQ=255)
                                                                                        aCTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGAc                                                             <  2:1443222/150‑1 (MQ=255)
                                                                                         cTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAATCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACg                                                            <  1:2490588/150‑1 (MQ=255)
                                                                                         cTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACg                                                            >  2:1810922/1‑150 (MQ=255)
                                                                                         cTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACg                                                            >  2:204433/1‑150 (MQ=255)
                                                                                          ttCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGt                                                           >  1:1351111/1‑150 (MQ=255)
                                                                                             gCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTg                                                        <  1:775616/150‑1 (MQ=255)
                                                                                              ccGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGc                                                       <  1:1467430/150‑1 (MQ=255)
                                                                                               cGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCa                                                      <  2:1578009/150‑1 (MQ=255)
                                                                                                gTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGAGGTTTGCAc                                                     <  2:735039/150‑1 (MQ=255)
                                                                                                gTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCAc                                                     >  1:386629/1‑150 (MQ=255)
                                                                                                 tGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACt                                                    >  1:789188/1‑150 (MQ=255)
                                                                                                     aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATg                                                >  2:857733/1‑150 (MQ=255)
                                                                                                     aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATg                                                <  1:596181/150‑1 (MQ=255)
                                                                                                     aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATg                                                <  1:2227739/150‑1 (MQ=255)
                                                                                                     aaCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATg                                                <  2:671108/150‑1 (MQ=255)
                                                                                                      aCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGc                                               <  2:1496529/150‑1 (MQ=255)
                                                                                                      aCCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGc                                               <  1:2421500/150‑1 (MQ=255)
                                                                                                                gCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGca                                     >  2:837393/1‑150 (MQ=255)
                                                                                                                gCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGca                                     <  1:1612201/150‑1 (MQ=255)
                                                                                                                 cGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGcac                                    >  1:2314234/1‑150 (MQ=255)
                                                                                                                 cGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGcac                                    >  2:2225033/1‑150 (MQ=255)
                                                                                                                    gtAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATg                                 <  1:2474865/150‑1 (MQ=255)
                                                                                                                     tAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATgg                                >  2:1306829/1‑150 (MQ=255)
                                                                                                                       atcgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGc                              >  1:1001810/1‑150 (MQ=255)
                                                                                                                       atcgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGc                              <  1:440072/150‑1 (MQ=255)
                                                                                                                        tcgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCa                             >  2:1208199/1‑150 (MQ=255)
                                                                                                                         cgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCaa                            >  1:2245798/1‑150 (MQ=255)
                                                                                                                          gccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAAt                           >  1:2234402/1‑150 (MQ=255)
                                                                                                                           ccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATa                          <  1:1810922/150‑1 (MQ=255)
                                                                                                                                gccgcGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCg                     >  2:2251778/1‑150 (MQ=255)
                                                                                                                                gccgcGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCg                     >  1:839192/1‑150 (MQ=255)
                                                                                                                                    ggaccggcagcGACGATGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATa                 >  1:647913/7‑150 (MQ=255)
                                                                                                                                       ccggcagcGACGATGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATAccc              >  1:1813475/4‑150 (MQ=255)
                                                                                                                                       ccggcagcGACGATGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATAccc              >  2:782172/4‑150 (MQ=255)
                                                                                                                                          gcagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCAt           <  2:322372/150‑1 (MQ=255)
                                                                                                                                          gcagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCAt           >  1:1658122/1‑150 (MQ=255)
                                                                                                                                          gcagcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCAt           <  1:1114801/150‑1 (MQ=255)
                                                                                                                                             gcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTAATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAg        >  1:105308/1‑150 (MQ=255)
                                                                                                                                               gACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGcc      <  1:687049/150‑1 (MQ=255)
                                                                                                                                                aCATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATGCGCCGGATACCCCATCAGCCa     >  1:698347/1‑150 (MQ=255)
                                                                                                                                                   tcgtcgCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGTg  <  2:323989/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGTG  >  NC_000913/3282781‑3283077

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: