Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,819 A→T 5.5% intergenic (+68/‑33) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,8190AT 5.5% 906.9 / 29.1 251intergenic (+68/‑33)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base A (124/112);  new base T (4/10);  total (129/122)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-02

Only 100 of 256 total aligned reads displayed.

GTGATAAAGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGTTGGTGAAACCCATTG  >  NC_000913/1689561‑1690080
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                     
gTGATAAAGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTAATCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCACGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                                                                                                                                                                                  >  1:219230/1‑279 (MQ=255)
    tAAAGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                                                                                                                                                                                    <  2:432978/274‑1 (MQ=255)
       aGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATgg                                                                                                                                                                                                                                                >  2:576327/1‑275 (MQ=255)
                                             ttgccGTCGAAGGGGATGCAACGTTGTGGAAAGGTTCTCAGAAGTTACCGCTTAAACCGTGTGAAGCAGCTTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCAACGCCGTTTAGAGCGTGTTTACAACAAGCTGTAATAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAATACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGt                                                                                                                                                                                                          <  4:576323/271‑1 (MQ=255)
                                                         gagATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATGGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                                                                                                                                                                                          <  1:179408/277‑1 (MQ=255)
                                                                                        gTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTa                                                                                                                                                                                                                                                   >  2:345211/1‑191 (MQ=255)
                                                                                                              tCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCa                                                                                                                                                                                                                                                                       <  6:704062/149‑1 (MQ=255)
                                                                                                              tCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCa                                                                                                                                                                                                                                                                       >  5:1102549/1‑149 (MQ=255)
                                                                                                                 gCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTc                                                                                                                                                                                                                                                                    >  7:978628/1‑149 (MQ=255)
                                                                                                                    tttATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGc                                                                                                                                                                                                                                                                 >  7:1464378/1‑149 (MQ=255)
                                                                                                                        ttGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGa                                                                                                                                                                                                                                                             >  7:437951/1‑149 (MQ=255)
                                                                                                                         tGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaa                                                                                                                                                                                                                                                            <  7:1631847/149‑2 (MQ=255)
                                                                                                                            cgccAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggt                                                                                                                                                                                                                                                         <  6:115577/149‑5 (MQ=255)
                                                                                                                            cgccAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggt                                                                                                                                                                                                                                                         <  8:1204443/149‑5 (MQ=255)
                                                                                                                             gccAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGTTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCa                                                                                                                                                                                                                                                        >  8:557232/1‑149 (MQ=255)
                                                                                                                                 aCGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                                                                                                                                                                                    >  6:415779/1‑149 (MQ=255)
                                                                                                                                  cGAATCACCGGTGACTGTCAAAGCCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCATATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTa                                                                                                                                                                                                                                                   <  8:1319409/149‑1 (MQ=255)
                                                                                                                                   gAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                                                                                                                                                                                  >  5:1176692/1‑149 (MQ=255)
                                                                                                                                       cACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                                                                                                                                                                                              >  6:588944/1‑149 (MQ=255)
                                                                                                                                               aCTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                                                                                                                                                                                      <  5:233594/149‑1 (MQ=255)
                                                                                                                                               aCTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                                                                                                                                                                                      <  8:1464378/149‑1 (MQ=255)
                                                                                                                                                cTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                                                                                                                                     <  6:1711549/149‑1 (MQ=255)
                                                                                                                                                 tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATg                                                                                                                                                                                                                                    >  6:302235/1‑149 (MQ=255)
                                                                                                                                                  gTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGa                                                                                                                                                                                                                                   <  7:645422/149‑1 (MQ=255)
                                                                                                                                                   tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGTTTATAACACCTTCAGGCGGCCa                                                                                                                                                                                                                                                                       <  5:1419327/112‑1 (MQ=255)
                                                                                                                                                   tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGaa                                                                                                                                                                                                                                  >  6:225858/1‑149 (MQ=255)
                                                                                                                                                   tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGaa                                                                                                                                                                                                                                  >  8:662598/1‑149 (MQ=255)
                                                                                                                                                     aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                                                                                                                                                                                                >  8:1544972/1‑149 (MQ=255)
                                                                                                                                                     aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                                                                                                                                                                                                >  5:1288672/1‑149 (MQ=255)
                                                                                                                                                     aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                                                                                                                                                                                                >  7:1311234/1‑149 (MQ=255)
                                                                                                                                                          ccACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaa                                                                                                                                                                                                                                                            >  7:1117987/1‑115 (MQ=255)
                                                                                                                                                            aCGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCt                                                                                                                                                                                                                                                               >  6:1199508/1‑111 (MQ=255)
                                                                                                                                                                 cGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGcc                                                                                                                                                                                                                                                                <  2:721192/105‑1 (MQ=255)
                                                                                                                                                                  gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAtt                                                                                                                                                                                                                                       >  7:470293/1‑129 (MQ=255)
                                                                                                                                                                  gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAtt                                                                                                                                                                                                                                       <  8:470293/129‑1 (MQ=255)
                                                                                                                                                                  gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGc                                                                                                                                                                                                                   <  5:480550/149‑1 (MQ=255)
                                                                                                                                                                         gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTTACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGc                                                                                                                                                                                                            >  6:735748/1‑149 (MQ=255)
                                                                                                                                                                         gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                                                                                                                                                                                      <  7:358317/107‑1 (MQ=255)
                                                                                                                                                                         gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGc                                                                                                                                                                                                            >  5:1203966/1‑149 (MQ=255)
                                                                                                                                                                         gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGc                                                                                                                                                                                                            >  7:537147/1‑149 (MQ=255)
                                                                                                                                                                          cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                                                                                                                                                                                               >  3:24933/1‑161 (MQ=255)
                                                                                                                                                                                  caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                                                                                                                                                                                                   >  7:1420642/1‑149 (MQ=255)
                                                                                                                                                                                   aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                                                                                                                                                                                                 >  6:938499/1‑102 (MQ=255)
                                                                                                                                                                                   aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAtt                                                                                                                                                                                                                                       >  6:1538531/1‑112 (MQ=255)
                                                                                                                                                                                   aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                                                                                                                                                                                                  >  8:1533295/1‑149 (MQ=255)
                                                                                                                                                                                   aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAAAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                                                                                                                                                                                                  >  8:384364/1‑149 (MQ=255)
                                                                                                                                                                                    acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                                                                                                                                                                                                 <  6:833794/149‑1 (MQ=255)
                                                                                                                                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCa                                                                                                                                                                                                        >  8:702415/1‑141 (MQ=255)
                                                                                                                                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAgg                                                                                                                                                                            >  4:217447/1‑169 (MQ=255)
                                                                                                                                                                                      aaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                                                                                                                                                                                               <  7:209050/149‑1 (MQ=255)
                                                                                                                                                                                        gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAg                                                                                                                                                                                                                                                                      >  4:708875/1‑76 (MQ=255)
                                                                                                                                                                                        gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTa                                                                                                                                                                                             >  5:90076/1‑149 (MQ=255)
                                                                                                                                                                                         cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc                                                                                                                                                                        <  2:217450/169‑1 (MQ=255)
                                                                                                                                                                                            tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaagaa                                                                                                                                                    <  4:399158/186‑1 (MQ=255)
                                                                                                                                                                                                  cTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCt                                                                                                                                                                                   <  6:1603780/149‑1 (MQ=255)
                                                                                                                                                                                                          aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc                                                                                                                                                                           <  5:415779/149‑1 (MQ=255)
                                                                                                                                                                                                           aaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc                                                                                                                                                                               >  7:1524140/1‑144 (MQ=255)
                                                                                                                                                                                                              aTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg                                                                                                                                                                       >  7:799849/1‑149 (MQ=255)
                                                                                                                                                                                                               ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                                                                                                                      >  5:1576158/1‑149 (MQ=255)
                                                                                                                                                                                                                 aaCNTCTCTTGCTAAGCTTATTAAAGGCTTATAACNCCTTCAGGNNGCCAGTCCGCCTGATTTCATTTTNTGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAAc                                                                                                                                           <  3:64131/174‑1 (MQ=255)
                                                                                                                                                                                                                 aaCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt                                                                                                                                                                    <  6:1454021/149‑1 (MQ=255)
                                                                                                                                                                                                                  aCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATg                                                                                                                                                                   <  7:1026808/149‑1 (MQ=255)
                                                                                                                                                                                                                     tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                                                                                                                                >  5:1016581/1‑149 (MQ=255)
                                                                                                                                                                                                                      ctctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtat                                                                                                                                                               >  7:83352/1‑149 (MQ=255)
                                                                                                                                                                                                                           gCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGc                                                                                                                                                                                                            <  8:299380/99‑1 (MQ=255)
                                                                                                                                                                                                                           gCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGc                                                                                                                                                                                                            >  7:299380/1‑99 (MQ=255)
                                                                                                                                                                                                                             tAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCa                                                                                                                                                        >  5:1233902/1‑149 (MQ=255)
                                                                                                                                                                                                                              aaGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCaa                     >  1:271073/1‑279 (MQ=255)
                                                                                                                                                                                                                                 cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                                                                                                                                                                                               <  1:43424/106‑1 (MQ=255)
                                                                                                                                                                                                                                 cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGc                                                                                                                                                         <  5:812942/144‑1 (MQ=255)
                                                                                                                                                                                                                                 cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAAGAATAAACCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaagaa                                                                                                                                                    <  5:1551406/149‑1 (MQ=255)
                                                                                                                                                                                                                                   tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAagaaga                                                                                                                                                  <  6:1118439/149‑1 (MQ=255)
                                                                                                                                                                                                                                   tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCaa                     >  2:249340/1‑274 (MQ=255)
                                                                                                                                                                                                                                    attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                                                                                                                                 <  5:623559/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                                                                                                                <  5:200201/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                                                                                                                <  6:1233902/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                                                                                                                >  7:61851/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTCTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCAGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                                                                                                                >  7:1151952/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                                                                                                            <  8:1240757/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGGCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                                                                                                            <  8:1158272/149‑1 (MQ=255)
                                                                                                                                                                                                                                           gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAAcc                                                                                                                                          >  8:817695/1‑149 (MQ=255)
                                                                                                                                                                                                                                            cTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAAccc                                                                                                                                         <  8:986888/149‑1 (MQ=255)
                                                                                                                                                                                                                                             ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                                                                                                                                     >  3:196991/1‑56 (MQ=255)
                                                                                                                                                                                                                                             ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCACATGGTATACAGGCAAGAAGATTGAAACCCa                                                                                                                                        >  7:737458/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ataACCCCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTCTGGATCATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                                                                                                                                      >  5:1347820/1‑149 (MQ=255)
                                                                                                                                                                                                                                                taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                                                                                                                     >  8:516488/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                                                                                                                                  <  6:1413381/119‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                                                                                                                                  >  5:1413381/1‑119 (MQ=255)
                                                                                                                                                                                                                                                  atgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCACTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaa                                                                                                                                                       <  5:939522/123‑1 (MQ=255)
                                                                                                                                                                                                                                                  acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCa                                                                                                                                                                     <  2:639664/115‑1 (MQ=255)
                                                                                                                                                                                                                                                  acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATcgcg                                            <  2:157896/236‑1 (MQ=255)
                                                                                                                                                                                                                                                    acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGaa                                                                                                                                 <  7:698251/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGaa                                                                                                                                 <  8:866888/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGaa                                                                                                                                 <  6:1035549/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTCGGCGTAGCCTGGACAGGCGGCGCATGGTATACCGGCAAGATGATTGAAACCCATCTCGTACACATGGTCGCGTAGGCGAACGTGCAACTGAACATGAAAGACCCTGACTCCAACCTGGAAGGGAGTTTTAAAAAATATCATAGCGGCGTATTCAGCAGCAAGTTGCAACTTTTGGTGAAACCCAttc  >  4:300027/1‑273 (MQ=255)
                                                                                                                                                                                                                                                       ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGAc                                                                                                                              >  5:1613338/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCa                                                                                                                                                                     <  4:729266/107‑1 (MQ=255)
                                                                                                                                                                                                                                                           ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtata                                                                                                                                                              <  2:740535/113‑1 (MQ=255)
                                                                                                                                                                                                                                                             cggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTCGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGCCGGCGCATGGTATACAGGCAAGAAGATTGAAAACCATCTCGAAGACATGGTc                                                                                                                        <  8:1062166/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              ggcCAGTCCGCCTGTTTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTcg                                                                                                                       >  7:127136/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                     
GTGATAAAGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGTTGGTGAAACCCATTG  >  NC_000913/1689561‑1690080

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