Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,914,648 C→T A185T (GCG→ACG)  glmU ← fused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,914,6480CT92.2% 166.2 / 5.1 66A185T (GCG→ACG) glmUfused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (20/39);  total (22/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

GTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAG  >  NC_000913/3914401‑3914895
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
gatAGCGTTCCACGCAGATCATAACGCCCTGTATCGCTCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCTCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTTAGCTTCgtcag                                                                                                                                                                                                                                                      <  2:208657/249‑1 (MQ=255)
   agcttccctctcataataaaaCGGGCTGGGTCGCTCAGCATAAAGCCTGCTAACGGCAGTTTTTCAGCCTGTTCGGACTGATAACCACGCTCCAGACGGGAGAGGTGCAGGCGGTTATTCACGTCTTCTCCTTCTCTTAAAAGTTGCGGATGAACGGCGTCGATTGTGCGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGACGTGAGCATTATTGTTTGTTAGCGTCgtcagcca                                                                                                                                                                                                                                                   <  2:352317/234‑1 (MQ=255)
           ccGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCa                                                                                                                                                                                                                                           <  1:165012/250‑1 (MQ=255)
              cAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCatat                                                                                                                                                                                                                                        <  1:256919/251‑1 (MQ=255)
                gATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCt                                                                                                                                                                                                                                      <  1:1657/251‑1 (MQ=255)
                  ccAAAACGCGCTGGATCGCGCAACATAACGCCTGCTAACAGCATTTTTTCAGCCTTGTCGGACTGATAAACTCGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAGGGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTgc                                                                                                                                                                                                                                    <  2:557691/250‑1 (MQ=255)
                    aaaaCGCGCTGGATCGCGCAGCAGAACGCCTGCGAACAGCAGTTTTTCATCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTgcgc                                                                                                                                                                                                                                  <  2:292846/251‑1 (MQ=255)
                     aaaCGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGcc                                                                                                                                                                                                                                 <  1:19609/251‑1 (MQ=255)
                         gctcTGGATCGCGCAGCATAACGCCTGCTAACAGTAGGTGTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACTGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTg                                                                                                                                                                                                                             <  2:310076/248‑1 (MQ=255)
                          cgcTGGATTGCGCATCATAACGCCTGCTAAAAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTgg                                                                                                                                                                                                                            <  2:454692/251‑1 (MQ=255)
                                   gcgcAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAGTGGTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGaa                                                                                                                                                                                                                   <  1:60777/251‑1 (MQ=255)
                                    cgcAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt                                                                                                                                                                                                                  <  1:451131/251‑1 (MQ=255)
                                        gcaTAACGCCTGCTACCAGCACTTTGTCAGCCTGTCCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAGCGCTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCTATAATGTCCGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCATCTTCGTCAGCCTGCGTTGCATATCTGCGCCGTTGGCAATCAGAATGcct                                                                                                                                                                                                              <  2:472799/251‑2 (MQ=255)
                                        gcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCg                                                                                                                                                                                                              <  1:291025/251‑1 (MQ=255)
                                          aTAACGCCTGCTACCAGCAGGTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGTGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCGTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGgt                                                                                                                                                                                                            <  2:403462/251‑1 (MQ=255)
                                          aTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTTCGGGTGAACGGCGACGATTTCACGCCCTTCCTGATAAGCCCGCGCCATAATGTCGGTGATGTAGTATTCGCCCCGAGCATTATTTTTGGTCAGCTTTGTCAGCCAGCGGTTCATATCTGCGCCGTTGGCAATCAGAATGCCggg                                                                                                                                                                                                            >  1:401331/1‑250 (MQ=255)
                                          aTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGgt                                                                                                                                                                                                            <  1:443639/251‑1 (MQ=255)
                                                    cTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTTGCTTTAACGTTTCGGGTGAACGGCGACGGTTTCACGCCCTTCCTGATACGCCAGCGCAATAAAGTCGGTGATGTAGTATTCGCCCTGTGGAATTTTGTTGGTCCGCTTCGTCAGCCAGCGGTTCATATTTGCGCCGTTGGCAATCAGAATACAGGTGTTgttatcc                                                                                                                                                                                                  >  1:541354/1‑245 (MQ=255)
                                                     tAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACACCTTCTACTTCGCTTAAACGTTGGGGGTGAACCGCGACGATTTCACCCCCTTCCTGATACGCCAGCGCCATAAAGTCGGGGGTGTTGTATTCGCCCTGGGCATTATTGTTGGTCCGCTTCATTAGCCCGCGTTTTATATCTTCGCCGTTGGCAATCAGAATACCGGGGGTGATCTCCt                                                                                                                                                                                                 >  1:50116/1‑251 (MQ=255)
                                                      aaCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTg                                                                                                                                                                                                <  1:170739/251‑1 (MQ=255)
                                                            agTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTg                                                                                                                                                                                          <  1:429636/251‑1 (MQ=255)
                                                             gttgttcatgttgtttggtcTGATTAACACCCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACACCTATGTCTTCGCTTAAACGTTGCGTATTAACGTCGACGATTTCTCGTCCTTCCTGATACGCCGTCGCAATAATGTCGGTTATGTAGTATTCGCCCTGAGCATTATTGTTTGTCAGCGTCGTCAGCCAGCGTTTCATATCTGCGACGTTGGCAATCAGAATGCCGGCTTTGATCTCCTGAATCTGa                                                                                                                                                                                         <  2:54336/232‑1 (MQ=255)
                                                             gTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGGTGAAAGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCCATAATGTCGGGGATGTAGTATTTGCCCTGAGGATTATTGTTGGGTAGGTTCGTCCGCCCGCGGTTTATATTTGCGCCGGTGGCAATCAGAATGCCGGGGTTGGTCTCCTGAATCTgg                                                                                                                                                                                         >  1:566829/1‑250 (MQ=255)
                                                                    aGCCTGCTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGGCGGTGATGTAGTATTCGCCCTGAGCATTATTGGTGGTCAGCTTCGTCAGCCAGCTTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACgctgct                                                                                                                                                                                  <  1:140159/251‑1 (MQ=255)
                                                                         gTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCtcgtcg                                                                                                                                                                             <  1:317563/251‑1 (MQ=255)
                                                                           tCGGACTGATAAACACGTTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCGGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGt                                                                                                                                                                           <  2:292538/250‑1 (MQ=255)
                                                                             ggACTGCTAAACACGCGCCAGACGGGAGATTTGCAGGCTGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATCCGCCAGCGCAATAATGTCGGTGATGGAGTATTCGCCCTGAGCATTATTTTTGGTCAGCTTCTGCAGCCAGCGTTTCATATCTGCGCCGTTGTCAATCAGAATGTCGGTTTTGATCTCCGGAATCTGACGCTGCTCGCCGGTgg                                                                                                                                                                         <  2:551373/251‑1 (MQ=255)
                                                                                     aaactcgtcccAGATGGGAGAGTTGCTGTCGGTTAGTCACGCTATCTTCTGTGCTTAACCGTTGCGGATGAACTGTGCCGATGTCACGCCCTTCCTGATACCCCAGCTCAATAATTTCGGTGATGTAGTAGTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCATCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTTTTGATTTCCTGAATCTTACGCTGCTCGTCGGTGGCATCtttt                                                                                                                                                                 <  1:346248/242‑2 (MQ=255)
                                                                                     aaaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTtg                                                                                                                                                                 <  1:346259/251‑1 (MQ=255)
                                                                                     aaaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCttt                                                                                                                                                                  <  1:263096/250‑1 (MQ=255)
                                                                                              ccAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt                                                                                                                                                                                                                  <  1:35891/193‑1 (MQ=255)
                                                                                                 gaagagagaggttcAGGCGGTTAATCACGTCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGAAGATTTCACGCCCTTCCAGATAAGCCAGCGCAATAATGTCGGTGATGTAGGATGCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTCTGTGCTCAACAATg                                                                                                                                                     <  2:498276/238‑1 (MQ=255)
                                                                                                    gggAGAGTTGCAGGCGGTTATTTACGCCTTCTTCTTCGCTTAAAAGTTGCGGGTGAAAGGCGACGATTTTACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGGTGTAGTATTTGCCCTGTGCATTTTTGTTGGTCAGCTTCGTCAGCCCGCGTTTTATATCTGCGCCGTTGGCAATCAGAATGTCGGTGGTTATCTCCTTAATCTAACGCTGCTCGTCGGTGGCCTCTTTGTGCGCAACAATTCCg                                                                                                                                                  >  1:546386/1‑251 (MQ=255)
                                                                                                     ggAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGt                                                                                                                                                                           <  1:394034/225‑1 (MQ=255)
                                                                                                            tGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTTGCAATCAGAATGCCGGTGTTGATCGCCTGAATCTGGCGTTGCGCGTCGGTGGCCTCTTTGTGCTCAACAATGCCGGTAACtt                                                                                                                                           <  1:145631/250‑1 (MQ=255)
                                                                                                                ggCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCg                                                                                                                                      <  2:439063/251‑1 (MQ=255)
                                                                                                                   ggTTATTCACGCCTTCTACTTCGCTTAACCGTTGCGGATGAACGGCGACGATTGCACGCCTCTACGGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTTAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTGCATATCTGCGCCGTTGGCAATCAGAATGCCGGTTTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGttt                                                                                                                                   <  2:26615/251‑1 (MQ=255)
                                                                                                                      tATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTc                                                                                                                                 >  1:528565/1‑250 (MQ=255)
                                                                                                                        ttCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGATGGCATCTTTGTGCTCAACAATGCCGGTAACTATGCCGTTTTCACg                                                                                                                              >  1:566146/1‑251 (MQ=255)
                                                                                                                             cccTTCTACCTCGCTTAATCGTTTGGTATGAACGTCTACGAATTCTCGTCCTGCTGGATACGCCAGGGACATTATTTCGGTGATGTAGTATTCGCCCTTAGCATGATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGACGGTGGCAACCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCGGCTTTGTGCTCAACAATGCCCGTAACTTTGCCGTTTTccctggtta                                                                                                                         <  2:188/250‑9 (MQ=255)
                                                                                                                                  tcttcTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGGCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTTAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTCGCCGTTTTCACGGTTGATa                                                                                                                       <  2:479698/245‑1 (MQ=255)
                                                                                                                                       accGCTTAATCGTTGCGTAGGAACGGTGACGATTGCACGCCCTTCCGTAATCGCTAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCTTTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATTTCTGCGCCGTTGGTAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCCCAACAATGCAGGTAAATTTGCCGTTTTCTCGGGTGATACGTCCa                                                                                                                 <  2:393613/248‑1 (MQ=255)
                                                                                                                                             tAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATCCGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCttt                                                                                                                                                                  <  1:298901/195‑1 (MQ=255)
                                                                                                                                             tAAACGTGGCGGATGAGCGGCGACGATTTCTCGCCCCGTCTGAGACGCCAGCGCAATAATGTCGGTAATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTTCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCGCCACAATGCCGGTAAATTTGCCGTTTTCACGGGTGATTCGTCCATAAcc                                                                                                            <  2:197591/249‑1 (MQ=255)
                                                                                                                                              aaaCGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGt                                                                                                         >  1:200606/1‑251 (MQ=255)
                                                                                                                                                acgatGCGGCTGAACGGCGACGATTACACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGTATTATTGTTTGTCAGCTTCGTCAGCCAGCGTTTCATATCAGCACCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGAAGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAGTGCCGGTAACTCTGCCGTTTTCACGGGGGATACGTCCATAACCGGTCg                                                                                                       >  2:85768/5‑251 (MQ=255)
                                                                                                                                                     gCGGATGAACGGGGTCGAGTACACGCCCTTCCTGATACGCCAGTGCAAGAATGTCGGTGATGTAGTATTCGCCCTGAGCATGATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATAGCTTCGCCGTTGGCAATCAGAATGCCGGTGTTGATATCCTGAATCTGACGCTGCTCGTCGGGGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGGCCATAACCGGGCTGatca                                                                                                  <  2:189650/251‑1 (MQ=255)
                                                                                                                                                     gCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGCTGTAGTATTCGCCCTGAGCATTCTTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGatc                                                                                                   >  1:516300/1‑250 (MQ=255)
                                                                                                                                                       ggATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt                                                                                                                                                                                                                  >  1:140169/1‑137 (MQ=255)
                                                                                                                                                       ggATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt                                                                                                                                                                                                                  <  2:140169/137‑1 (MQ=255)
                                                                                                                                                                   acgaTTTCACGCCCTTCCTGATACGCCAGCGCAATATTGTCGGAGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGATGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAAAATTGCCGGGAACTTTGCCGTTTTTCCGGGTGTTCAGTCCATCACCTGTAGGATCATCCCGttttcccgt                                                                                    >  2:382028/1‑245 (MQ=255)
                                                                                                                                                                      aTTTCACGCCCGTCCTGATACGACAGCGCAATAATGTCGGTGATGTAGTATTCGCCATGGTCATTATTGTTGGTCAGCTTCGTCAGTCAGCGTTTCTGCTCTTCGCCGGTGGTAAACGTGATGCTGTTGTTGATCACCTGATTCTGTCGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTACCGTTGCCTTTTTCACTGACGATACGGCCATAACCGGtgtttgctgccagtttccccgtcag                                                                                 <  2:568557/251‑25 (MQ=255)
                                                                                                                                                                          cACGCCCTTCCTGCTACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCCTATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGACCTCCTGAATCTGACCCTGCTCGTCGATGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGGGATACGTCCATCACCGGTCGGAACATCCCGTTTCACCGtacgcaa                                                                              >  2:486721/1‑244 (MQ=255)
                                                                                                                                                                          cACGCCCTGCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTcagcag                                                                              <  1:77427/250‑1 (MQ=255)
                                                                                                                                                                           aCGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGGTGTAGTATTCGCCCTGAGCCTTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGCCGCTGCTCGTCGAGGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTAGTACGTCCATAACCGGCCGGATCATCCCGATTCACCGACAGCAGAc                                                                            >  2:438982/1‑251 (MQ=255)
                                                                                                                                                                                        aTACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTTTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGATAACTTTGCCGTTTTCACGGGTGATACGTCCCTAACCGGTCGGACCATCCCGTTTAACCGTCAGCAGAACAATGCCACCCt                                                                >  1:545590/1‑251 (MQ=255)
                                                                                                                                                                                                        atGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACg                                                >  1:433833/1‑251 (MQ=255)
                                                                                                                                                                                                         tGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGc                                               >  1:306441/1‑251 (MQ=255)
                                                                                                                                                                                                           tCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGc                                               >  1:31327/1‑249 (MQ=255)
                                                                                                                                                                                                                     gttaTTCGCCCTGGGCATTATTGTTGGTCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCTTCAGCAGACCAATGCCGCCCTTCGGTTTAGCATCACGCAGACGCTGGAg                                    <  2:553877/248‑1 (MQ=255)
                                                                                                                                                                                                                                tGAGCATTATTGTTGGTCAGCTTCGCCAGCCATCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGAc                           >  1:204680/1‑248 (MQ=255)
                                                                                                                                                                                                                                       tattGTTGGTCGGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGc                 >  2:168384/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tattGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGACTGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCACCAATGCCGGTAACTTTGCCGTTTTCACGGGGGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCGGAGGTTTAGCATCACGCAGCCGCTGGCGTGTTGCGACagcgagccgc                 >  2:187787/1‑243 (MQ=255)
                                                                                                                                                                                                                                           gttgttCAGCTTCGTCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTTAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATACCGATAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCAGCCAGTTTCGCCGTCAGCAGACCAATGCCACCCTGCGGGTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCa             <  1:217873/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ggTCAGCTTCGTCAACCAGCGTTTCATATCTGCGCCTTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCGCGTCGGTGGCATCTTTGTGCTCACCAATGCCGGTAACTTTGCCGTTTTCACGGGGGATCCGTCCCTAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGCTTTAGCATCACGCAGCCGCTGGAGTGATTCGAccacgttccgcggcacgc          >  2:36163/1‑234 (MQ=255)
                                                                                                                                                                                                                                                gcagtttCGGCAGCCAGCGTTTCATATCTTCGCCGTTGGCAAGCAGAATTCCGGTTTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACTCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCg        <  2:513402/246‑1 (MQ=255)
                                                                                                                                                                                                                                                  atgttCGTCAGCCAGCGTTTCATATCTGCGCCGGTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGcc      <  2:566829/248‑1 (MQ=255)
                                                                                                                                                                                                                                                       cGCCAGCCAGCGTTTCATATCTGCGCCGGTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTag  <  1:187787/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
GTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCATCCAGTTTCACCGTCAGCAGACCAATGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAG  >  NC_000913/3914401‑3914895

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