Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,230,687 C→A 8.9% L22M (CTG→ATG) ‡ yohJ → UPF0299 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,230,6870CA 8.9% 310.6 / 13.4 111L22M (CTG→ATG) ‡yohJUPF0299 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (45/55);  new base A (8/3);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 115 total aligned reads displayed.

AATAAGCGCTGATTTCTTAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGATTCGCTATA  >  NC_000913/2230539‑2230834
                                                                                                                                                    |                                                                                                                                                   
aaTAAGCGCTGATTTCTTAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACAATGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATgcctgcct                                                                                                                                                    >  2:1039331/1‑150 (MQ=255)
      cgcTGATTTCTTAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGACTTATGCCTGCCTGTATGc                                                                                                                                              >  2:521201/1‑150 (MQ=255)
        cTGATTTCTTAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAg                                                                                                                                            <  1:895409/150‑1 (MQ=255)
            tttCTTAATGTGATCGGTAGCACGTTTTACGAATTTATTGGATGATGAATCCATCTCCTATGGGGTGTTGATTTTGTGTTAGCCACTGAACATTATCTGGGAATTTTTACGCGCTTTGGGGCTGATTCATGCCTGGCGGTATGCAgggct                                                                                                                                        >  2:699353/1‑147 (MQ=255)
                ttAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGCCACTGAACATTATCTGGCAATATTTACAGGCTTTCGTCCTGATTTATGACTGCATGTATGCAGGCATTTTc                                                                                                                                    >  2:541128/1‑150 (MQ=255)
                ttAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTc                                                                                                                                    >  2:337855/1‑150 (MQ=255)
                   aTGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCAtt                                                                                                                                 <  1:427959/150‑1 (MQ=255)
                   aTGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCATACAGGCAGGCAtaaatcag                                                                                                                                 <  2:30691/150‑8 (MQ=255)
                     gtgATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGc                                                                                                                               >  1:1084695/1‑150 (MQ=255)
                         tCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTcc                                                                                                                           <  2:474059/150‑1 (MQ=255)
                         tCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTcc                                                                                                                           <  2:1157276/150‑1 (MQ=255)
                           ggTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCt                                                                                                                         >  2:1170604/1‑150 (MQ=255)
                            gTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGTCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTg                                                                                                                        >  2:324873/1‑150 (MQ=255)
                            gTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTg                                                                                                                        <  1:180899/150‑1 (MQ=255)
                            gTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTg                                                                                                                        >  1:780880/1‑150 (MQ=255)
                             tAGCACGTTTTACGAATTAATTGTATGATGAATCCCTCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGc                                                                                                                       <  2:626821/150‑1 (MQ=255)
                             tAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGc                                                                                                                       <  2:1516073/150‑1 (MQ=255)
                             tAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGc                                                                                                                       <  2:996584/150‑1 (MQ=255)
                               gCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTACTGATTTATGCCTGCATGTATGCAGGCATTTTCATTGCTTCCCTGCTa                                                                                                                     >  1:62028/1‑150 (MQ=255)
                                 aCGTTTTACGAATTAATTGTATGATGAATCCATCGCATCTGGGGTGTTGTTTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCTGGCATTTTCATTGCTTCCCTGCTAcc                                                                                                                   <  1:1322761/150‑1 (MQ=255)
                                  cGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCg                                                                                                                  >  2:246528/1‑150 (MQ=255)
                                  cGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCATACAGGCAGGCAtaaatca                                                                                                                                  <  1:894422/134‑7 (MQ=255)
                                   gTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCgg                                                                                                                 <  2:801943/150‑1 (MQ=255)
                                    ttttACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGt                                                                                                                <  1:1419868/150‑1 (MQ=255)
                                       tACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAAc                                                                                                             <  2:483027/150‑1 (MQ=255)
                                       tACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAAc                                                                                                             <  1:358934/150‑1 (MQ=255)
                                        aCGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAAcc                                                                                                            >  2:718976/1‑150 (MQ=255)
                                        aCGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAAcc                                                                                                            >  2:357442/1‑150 (MQ=255)
                                         cGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCa                                                                                                           >  2:728647/1‑150 (MQ=255)
                                              taattGTATTATGAATACATCTCATCTGGTTTTTTGATTATGAGTAAGACACTTAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCg                                                                                                      <  2:997167/150‑1 (MQ=255)
                                              taattGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCg                                                                                                      >  1:949276/1‑150 (MQ=255)
                                              taattGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCg                                                                                                      >  1:374761/1‑150 (MQ=255)
                                                 ttGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGgc                                                                                                   <  2:780880/150‑1 (MQ=255)
                                                  tGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACTCGCTTTAGTCCTGATTTATGCCTGCCTTTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGgca                                                                                                  >  1:88043/1‑150 (MQ=255)
                                                   gTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATTCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGgcac                                                                                                 >  1:287083/1‑149 (MQ=255)
                                                      tgatgaATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGcat                                                                                              >  1:690634/1‑150 (MQ=255)
                                                           aaTCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCg                                                                                         >  2:450155/1‑150 (MQ=255)
                                                           aaTCCATCTCATCTGGGGTGTGGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCg                                                                                         <  2:663610/150‑1 (MQ=255)
                                                            aTCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCgg                                                                                        <  1:337855/150‑1 (MQ=255)
                                                             tCCATCTCATCTGGGGTGTTGATTATGTGTAAGACCCTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCggg                                                                                       <  2:494863/150‑1 (MQ=255)
                                                               catctcatctGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGAt                                                                                     >  1:779904/1‑150 (MQ=255)
                                                                      actGGGGTGATGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCATTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATc                                                                              <  1:280155/149‑1 (MQ=255)
                                                                       ctGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTTGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATcc                                                                             <  2:822611/150‑1 (MQ=255)
                                                                         ggggTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTg                                                                           >  1:990385/1‑150 (MQ=255)
                                                                           ggTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGtt                                                                         >  2:1222678/1‑150 (MQ=255)
                                                                           ggTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGtt                                                                         >  1:856651/1‑150 (MQ=255)
                                                                           ggTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGtt                                                                         >  1:883245/1‑150 (MQ=255)
                                                                            gtgtTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGttt                                                                        <  2:416235/150‑1 (MQ=255)
                                                                              gtTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGt                                                                      >  1:1515630/1‑150 (MQ=255)
                                                                               ttGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTc                                                                     <  1:730499/150‑1 (MQ=255)
                                                                                     aTGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCctgctg                                                               >  2:234892/1‑150 (MQ=255)
                                                                                         ttAAGACACTGAACATTATCTTGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCt                                                           <  2:62028/149‑1 (MQ=255)
                                                                                            aGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGc                                                        <  2:1652005/150‑1 (MQ=255)
                                                                                             gACACTTAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCa                                                       <  2:949276/150‑1 (MQ=255)
                                                                                              acacTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGATGGCCTTGCaa                                                      >  1:303372/1‑150 (MQ=255)
                                                                                               cacTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCaaa                                                     >  1:1208619/1‑150 (MQ=255)
                                                                                                acTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAAt                                                    <  2:883245/150‑1 (MQ=255)
                                                                                                   gAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAAttct                                                 >  2:685254/1‑150 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTCTGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCgg                                            <  1:240470/150‑1 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCgg                                            <  1:1039331/150‑1 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCgg                                            <  1:1588501/150‑1 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCgg                                            <  1:685254/150‑1 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCgg                                            <  2:88043/150‑1 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCgg                                            <  2:990385/150‑1 (MQ=255)
                                                                                                        ttATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTAGCCTTGCAAATTCTTCCgg                                            <  1:246528/150‑1 (MQ=255)
                                                                                                          aTCTGGCAATATTAACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCa                                          <  2:201214/150‑1 (MQ=255)
                                                                                                             tGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCAATCCGGGCAGCATCATCGGGATGCTGATCCTGTTAGTCCTGCTGGCCTTGCAAATTCTTCCGGCaaaa                                       >  2:947945/1‑150 (MQ=255)
                                                                                                               gCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCAATCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATg                                     <  1:139683/150‑1 (MQ=255)
                                                                                                                cAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATgg                                    >  2:409577/1‑150 (MQ=255)
                                                                                                                 aaTATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATggg                                   >  1:1449775/1‑150 (MQ=255)
                                                                                                                  atatTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCAGGCAAAATgggg                                  >  2:1549255/1‑149 (MQ=255)
                                                                                                                   tatTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTc                                 <  2:1463502/150‑1 (MQ=255)
                                                                                                                   tatTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTc                                 <  2:1319242/150‑1 (MQ=255)
                                                                                                                   aaattACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTc                                 <  2:755898/147‑1 (MQ=255)
                                                                                                                         cgcgCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCg                           >  1:1030866/1‑150 (MQ=255)
                                                                                                                           cgcTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCggg                         <  2:1669352/150‑1 (MQ=255)
                                                                                                                            gcTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCgggg                        >  2:1329098/1‑150 (MQ=255)
                                                                                                                            gcTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTACCTGCTACCGGTAACCATTCCGGGCAGCATCATCTGTATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCgggg                        <  2:902578/150‑1 (MQ=255)
                                                                                                                              tttCGTCCTGATTTATGTCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCGTCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTg                      <  1:324873/150‑1 (MQ=255)
                                                                                                                              tttCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTg                      >  1:1206982/1‑150 (MQ=255)
                                                                                                                                 cGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTa                   >  1:78712/1‑150 (MQ=255)
                                                                                                                                  gTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTAc                  >  2:1428851/1‑150 (MQ=255)
                                                                                                                                   tCCTGACTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACg                 <  1:521201/150‑1 (MQ=255)
                                                                                                                                    ccTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTTCCGGTAACCATTCCGGGCAGCATCATCGGGCTGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCATAATGGGTCATTCCGGGGTGCTACGt                >  1:1380927/1‑150 (MQ=255)
                                                                                                                                     cTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTa               <  1:1683208/150‑1 (MQ=255)
                                                                                                                                      atgaaaaTGCCTGCATACAGGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTAc              >  1:1330737/7‑150 (MQ=255)
                                                                                                                                       tgaaaaTGCCTGCATACATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACt             >  2:1504211/6‑150 (MQ=255)
                                                                                                                                       gATTTATTCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACt             <  2:287083/150‑1 (MQ=255)
                                                                                                                                       gATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCTGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACt             >  1:110952/1‑150 (MQ=255)
                                                                                                                                       gATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAAGCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCGACGTAcg             >  1:1055301/1‑149 (MQ=255)
                                                                                                                                         aaaaTGCCTGCATACAGGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGa           >  2:561306/4‑150 (MQ=255)
                                                                                                                                           tATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATgg                                    >  1:601463/1‑123 (MQ=255)
                                                                                                                                           tATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGAtt         >  2:534607/1‑150 (MQ=255)
                                                                                                                                           tATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGAtt         <  2:1449775/150‑1 (MQ=255)
                                                                                                                                            aTGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGATTc        <  2:1208619/150‑1 (MQ=255)
                                                                                                                                            aTGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGATTc        <  1:450155/150‑1 (MQ=255)
                                                                                                                                            aTGCCTGCATACAGGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACt             <  1:1552635/145‑1 (MQ=255)
                                                                                                                                             tgcctgccTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGATTCg       >  2:1381643/1‑150 (MQ=255)
                                                                                                                                               cctgcctgTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGATTCGCt     >  2:797465/1‑150 (MQ=255)
                                                                                                                                                  gcctgTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGGCAATCCGGGGTGCTACGTACTGATTCGCtata  <  1:718976/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATAAGCGCTGATTTCTTAATGTGATCGGTAGCACGTTTTACGAATTAATTGTATGATGAATCCATCTCATCTGGGGTGTTGATTATGAGTAAGACACTGAACATTATCTGGCAATATTTACGCGCTTTCGTCCTGATTTATGCCTGCCTGTATGCAGGCATTTTCATTGCTTCCCTGCTACCGGTAACCATTCCGGGCAGCATCATCGGGATGCTGATCCTGTTTGTCCTGCTGGCCTTGCAAATTCTTCCGGCAAAATGGGTCAATCCGGGGTGCTACGTACTGATTCGCTATA  >  NC_000913/2230539‑2230834

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