Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,306 A→T 9.0% intergenic (‑243/+29) rnc ← / ← lepB RNase III/leader peptidase (signal peptidase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,704,3060AT 9.0% 440.6 / 25.0 131intergenic (‑243/+29)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base A (62/57);  new base T (3/9);  total (65/66)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 135 total aligned reads displayed.

CCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGC  >  NC_000913/2704158‑2704453
                                                                                                                                                    |                                                                                                                                                   
ccGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATa                                                                                                                                                     >  1:174205/1‑149 (MQ=255)
 cGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATaa                                                                                                                                                    >  3:540846/1‑149 (MQ=255)
 cGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATaa                                                                                                                                                    >  1:659643/1‑149 (MQ=255)
 cGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCCATACACAACGAAGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATaa                                                                                                                                                    >  2:785395/1‑149 (MQ=255)
     atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCTAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                <  2:22171/149‑1 (MQ=255)
     atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                <  2:134185/149‑1 (MQ=255)
     atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                <  3:470552/149‑1 (MQ=255)
        tGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACCACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                                                                             <  3:333066/149‑1 (MQ=255)
          gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACCCAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACaa                                                                                                                                           <  4:707477/149‑1 (MQ=255)
          gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACaa                                                                                                                                           <  4:381281/149‑1 (MQ=255)
          gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACaa                                                                                                                                           >  4:450019/1‑149 (MQ=255)
           ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                          <  2:740160/149‑1 (MQ=255)
           ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                          <  4:823840/149‑1 (MQ=255)
           ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                          <  2:944463/149‑1 (MQ=255)
           ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTACCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCgcccac                                                                                                                                          <  2:401092/149‑6 (MQ=255)
            gACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACg                                                                                                                                         <  3:899116/149‑1 (MQ=255)
             aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGt                                                                                                                                        <  1:797897/149‑1 (MQ=255)
                aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                                                                     <  4:353099/149‑1 (MQ=255)
                aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                                                                     >  3:72857/1‑149 (MQ=255)
                aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaa                                                                                                                                     >  4:766843/1‑148 (MQ=255)
                 aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                    <  1:76180/149‑1 (MQ=255)
                 aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGAGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGACACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                    <  3:708283/149‑1 (MQ=255)
                   gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                >  1:60790/1‑135 (MQ=255)
                   gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                <  2:60790/135‑1 (MQ=255)
                    cTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGaa                                                                                                                                 <  1:560926/149‑1 (MQ=255)
                     tGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAg                                                                                                                                >  1:339146/1‑149 (MQ=255)
                       tttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAt                                                                                                                              <  3:209626/149‑1 (MQ=255)
                        ttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATg                                                                                                                             <  1:492887/149‑1 (MQ=255)
                           gtgtGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCt                                                                                                                          <  3:1068224/149‑1 (MQ=255)
                            tgtgCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTa                                                                                                                         <  4:167077/149‑1 (MQ=255)
                                 cGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                                                                    <  3:1080938/149‑1 (MQ=255)
                                    ggCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGa                                                                                                                 <  1:346598/149‑1 (MQ=255)
                                     gCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGTGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACACCGTGAACGAATATGGCTATTAATGGAt                                                                                                                >  2:1064457/1‑149 (MQ=255)
                                       aaCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                                                              >  1:966227/1‑149 (MQ=255)
                                       aaCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                                                              >  3:62550/1‑149 (MQ=255)
                                          cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                           >  1:361184/1‑149 (MQ=255)
                                          cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                           <  4:62550/149‑1 (MQ=255)
                                             gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                        >  4:207212/1‑149 (MQ=255)
                                             gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                        <  2:532711/149‑1 (MQ=255)
                                                gTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                     >  1:455160/1‑149 (MQ=255)
                                                   ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                  >  3:283071/1‑149 (MQ=255)
                                                    cAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACt                                                                                                 <  3:763339/149‑1 (MQ=255)
                                                      aCCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                                                                               <  2:695036/149‑1 (MQ=255)
                                                        cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                             >  3:950283/1‑149 (MQ=255)
                                                         ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgc                                                                                            <  3:598216/149‑1 (MQ=255)
                                                         ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgc                                                                                            >  3:43876/1‑149 (MQ=255)
                                                            aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                          >  1:777111/1‑100 (MQ=255)
                                                            aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCa                                                                                         <  3:207212/149‑1 (MQ=255)
                                                              cTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                                                                       <  4:584850/149‑1 (MQ=255)
                                                               ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAc                                                                                      >  1:238973/1‑149 (MQ=255)
                                                                    gggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                     <  2:970074/129‑1 (MQ=255)
                                                                     ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCg                                                                                >  3:655719/1‑149 (MQ=255)
                                                                      gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                                                               <  2:361184/149‑1 (MQ=255)
                                                                           ttCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                          >  1:95288/1‑149 (MQ=255)
                                                                                  cacaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttc                                                                   >  1:206154/1‑149 (MQ=255)
                                                                                   acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                  >  1:685361/1‑149 (MQ=255)
                                                                                   acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                  >  1:233894/1‑149 (MQ=255)
                                                                                    caACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTACTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                                                                 <  1:722720/149‑1 (MQ=255)
                                                                                     aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                                >  4:334404/1‑149 (MQ=255)
                                                                                      acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCCCGTTGTCGCCATAGCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCAATTCGCCTTCTTGc                                                               <  2:156069/149‑1 (MQ=255)
                                                                                           gggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                          >  3:302688/1‑149 (MQ=255)
                                                                                           gggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                          <  4:72857/149‑1 (MQ=255)
                                                                                              ggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGTTGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa                                                       <  1:945321/149‑1 (MQ=255)
                                                                                               gATGCCGTTAGCCACGGGAGATTTATCTCATAAATAATTAACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGCTGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAg                                                      <  4:910385/149‑1 (MQ=255)
                                                                                                     gTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATc                                                >  4:224976/1‑149 (MQ=255)
                                                                                                          ccACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGAt                                           <  3:1023942/149‑1 (MQ=255)
                                                                                                           cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATa                                          >  2:123153/1‑149 (MQ=255)
                                                                                                           cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATa                                          <  2:1097958/149‑1 (MQ=255)
                                                                                                                gagaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGt                                     >  2:584725/1‑149 (MQ=255)
                                                                                                                    tttATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGcc                                 >  3:721292/1‑149 (MQ=255)
                                                                                                                     ttATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGccc                                <  1:659817/149‑1 (MQ=255)
                                                                                                                        tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                             >  3:612509/1‑149 (MQ=255)
                                                                                                                          tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccg                           >  1:524034/1‑149 (MQ=255)
                                                                                                                            ataaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                         >  1:191739/1‑149 (MQ=255)
                                                                                                                            ataaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                         <  1:647574/149‑1 (MQ=255)
                                                                                                                             tatataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                        >  4:184077/4‑149 (MQ=255)
                                                                                                                             taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                        <  2:886255/149‑1 (MQ=255)
                                                                                                                             taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                        <  4:43876/149‑1 (MQ=255)
                                                                                                                              aaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCa                       >  2:733393/1‑149 (MQ=255)
                                                                                                                                 taatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGAt                    <  1:1064457/149‑1 (MQ=255)
                                                                                                                                   atTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTc                  <  2:785727/149‑1 (MQ=255)
                                                                                                                                   atTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTc                  <  1:792742/149‑1 (MQ=255)
                                                                                                                                    ttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGACCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCg                 >  2:961789/1‑149 (MQ=255)
                                                                                                                                    ttCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCg                 >  2:595031/1‑149 (MQ=255)
                                                                                                                                       aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc                                                     >  1:289564/1‑110 (MQ=255)
                                                                                                                                       aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt              >  4:575362/1‑149 (MQ=255)
                                                                                                                                       aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt              >  4:590707/1‑149 (MQ=255)
                                                                                                                                       aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt              >  4:914075/1‑149 (MQ=255)
                                                                                                                                       aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt              >  4:281903/1‑149 (MQ=255)
                                                                                                                                        cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTc             <  1:123153/149‑1 (MQ=255)
                                                                                                                                        cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTCAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCCAAGCTCTTCCAGATAGCCGTTGCCCGACCGACCAGATTCCCTTc             <  2:238973/149‑1 (MQ=255)
                                                                                                                                         gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc            <  2:524034/149‑1 (MQ=255)
                                                                                                                                         gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc            <  2:191739/149‑1 (MQ=255)
                                                                                                                                          ttGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCg           >  4:764022/1‑149 (MQ=255)
                                                                                                                                            gTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGc         >  3:893416/1‑149 (MQ=255)
                                                                                                                                               gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGcaca      >  4:1047776/1‑149 (MQ=255)
                                                                                                                                                 cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaaa    <  3:764022/149‑1 (MQ=255)
                                                                                                                                                  aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAg   >  1:1083091/1‑149 (MQ=255)
                                                                                                                                                   tAACGGCGACAATGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGc  >  2:520718/1‑149 (MQ=255)
                                                                                                                                                   tAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGc  <  2:138947/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGC  >  NC_000913/2704158‑2704453

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