Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,082 G→T 5.3% Q216H (CAG→CAT aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,0820GT 5.3% 473.3 / 11.0 131Q216H (CAG→CATaroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base G (66/57);  new base T (3/4);  total (69/62)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

Only 100 of 136 total aligned reads displayed.

CAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAA  >  NC_000913/1788812‑1789338
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                
   aTTACTTTGACAGTAAATGGGCTGGGCGGCCAAACCGCGACCGAGTTCCGCGATATGGTGCCCGGTCAGTTTGTTGCTGATTTAATCGGTTGTGGCGCGATTGGCGCACGTACTACCGACAGTCAGAGCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAACAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTAACCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCa                                                                                                                                                                                                                                                           <  4:75408/274‑1 (MQ=255)
                            ggCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                  >  2:389324/1‑275 (MQ=255)
                                    aaccttgACCGAGTTCTTCGATATGGTTGCCGGTCAGTTTATTGCTCATTGTATCATTCGGGGCGGGATTGGCGAGCGTACAACCGAAAGGCTGATCCACCGGGATATGGCTTCGGCACTGCCCTCTCCGGAAGGTTTTAAGAATGGATCCGGTGGCAAGACGTGGATTGTTGTGGATGCTAGCCGCGCAGCCCGCTGCAGCCATTTGTTCCTCTGGCCTGGCAAAAAGGGTCAGATGACCATCTATCAGACCAGCGGCAACCTGTATGGcta                                                                                                                                                                                                                            <  4:150653/267‑3 (MQ=255)
                                              gnattcgtgattatGGTGACCGGTCAGGTTATTGATTATTTAATCAGTTGGGGCTGGATTGTCGCACGTACTACCGAAAGTCAGAGCCACCGCGAAATGGCTTCGGCANTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                                                                                                                                                            <  3:108744/236‑1 (MQ=255)
                                                                   ggTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                                                                                                                                 <  1:95280/237‑1 (MQ=255)
                                                                                   aTTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                                                                                                                                           >  1:67159/1‑275 (MQ=255)
                                                                                   aTTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                                                                                                                                           >  2:229472/1‑275 (MQ=255)
                                                                                   aTTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                                                                                                                                           >  2:248058/1‑275 (MQ=255)
                                                                                               ggggCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGc                                                                                                                                                               >  2:513367/1‑275 (MQ=255)
                                                                                                 ggCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGa                                                                                                                                                             >  2:490498/1‑275 (MQ=255)
                                                                                                           gggcACGTCTTACCGAAAGTCAGATCAACGGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCa                                                                                                                                                    <  4:490497/272‑1 (MQ=255)
                                                                                                                      acCGAAGGTCAGATCAACCGCGAAAGGGGTTCGGCAACGTCCGGTCCGGTAGGTTTTAAAAATGGTACAGATGGAAATACGCGGATGGCTGTGGAGGCTATACGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGGCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGAtt                                                                                                                                         <  4:121381/274‑1 (MQ=255)
                                                                                                                          aaaGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATTATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCt                                                                                                                                    >  2:18032/1‑275 (MQ=255)
                                                                                                                                 gATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                                                                                                                                                                                                                                            >  5:860912/1‑148 (MQ=255)
                                                                                                                                 gATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCTTGCGATACGCTGACCGAGTTTGAGTTACCttgaactc                                                                                                                             >  4:474108/1‑268 (MQ=255)
                                                                                                                                   tCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCAt                                                                                                                                                                                                                                                          >  6:848269/1‑148 (MQ=255)
                                                                                                                                   tCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCAt                                                                                                                                                                                                                                                          >  6:185156/1‑148 (MQ=255)
                                                                                                                                    ccACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                                                                                                                                                                                                                                         <  6:860912/148‑1 (MQ=255)
                                                                                                                                     cACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCATCTGACCAttt                                                                                                                                                                                                                                                        <  8:92131/148‑3 (MQ=255)
                                                                                                                                        cgcgAAATGGTTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt                                                                                                                      <  4:389323/275‑1 (MQ=255)
                                                                                                                                         gcgAAATGGCTTCGGCACTCTCATGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCGATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt                                                                                                                      <  1:45911/274‑1 (MQ=255)
                                                                                                                                           gAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                                                                                                                                                                                                                                  <  7:259477/148‑1 (MQ=255)
                                                                                                                                           gAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                                                                                                                                                                                                                                  <  5:14847/148‑1 (MQ=255)
                                                                                                                                            aaaTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACCTCGGGTGGTGGa                                                                                                                   >  2:352190/1‑274 (MQ=255)
                                                                                                                                                ggCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                 >  8:732482/1‑128 (MQ=255)
                                                                                                                                                ggCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                 <  7:732482/128‑1 (MQ=255)
                                                                                                                                                    tCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAg                                                                                                                                                                                                                                                                  >  6:603740/1‑123 (MQ=255)
                                                                                                                                                    tCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAg                                                                                                                                                                                                                                                                  <  5:603740/123‑1 (MQ=255)
                                                                                                                                                         aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                                                                                                                                                                                                                                    <  6:290939/148‑1 (MQ=255)
                                                                                                                                                              ccTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATgg                                                                                                                                                                                                                               <  5:185156/148‑1 (MQ=255)
                                                                                                                                                                    cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAag                                                                                          >  2:213640/1‑275 (MQ=255)
                                                                                                                                                                     ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatt                                                                                                                                                                                                                        <  6:50025/148‑1 (MQ=255)
                                                                                                                                                                       tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGTGTTTGATTTACCTGAACATCGGGTGGTTGATTTCAGCCACGGTAACCGCCAGAaacag                                                                                       >  4:236491/1‑275 (MQ=255)
                                                                                                                                                                       tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAagcag                                                                                       >  1:110973/1‑275 (MQ=255)
                                                                                                                                                                           ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                                                                                                                                      >  8:864148/1‑112 (MQ=255)
                                                                                                                                                                            tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                 >  7:559221/1‑100 (MQ=255)
                                                                                                                                                                            tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                 <  8:559221/100‑1 (MQ=255)
                                                                                                                                                                               aaaaaTGTGACCGATGGCAATACGCGGATTGCTGTGGAGGCTTTCAGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATTATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCg                                                                               <  4:18032/275‑1 (MQ=255)
                                                                                                                                                                               aaaaaTGGTACCGATGGCAATACGCGGAGTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCg                                                                               <  2:474109/275‑1 (MQ=255)
                                                                                                                                                                                  aaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGt                                                                                                      <  2:170086/249‑1 (MQ=255)
                                                                                                                                                                                   aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                  >  5:480596/1‑124 (MQ=255)
                                                                                                                                                                                   aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                  <  6:480596/124‑1 (MQ=255)
                                                                                                                                                                                    tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATTACCATCTATCAGACCAGCGGCaa                                                                                                                                                                                                                                        >  5:193120/1‑117 (MQ=255)
                                                                                                                                                                                    tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                                                                                                                                                                                                                                        <  6:193120/117‑1 (MQ=255)
                                                                                                                                                                                     ggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGCTGATATCGCGGCAGCCTGCGATCCCCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGTTTTCAGCCACGGTAACTGCCAGAGGCATCCCCGTCGCCAGtt                                                                         >  4:193480/1‑275 (MQ=255)
                                                                                                                                                                                       tACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta                                                                                                                                                                                                                    <  8:14394/134‑1 (MQ=255)
                                                                                                                                                                                            aTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAAcc                                                                                                                                                                                                 >  6:473837/1‑148 (MQ=255)
                                                                                                                                                                                                   tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                                                                                                                                                            >  2:544477/1‑98 (MQ=255)
                                                                                                                                                                                                    aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                                                                                                                                                                                                                                   >  8:304460/1‑90 (MQ=255)
                                                                                                                                                                                                    aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                                                                                                                                                                                                                                   <  7:304460/90‑1 (MQ=255)
                                                                                                                                                                                                     cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                    >  5:92066/1‑136 (MQ=255)
                                                                                                                                                                                                       cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                                                                                                                                                                                                                                        <  4:544476/98‑1 (MQ=255)
                                                                                                                                                                                                         gATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                                                                                                                                      >  7:588544/1‑82 (MQ=255)
                                                                                                                                                                                                         gATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                                                                                                                                      <  8:588544/82‑1 (MQ=255)
                                                                                                                                                                                                          aTTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCg                                                                                                                                                                                   <  6:191203/148‑1 (MQ=255)
                                                                                                                                                                                                           ttGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCga                                                                                                                                                                                  <  5:848269/148‑1 (MQ=255)
                                                                                                                                                                                                                gtgGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                                                                                                                                      <  5:445749/139‑1 (MQ=255)
                                                                                                                                                                                                                  ggATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                 <  6:1125037/62‑1 (MQ=255)
                                                                                                                                                                                                                  ggATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                 >  5:1125037/1‑62 (MQ=255)
                                                                                                                                                                                                                   gATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCAAGGTAACTGCCAGAAGCAGCACCGTCGCCCGTTAGCAGTTTGTGAGGATATTTGTCAGCaaa                                            >  2:118403/1‑274 (MQ=255)
                                                                                                                                                                                                                    aTGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCATCTGACCAtt                                                                                                                                                                                                                                                         >  7:978552/1‑67 (MQ=255)
                                                                                                                                                                                                                     tGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGatat                                                                                                                                                                            <  7:71323/144‑1 (MQ=255)
                                                                                                                                                                                                                        tATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt                                                                                                                                                                                                                     <  7:483713/100‑1 (MQ=255)
                                                                                                                                                                                                                         acccGCGCCGCCCCCGCCGGCCATATGTTCCTCCCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCaa                                     <  4:513366/273‑1 (MQ=255)
                                                                                                                                                                                                                           ccGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg                                                                                                                                                                                                                  <  5:492177/100‑1 (MQ=255)
                                                                                                                                                                                                                           ccGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg                                                                                                                                                                                                                  >  6:492177/1‑100 (MQ=255)
                                                                                                                                                                                                                              cgcAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                  <  6:1037713/81‑1 (MQ=255)
                                                                                                                                                                                                                                      gcgcCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGc                         <  2:412448/274‑1 (MQ=255)
                                                                                                                                                                                                                                      gcgcCAGCCATAAGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                    <  7:963651/103‑1 (MQ=255)
                                                                                                                                                                                                                                      gcgcCAGCCATAAGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                    >  8:963651/1‑103 (MQ=255)
                                                                                                                                                                                                                                            gccaGATGTTCCTCTAGCCAAACAAAAATAGACAAAGGACCAGTTAAAGGAACAGGGGCAACCCGTATGGCAACATTATTAGGCGTGGCGGCAAAAAACCGAATTATCATGCGGATGATATAGCCGCAGGCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTGGTCAGCAAATCCGCAATGGCTCTACGGCGATTGc                   <  4:155581/274‑1 (MQ=255)
                                                                                                                                                                                                                                                 atGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg                                                                                                                                                                                                                  <  6:99477/78‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGa                                                                                                                                                                                               >  8:200545/1‑94 (MQ=255)
                                                                                                                                                                                                                                                    ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGa                                                                                                                                                                                               <  7:200545/94‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGa                                                                                                                                           <  6:856742/146‑1 (MQ=255)
                                                                                                                                                                                                                                                         ctcGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCa                                                                                                                                                                                        >  1:13569/1‑96 (MQ=255)
                                                                                                                                                                                                                                                         ctcGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATcg                                                                                                                                                                          >  8:361876/1‑110 (MQ=255)
                                                                                                                                                                                                                                                            gCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGaaa  >  4:444991/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             ccAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                                                                                                                                      >  7:78691/1‑94 (MQ=255)
                                                                                                                                                                                                                                                             ccAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                                                                                                                                      <  8:78691/94‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCa                                                                                                                                                                     >  5:257792/1‑111 (MQ=255)
                                                                                                                                                                                                                                                             acAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGcc                                                                                                                                                                                    <  3:13569/95‑1 (MQ=255)
                                                                                                                                                                                                                                                              cAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACa                                                                                                                               >  7:1099813/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCa                                                                                                                                                                                        >  5:706497/1‑90 (MQ=255)
                                                                                                                                                                                                                                                               aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGa                                                                                                                                  >  6:296232/1‑144 (MQ=255)
                                                                                                                                                                                                                                                               aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGa                                                                                                                                  <  5:296232/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                gACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat                                                                                                                                                                                                                   <  5:485842/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAcc                                                                                                                                     <  5:578548/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCt                                                                                                                            >  8:636058/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATTGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTg                                                                                                                                            >  5:502769/2‑127 (MQ=255)
                                                                                                                                                                                                                                                                      gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAATCCGAAGTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTg                                                                                                                                            <  6:502769/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatga                                                                                                                                                                               >  6:1024710/2‑92 (MQ=255)
                                                                                                                                                                                                                                                                      gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatga                                                                                                                                                                               <  5:1024710/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatg                                                                                                                                                                                >  6:391235/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                      aaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatg                                                                                                                                                                                <  5:391235/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtgg                                                                                                                       >  5:986174/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       aTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCa                                                                                                                                                    >  7:879740/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                       aTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt                                                                                                                      >  8:288159/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                        tGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggtg                                                                                                                     >  6:633271/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              gATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTa                                                                                                                                       >  6:699402/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                
CAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAA  >  NC_000913/1788812‑1789338

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