Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,308 C→T 6.6% intergenic (‑245/+27) 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,3080CT 6.6% 610.9 / 21.5 191intergenic (‑245/+27)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base C (94/84);  new base T (4/9);  total (98/93)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 191 total aligned reads displayed.

TAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATCC  >  NC_000913/2704074‑2704564
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                
    taaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                                                                                                                                                                                            <  4:240540/237‑1 (MQ=255)
                 agagacagacacGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcg                                                                                                                                                                                                          <  2:146128/266‑1 (MQ=255)
                                 cAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                                                                                                                                         <  4:65496/275‑1 (MQ=255)
                                                  gTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTTTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGa                                                                                                                                                                                                                                                 <  4:21553/202‑1 (MQ=255)
                                                  gTTCTGCTAACTGCTTCGCTGCAGTTTAAGTTTACCGGTATATGGGGACCAACGCTGGTTCGTGTGCCGTGGCAAACTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                                                                                                                                         <  2:18715/274‑1 (MQ=255)
                                                   ttCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGCCCGGTCGGCCATTAGCCTTCTTGCTTATCGa                                                                                                                                                                       >  2:626528/1‑275 (MQ=255)
                                                       gCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACt                                                                                                                                                                                                                <  2:369703/230‑1 (MQ=255)
                                                        cTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCAATTCGCCTTCTTGCTTATCGAAGCta                                                                                                                                                                  >  1:53818/1‑274 (MQ=255)
                                                         tAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACg                                                                                                                                                                                                                                                  <  2:344954/194‑1 (MQ=255)
                                                                            tcagtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACCACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGGGACTTAAGCACAGACCAGTCGGCCATTCGCCTTATTTGTTAACGAAGCCCAACCAAATTCCCGTTGccc                                                                                                                                               >  4:394535/1‑274 (MQ=255)
                                                                                        tataTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACgg                                                                                                                                                                                                                                                                <  6:420964/149‑1 (MQ=255)
                                                                                         atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCTAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                                                                                                                               <  6:22171/149‑1 (MQ=255)
                                                                                            tGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACCACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                                                                                                                                                                                            <  7:333066/149‑1 (MQ=255)
                                                                                              gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACCCAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACaa                                                                                                                                                                                                                                                          <  8:707477/149‑1 (MQ=255)
                                                                                              gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACaa                                                                                                                                                                                                                                                          <  8:381281/149‑1 (MQ=255)
                                                                                              gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACaa                                                                                                                                                                                                                                                          >  8:450019/1‑149 (MQ=255)
                                                                                               ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                                                                                                                                         <  6:944463/149‑1 (MQ=255)
                                                                                               ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                                                                                                                                         <  6:977877/149‑1 (MQ=255)
                                                                                                gACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                                                                                                                                                                                        >  7:584850/1‑149 (MQ=255)
                                                                                                  ccAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGc                                                                                                                        >  2:401449/1‑275 (MQ=255)
                                                                                                    aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                                                                                                                                                                                    >  7:72857/1‑149 (MQ=255)
                                                                                                    aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                                                                                                                                                                                    <  8:353099/149‑1 (MQ=255)
                                                                                                    aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaa                                                                                                                                                                                                                                                    >  8:766843/1‑148 (MQ=255)
                                                                                                        cTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGaa                                                                                                                                                                                                                                                <  6:174205/149‑1 (MQ=255)
                                                                                                        cTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGaa                                                                                                                                                                                                                                                <  5:560926/149‑1 (MQ=255)
                                                                                                         tGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAg                                                                                                                                                                                                                                               >  5:339146/1‑149 (MQ=255)
                                                                                                           tttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAt                                                                                                                                                                                                                                             <  7:209626/149‑1 (MQ=255)
                                                                                                           tttCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGAt                                                                                                                                                                                                                                             <  6:12820/149‑1 (MQ=255)
                                                                                                               gtgtGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCt                                                                                                                                                                                                                                         <  7:1068224/149‑1 (MQ=255)
                                                                                                                      gTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                                                                                                                                                                                  >  8:598216/1‑149 (MQ=255)
                                                                                                                       tGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATgg                                                                                                                                                                                                                                 >  7:650043/1‑149 (MQ=255)
                                                                                                                          cAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACAGAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATg                                                                                                                                                                                                                              >  5:866916/1‑149 (MQ=255)
                                                                                                                           aaCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTaa                                                                                                                                                                                                                                    <  2:397018/142‑1 (MQ=255)
                                                                                                                              cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                                                                                                                                          >  5:361184/1‑149 (MQ=255)
                                                                                                                                 gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                                                                                                                                       >  8:207212/1‑149 (MQ=255)
                                                                                                                                 gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                                                                                                                                       >  8:902627/1‑149 (MQ=255)
                                                                                                                                    gTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                                                                                                                                    >  5:455160/1‑149 (MQ=255)
                                                                                                                                     tGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCAACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCg                                                                                                                                                                                                                   >  6:658552/1‑149 (MQ=255)
                                                                                                                                       ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                                                                                                                                 >  8:167201/1‑149 (MQ=255)
                                                                                                                                       ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                                                                                                                                 >  8:763339/1‑149 (MQ=255)
                                                                                                                                       ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                                                                                                                                 >  7:283071/1‑149 (MQ=255)
                                                                                                                                            cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                                                                                                                                            >  8:277534/1‑149 (MQ=255)
                                                                                                                                            cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                                                                                                                                            >  7:950283/1‑149 (MQ=255)
                                                                                                                                             ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgc                                                                                                                                                                                                           >  7:43876/1‑149 (MQ=255)
                                                                                                                                               aaaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAACGCCCCAGTAACGGCCGTCGCCGCTGTTGTCGCGGtttgcc                                                                           >  4:552271/1‑271 (MQ=255)
                                                                                                                                                aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                                                                                                                                         <  6:777111/100‑1 (MQ=255)
                                                                                                                                                  cTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                                                                                                                                                                                      <  8:584850/149‑1 (MQ=255)
                                                                                                                                                   ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAc                                                                                                                                                                                                     <  5:77344/149‑1 (MQ=255)
                                                                                                                                                     cgtgggAACAGTCTTTACACAACTACTGAGTCAGCCGTTAGCAACGTGTAACTTATTTAATAATAAGTTCCAGTCGTCGGAATAATGTCGACAACATTTACGAATCTGTCTATTATACGATTCTGCCAATGCGACTGAGGCACAGACTATACGGCTATTCGCCTTCTTGCTTATCGAAGCACATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCCCAAAGCCACAGTAACGGCTGTCCGCGCTGTTGTCGTGGTTGCAGcccct                                                                      >  4:187090/4‑272 (MQ=255)
                                                                                                                                                     cGGGGGCATATTCTATACACAACGACGGGGGCTGTCGTTAGCCGCGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCANTGCGACTTAAGCGCAGACCAGTCGGCCCTTCTCCTTCTTGCTTATCGACGCTCATCAAGATAGCCGTTGNCCGACCGACCATATTCGCTTCCGGCCCCCACCCCCCCTAACGGCTNCNCCCCCTCTTGTCGGGGttattcccacg                                                                      >  3:57580/1‑265 (MQ=255)
                                                                                                                                                        gggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                                                                                                                                    >  5:970074/1‑129 (MQ=255)
                                                                                                                                                        gggAATATGCTAAACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATg                                                                  <  2:242362/275‑1 (MQ=255)
                                                                                                                                                         ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCg                                                                                                                                                                                               >  7:655719/1‑149 (MQ=255)
                                                                                                                                                         ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAg                                                                                                                                                                     <  4:697022/175‑1 (MQ=255)
                                                                                                                                                         ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCcg                                                                                               <  4:365023/245‑1 (MQ=255)
                                                                                                                                                          gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACtt                                                                                                                                                                                                               <  2:194054/132‑1 (MQ=255)
                                                                                                                                                          gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                                                                                                                                                                              <  6:361184/149‑1 (MQ=255)
                                                                                                                                                             tatTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTa                                                             <  2:574060/275‑1 (MQ=255)
                                                                                                                                                               ttCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                                                                                                                                         >  5:95288/1‑149 (MQ=255)
                                                                                                                                                                      cacaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttc                                                                                                                                                                                  >  5:206154/1‑149 (MQ=255)
                                                                                                                                                                       acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACg                                                                                                                                                                                                                                                  >  2:751859/1‑84 (MQ=255)
                                                                                                                                                                       acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                                                                                                                                 >  5:233894/1‑149 (MQ=255)
                                                                                                                                                                       acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                                                                                                                                 >  5:685361/1‑149 (MQ=255)
                                                                                                                                                                        caACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTACTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                                                                                                                                                                                <  5:722720/149‑1 (MQ=255)
                                                                                                                                                                         aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                                                                                                                                               >  8:334404/1‑149 (MQ=255)
                                                                                                                                                                         aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaa                                                                                                                    >  1:144542/1‑208 (MQ=255)
                                                                                                                                                                         aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCggag                                                  >  2:435918/1‑274 (MQ=255)
                                                                                                                                                                             acgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAGCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATc                                                                                                                                                               <  2:716455/161‑1 (MQ=255)
                                                                                                                                                                              cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCACGACCGACCAGATTCGCTTCCGGCACAAAGCCCAAGTAACGGCTGTCAGCGCTGGTGCCGGGGTCTGCGCCCCACATTGTGTATGGTCCCGGAGGAACag                                            >  4:73025/1‑274 (MQ=255)
                                                                                                                                                                              cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCAAAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGcc                                                                                                                <  2:679514/207‑1 (MQ=255)
                                                                                                                                                                                        cGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCTCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGTTTATCGAATCTCATCCAGATAGCCGTTGCCCGACCGACCTGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTTGCCCTTCATGAAGTATTTTCCTGGCGGAACAATCTAGGTTGc                                  >  1:140943/1‑275 (MQ=255)
                                                                                                                                                                                               cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATa                                                                                                                                                         <  6:1097958/149‑1 (MQ=255)
                                                                                                                                                                                                        tttATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGcc                                                                                                                                                >  7:721292/1‑149 (MQ=255)
                                                                                                                                                                                                         ttATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGccc                                                                                                                                               <  5:659817/149‑1 (MQ=255)
                                                                                                                                                                                                          tATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGgctg                <  2:223522/275‑1 (MQ=255)
                                                                                                                                                                                                            tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                                                                                                                                            >  7:612509/1‑149 (MQ=255)
                                                                                                                                                                                                              tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccg                                                                                                                                          >  5:524034/1‑149 (MQ=255)
                                                                                                                                                                                                                ataaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                                                                                                                                        >  5:191739/1‑149 (MQ=255)
                                                                                                                                                                                                                ataaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCcgaccgac                                                                                                                                        <  5:647574/149‑1 (MQ=255)
                                                                                                                                                                                                                 taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                                                                                                                                       <  6:886255/149‑1 (MQ=255)
                                                                                                                                                                                                                  aaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCa                                                                                                                                      >  6:733393/1‑149 (MQ=255)
                                                                                                                                                                                                                     taatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGAt                                                                                                                                   <  5:1064457/149‑1 (MQ=255)
                                                                                                                                                                                                                       atTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTc                                                                                                                                 <  5:792742/149‑1 (MQ=255)
                                                                                                                                                                                                                        ttacagTTGGCGCCAAATCGGCGAAAACGTGAACGATGATGGCTATTAATGGATGCCCCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGCTAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATcc  <  4:254303/270‑1 (MQ=255)
                                                                                                                                                                                                                        ttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGACCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCg                                                                                                                                >  6:961789/1‑149 (MQ=255)
                                                                                                                                                                                                                        ttCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCg                                                                                                                                >  6:595031/1‑149 (MQ=255)
                                                                                                                                                                                                                        ttCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATcc  <  4:435918/275‑1 (MQ=255)
                                                                                                                                                                                                                           aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc                                                                                                                                                                    >  5:289564/1‑110 (MQ=255)
                                                                                                                                                                                                                           aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt                                                                                                                             >  8:914075/1‑149 (MQ=255)
                                                                                                                                                                                                                            cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTc                                                                                                                            <  5:123153/149‑1 (MQ=255)
                                                                                                                                                                                                                            cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTCAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCCAAGCTCTTCCAGATAGCCGTTGCCCGACCGACCAGATTCCCTTc                                                                                                                            <  6:238973/149‑1 (MQ=255)
                                                                                                                                                                                                                             gTTGTCGCCATAACGGCGACAACGTGAACGAAGCTGACTATTAATGGATGCCACCAATGCGACTTAGGCGCAAACCGGTTGGCCATTCTTCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCTGATTCGCTTcc                                                                                                                           >  8:569782/1‑149 (MQ=255)
                                                                                                                                                                                                                             gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                                                                                                                                                                             >  1:174978/1‑83 (MQ=255)
                                                                                                                                                                                                                             gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc                                                                                                                           <  6:524034/149‑1 (MQ=255)
                                                                                                                                                                                                                                gTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGc                                                                                                                        >  7:893416/1‑149 (MQ=255)
                                                                                                                                                                                                                                   gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGcaca                                                                                                                     >  8:1047776/1‑149 (MQ=255)
                                                                                                                                                                                                                                    ccATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaa                                                                                                                    >  6:570357/1‑149 (MQ=255)
                                                                                                                                                                                                                                     cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaaa                                                                                                                   <  7:764022/149‑1 (MQ=255)
                                                                                                                                                                                                                                      aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAg                                                                                                                  >  7:194782/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGc                                                                                                                 <  6:138947/149‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                
TAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATCC  >  NC_000913/2704074‑2704564

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