Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,821,489 T→A 5.2% I21N (ATT→AAT)  gmk → guanylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,821,4890TA 5.2% 1711.9 / 57.0 462I21N (ATT→AAT) gmkguanylate kinase
Reads supporting (aligned to +/- strand):  ref base T (214/224);  new base A (11/13);  total (225/237)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 473 total aligned reads displayed.

GCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCG  >  NC_000913/3821344‑3821637
                                                                                                                                                 |                                                                                                                                                    
gCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCa                                                                                                                                                   >  2:944985/1‑149 (MQ=255)
gCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCa                                                                                                                                                   <  3:1117506/149‑1 (MQ=255)
       cTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtgga                                                                                                                                            <  1:1178285/149‑4 (MQ=255)
        tACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTttat                                                                                                                                           >  3:504958/1‑149 (MQ=255)
        tACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtggat                                                                                                                                           <  3:6514/149‑5 (MQ=255)
                          cGTGTATAATGAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGt                                                                                                                         <  3:200668/149‑1 (MQ=255)
                            tgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       <  1:520242/149‑1 (MQ=255)
                                taataaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGtt                                                                                                                           >  3:557249/1‑141 (MQ=255)
                                   taaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGtt                                                                                                                           >  4:893057/1‑138 (MQ=255)
                                   taaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtg                                                                                                                                              >  1:231767/1‑118 (MQ=255)
                                    aaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCa                                                                                                               >  2:1271343/1‑149 (MQ=255)
                                     aGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCttt                                                                                                                                             <  1:130500/118‑1 (MQ=255)
                                       cTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGt                                                                                                            >  2:655071/1‑149 (MQ=255)
                                        tCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGtt                                                                                                           <  2:685772/149‑1 (MQ=255)
                                              gTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtgga                                                                                                                                            <  4:707097/110‑4 (MQ=255)
                                              gTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtgga                                                                                                                                            >  3:707097/1‑107 (MQ=255)
                                                aGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTttatt                                                                                                                                          <  4:269694/110‑1 (MQ=255)
                                                aGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaa                                                                                                                                      <  3:761525/114‑1 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAg                                                                                                                                                  >  1:1352047/1‑100 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCttt                                                                                                                                             <  4:208503/105‑1 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGtt                                                                                                                           <  3:218037/123‑1 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       <  4:1317024/127‑1 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                                 <  1:727712/149‑1 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                                 >  2:278588/1‑149 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGc                                                                                                                                                >  2:1355770/1‑102 (MQ=255)
                                                  gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtgga                                                                                                                                            >  4:352739/1‑103 (MQ=255)
                                                    tttATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaaa                                                                                                                                     >  2:520230/1‑111 (MQ=255)
                                                     ttATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       <  4:624212/124‑1 (MQ=255)
                                                             cTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAg                                                                                                                                                  <  2:757595/89‑1 (MQ=255)
                                                             cTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       >  3:960899/1‑116 (MQ=255)
                                                               aaTCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAgg                                                                                                                                                 <  4:1247338/88‑1 (MQ=255)
                                                                 tCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAAccc                                                                                                                                  >  4:1197008/1‑101 (MQ=255)
                                                                   acaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGtt                                                                                                                           <  4:220627/106‑1 (MQ=255)
                                                                   acaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcg                                                                                >  2:220869/1‑149 (MQ=255)
                                                                     aTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGt                                                                                                                         >  2:323662/1‑106 (MQ=255)
                                                                     aTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGt                                                                                                                         <  1:323662/106‑1 (MQ=255)
                                                                      tACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGt                                                                             >  3:841830/1‑149 (MQ=255)
                                                                       aCGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTa                                                                                                                        <  3:1080384/105‑1 (MQ=255)
                                                                         gAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTtta                                                                                                                                            <  2:583940/83‑1 (MQ=255)
                                                                           aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCt                                                                          >  1:928075/1‑147 (MQ=255)
                                                                           aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTgg                                                                        >  2:153832/1‑149 (MQ=255)
                                                                            aGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                       <  2:440167/149‑1 (MQ=255)
                                                                               tACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                    <  4:841830/149‑1 (MQ=255)
                                                                                  tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCtt                                                                                                                                              <  1:671396/72‑1 (MQ=255)
                                                                                  tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                         <  4:534352/141‑1 (MQ=255)
                                                                                  tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGt                                                                  <  3:338764/148‑1 (MQ=255)
                                                                                  tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                                 >  1:849514/1‑149 (MQ=255)
                                                                                   cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTcc                                                                >  4:1134877/1‑149 (MQ=255)
                                                                                    aTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACAccc                                                                                                                <  1:786376/100‑1 (MQ=255)
                                                                                    aTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCacacac                                                                                             >  3:1190869/1‑119 (MQ=255)
                                                                                    aTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacaccaca                                                                                         <  1:112198/123‑1 (MQ=255)
                                                                                    aTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCa                                                               >  3:1141130/1‑149 (MQ=255)
                                                                                    aTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCa                                                               >  4:43187/1‑149 (MQ=255)
                                                                                        cTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTGTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGt                                                           >  1:460335/1‑149 (MQ=255)
                                                                                         tCAAGGCACGCTTTATATTGTTTCTGCCCCCGGTGGCGCGGGTCAATCCAGCCTGATTCAGGCTTTATTAAAAACCGAGCCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACCCTCCAACCGCGTCCTGATGAAGTCCACGGTg                                                          >  3:904739/1‑149 (MQ=255)
                                                                                         tCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGa                                                                                                                     >  4:583447/1‑90 (MQ=255)
                                                                                         tCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                    >  4:1451636/1‑139 (MQ=255)
                                                                                         tCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTg                                                          <  4:1039056/149‑1 (MQ=255)
                                                                                           aaGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGaa                                                        <  3:1301115/149‑1 (MQ=255)
                                                                                             ggCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATTCAGCCTGATTCAGGCTTTATTaaaaa                                                                                                                                     <  3:1230766/70‑1 (MQ=255)
                                                                                              gCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGtt                                                                                                     <  3:284958/101‑1 (MQ=255)
                                                                                                  gCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                         <  3:253927/125‑1 (MQ=255)
                                                                                                   cTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTAttt                                                >  3:1024508/1‑149 (MQ=255)
                                                                                                    tttATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                       >  1:1001428/1‑125 (MQ=255)
                                                                                                     ttATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCt                                              <  4:295083/149‑1 (MQ=255)
                                                                                                      tataTTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                             >  2:1424808/1‑149 (MQ=255)
                                                                                                       atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCttt                                            <  4:798933/149‑1 (MQ=255)
                                                                                                       atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGCGAACATTATTTCttt                                            <  4:1084827/149‑1 (MQ=255)
                                                                                                         atTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCaa                                                                                                                                <  4:1094665/63‑1 (MQ=255)
                                                                                                              ttCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATc                                     <  1:1424808/149‑1 (MQ=255)
                                                                                                               tCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCa                                    <  4:495915/149‑1 (MQ=255)
                                                                                                                  gCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatga                                 <  1:153832/149‑1 (MQ=255)
                                                                                                                    ccccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcacaca                                                                                              >  3:1334729/1‑86 (MQ=255)
                                                                                                                    ccccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgatg                               >  3:920735/1‑149 (MQ=255)
                                                                                                                     cccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                             <  2:843039/134‑1 (MQ=255)
                                                                                                                      ccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGaa                             >  3:1390699/1‑149 (MQ=255)
                                                                                                                      ccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGaa                             <  2:1109121/149‑1 (MQ=255)
                                                                                                                              gcgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTaaa                       <  3:1274614/147‑1 (MQ=255)
                                                                                                                               cgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGaa                    <  4:130515/149‑1 (MQ=255)
                                                                                                                                gggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                                 <  4:1365822/103‑1 (MQ=255)
                                                                                                                                gggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGaaa                   <  4:920735/149‑1 (MQ=255)
                                                                                                                                gggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGaaa                   >  3:1459026/1‑149 (MQ=255)
                                                                                                                                 ggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTGAAGAAAt                  >  2:105597/1‑149 (MQ=255)
                                                                                                                                  gTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAcc                                                                                 <  3:1092667/85‑1 (MQ=255)
                                                                                                                                  gTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                         <  3:422320/93‑1 (MQ=255)
                                                                                                                                  gTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgat                                >  1:1238587/1‑134 (MQ=255)
                                                                                                                                  gTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAtt                           >  2:895780/1‑139 (MQ=255)
                                                                                                                                   tAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGa                <  1:606769/149‑1 (MQ=255)
                                                                                                                                      aTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                         >  1:279530/1‑89 (MQ=255)
                                                                                                                                       tCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAAt                  >  3:468117/1‑143 (MQ=255)
                                                                                                                                        ccAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatga                                 <  2:1192893/127‑1 (MQ=255)
                                                                                                                                          aGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCag         >  3:599581/1‑149 (MQ=255)
                                                                                                                                            ccTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCGGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTg                                                          >  1:1075774/1‑98 (MQ=255)
                                                                                                                                            ccTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCttt                                            >  4:1146658/1‑112 (MQ=255)
                                                                                                                                            ccTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                             <  3:1305605/111‑1 (MQ=255)
                                                                                                                                             cTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAgaga      >  1:983461/1‑149 (MQ=255)
                                                                                                                                              tGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAAt                                      <  1:391659/116‑1 (MQ=255)
                                                                                                                                              tGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAAt                                      >  2:391659/1‑116 (MQ=255)
                                                                                                                                              tGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCa          <  4:1291005/144‑1 (MQ=255)
                                                                                                                                                 ttCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcaca                                                                                                <  2:1312669/55‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCG  >  NC_000913/3821344‑3821637

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