Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,275,154 T→A 5.0% intergenic (‑113/‑317) yjcB ← / → yjcC putative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,275,1540TA 5.0% 433.3 / 6.3 144intergenic (‑113/‑317)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (61/75);  new base A (3/5);  total (64/80)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

GCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAAT  >  NC_000913/4275006‑4275302
                                                                                                                                                    |                                                                                                                                                    
gCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtg                                                                                                                                                     >  1:1175907/1‑150 (MQ=255)
gCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtg                                                                                                                                                     >  1:807523/1‑150 (MQ=255)
 cGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgt                                                                                                                                                    >  2:1499433/1‑150 (MQ=255)
   aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                                                                                                                                  <  1:1947692/150‑1 (MQ=255)
    agaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttg                                                                                                                                                 <  2:1121525/150‑1 (MQ=255)
      aagaaCCACGCCAGTGGTGAGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACTGTAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTgtttgttt                                                                                                                                              <  1:1391105/150‑1 (MQ=255)
       tgaaACACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGACAAGAGCCGAGGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTgtttgttt                                                                                                                                              <  2:514821/149‑1 (MQ=255)
       agaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTgtttgttt                                                                                                                                              >  1:732064/1‑150 (MQ=255)
        gaaCCACGCCAGTGGTGAGGGTACCCATCAAACAATTCCTCAGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGTAGTCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTc                                                                                                                                             <  1:267021/150‑1 (MQ=255)
         caCCCCGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCTCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCa                                                                                                                                            <  1:183226/149‑1 (MQ=255)
          aCCACGCCAGTGGTGAGGGTAGCCATCAACCAATTCCTCTGAAAAGAGCCGAAGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGCTCACAAAAATGAATCAGGCAGTTTGTTTCAt                                                                                                                                           >  1:1022932/1‑150 (MQ=255)
          aCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAt                                                                                                                                           >  1:152209/1‑150 (MQ=255)
           ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAtt                                                                                                                                          >  1:198845/1‑150 (MQ=255)
              cGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAAt                                                                                                                                       <  2:1567923/150‑1 (MQ=255)
              cGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAAt                                                                                                                                       >  1:53701/1‑150 (MQ=255)
              cGCCAGTGGTGAGGGTAGCCATAAAACAATTACTCTGAAAAGAGCCGATGACCTGGACAGCGCGGCGTAGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAAAGAAACACTCTGTATGTTTCATTAAt                                                                                                                                       >  1:744532/1‑150 (MQ=255)
               gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                      <  2:942219/150‑5 (MQ=255)
                  aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtg                                                                                                                                   <  1:446921/150‑1 (MQ=255)
                  aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGACGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttgtgt                                                                                                                                   <  2:121680/150‑8 (MQ=255)
                   gtggtgAGGGGAGCCATCACACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgt                                                                                                                                  <  1:152012/150‑1 (MQ=255)
                    tggtgAGGGTAGCCATCAAACAATTACTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgtg                                                                                                                                 >  1:71887/1‑150 (MQ=255)
                     ggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTGTGTTTCATTAATTTTGTGa                                                                                                                                >  1:337248/1‑150 (MQ=255)
                       tgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAAc                                                                                                                              >  1:1430207/1‑150 (MQ=255)
                         aGGGTAGCCATCACACAATTCCTCTGAAAAGAGCCGATGCCCTGGACATCTCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTTCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACta                                                                                                                            <  2:642996/150‑1 (MQ=255)
                          gggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACtat                                                                                                                           >  1:106998/1‑150 (MQ=255)
                           ggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACtata                                                                                                                          <  1:420297/150‑1 (MQ=255)
                           ggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACtata                                                                                                                          <  1:818045/150‑1 (MQ=255)
                           ggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACtata                                                                                                                          >  2:1035292/1‑150 (MQ=255)
                             tAGACATCAAACAATTCCTCGGAAAAGAGCCGTTGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATc                                                                                                                        <  2:386333/150‑1 (MQ=255)
                              aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATca                                                                                                                       <  2:1976916/150‑1 (MQ=255)
                                 cATCAAACAATTCCTCTGAAAAGAGCCGATGCACTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAACCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACaa                                                                                                                    >  2:621710/1‑150 (MQ=255)
                                   tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAatt                                                                                                                  >  1:725872/1‑150 (MQ=255)
                                             ccTCTGACACGTGGCGATGCCCTGGTCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACTCTCTGGCGTTCGCGGGGCTTGCTGAAGTCAAAAATGAAACACTCAGTTTGTTTCATTTATTTTGTGAACTTTATCACTATAGTttgtttgt                                                                                                        <  2:1086614/150‑1 (MQ=255)
                                              ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACGCAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                                       >  1:122088/1‑150 (MQ=255)
                                              ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                                       >  1:1197878/1‑150 (MQ=255)
                                              ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTa                                                                                                      >  2:1337897/1‑150 (MQ=255)
                                               tctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTa                                                                                                      <  1:1261680/150‑1 (MQ=255)
                                                gtGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAg                                                                                                     <  1:933932/149‑1 (MQ=255)
                                                 tGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGc                                                                                                    >  1:1800471/1‑150 (MQ=255)
                                                 tGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGc                                                                                                    <  2:631638/150‑1 (MQ=255)
                                                  gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGcc                                                                                                   >  1:625327/1‑150 (MQ=255)
                                                  gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGcc                                                                                                   >  2:1515718/1‑150 (MQ=255)
                                                    aaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCat                                                                                                 >  1:1130924/1‑150 (MQ=255)
                                                    aaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTGCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCat                                                                                                 >  1:996446/1‑150 (MQ=255)
                                                         gCCGAAGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTTGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGGTTGTTAGCCATATTAg                                                                                            <  2:1670275/150‑1 (MQ=255)
                                                           cGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGc                                                                                          >  2:816580/1‑150 (MQ=255)
                                                            gATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGcc                                                                                         >  1:1239757/1‑150 (MQ=255)
                                                               gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTg                                                                                      <  1:703150/150‑1 (MQ=255)
                                                                       cAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAt                                                                              >  2:1343574/1‑150 (MQ=255)
                                                                        aGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAtt                                                                             <  2:1935035/150‑1 (MQ=255)
                                                                         gcgcGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTg                                                                            >  2:33542/1‑150 (MQ=255)
                                                                          cgcgGCGGTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCCCAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGc                                                                           <  1:274412/150‑1 (MQ=255)
                                                                                gttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACa                                                                     <  1:648047/150‑1 (MQ=255)
                                                                                gttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACa                                                                     >  1:215600/1‑150 (MQ=255)
                                                                                   gttCACGGGAGGTGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAt                                                                  <  2:1637448/150‑1 (MQ=255)
                                                                                     tCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTa                                                                <  2:775553/150‑1 (MQ=255)
                                                                                     tCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTa                                                                <  2:1633789/150‑1 (MQ=255)
                                                                                       ccGGGAGGCGTGGTTCACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACATTTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTAGAGATTAtg                                                              <  2:927670/149‑1 (MQ=255)
                                                                                       aCGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtg                                                              >  1:681329/1‑150 (MQ=255)
                                                                                        cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                             <  2:566037/150‑1 (MQ=255)
                                                                                         gggAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgtg                                                            <  1:897546/150‑1 (MQ=255)
                                                                                         gggAGGCGTGATACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATCAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgtg                                                            <  1:1880882/150‑1 (MQ=255)
                                                                                             ggCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGttt                                                        >  1:1033378/1‑150 (MQ=255)
                                                                                              gCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTACATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtttt                                                       <  2:993026/150‑1 (MQ=255)
                                                                                                gTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                                     <  2:328004/150‑1 (MQ=255)
                                                                                                   gTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGt                                                  >  2:1170043/1‑150 (MQ=255)
                                                                                                    tACACTCTGTCTTTCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTg                                                 <  2:1057575/150‑1 (MQ=255)
                                                                                                     acacTCTTGCTATCGCGTGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATCTCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTgg                                                <  2:378177/150‑1 (MQ=255)
                                                                                                     acacTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACATATTATGTGGTTTTTCAGTgg                                                >  2:1097595/1‑150 (MQ=255)
                                                                                                       acTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCa                                              >  2:355841/1‑150 (MQ=255)
                                                                                                         tctGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                                            <  1:1270407/150‑1 (MQ=255)
                                                                                                            ggCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAg                                         >  1:2000246/1‑150 (MQ=255)
                                                                                                             gCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAgg                                        <  1:415693/150‑1 (MQ=255)
                                                                                                               tATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                      >  1:1389589/1‑150 (MQ=255)
                                                                                                                  cgcgGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCata                                   >  2:443380/1‑150 (MQ=255)
                                                                                                                      gggCTTGCAGAACACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCtt                               <  2:1430207/150‑1 (MQ=255)
                                                                                                                        gCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                             <  1:1643686/150‑1 (MQ=255)
                                                                                                                          ttGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccg                           <  2:1996862/150‑1 (MQ=255)
                                                                                                                            gCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgcc                         <  2:144858/150‑1 (MQ=255)
                                                                                                                            gCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgcc                         >  2:163854/1‑150 (MQ=255)
                                                                                                                             cAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgccg                        <  2:1800471/150‑1 (MQ=255)
                                                                                                                                 caaaattaaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctc                    >  2:736864/8‑150 (MQ=255)
                                                                                                                                 acacAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctc                    >  1:847117/1‑150 (MQ=255)
                                                                                                                                  cacaAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctct                   <  2:1389589/150‑1 (MQ=255)
                                                                                                                                  cacaAAAATGAAACACTCTCTATAAATCATTAATTTTGAGAACAATATCACAATTGATTGTTAGTTAGCCATATTAGGACGTGACTATTATTGCTGAACAGATAATTTCGTTTTTCAGTGGCATTAAGGGCATATCTACCCGCACCctct                   >  2:1684698/1‑150 (MQ=255)
                                                                                                                                   acaAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTg                  <  2:1239757/150‑1 (MQ=255)
                                                                                                                                   acaAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTg                  <  2:1867447/150‑1 (MQ=255)
                                                                                                                                    caAAAATGAAACACACTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGc                 >  1:574445/1‑150 (MQ=255)
                                                                                                                                       aaaTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAtt              <  1:385102/150‑1 (MQ=255)
                                                                                                                                        aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTc             >  1:832234/1‑150 (MQ=255)
                                                                                                                                         aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc            <  1:1499433/150‑1 (MQ=255)
                                                                                                                                         aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGGTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc            <  1:848276/150‑1 (MQ=255)
                                                                                                                                           gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCAGCCGCCTCTGCATTCCtt          >  1:743472/1‑149 (MQ=255)
                                                                                                                                           gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTATGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGGTTTTCAGTGGCAGTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTg          <  2:1630601/150‑1 (MQ=255)
                                                                                                                                              acacTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAg       <  2:1175907/150‑1 (MQ=255)
                                                                                                                                              acacTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAg       <  1:1188421/150‑1 (MQ=255)
                                                                                                                                                 ctctGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGaa    <  2:725872/150‑1 (MQ=255)
                                                                                                                                                  tatgtttgttTCATTAATTTTGTGAACTTTATCACAATTGATTGTTTGTTAGCAATATTAGGCCGTGACTTTTATTGATGTACAGCTTATGTGGTTTTTCATAGGCAATAAGGGCATATCTTCCCGCCTCCTCTGCATTCCAGTAGTaaa   >  1:1493980/3‑150 (MQ=255)
                                                                                                                                                   ctGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaat  <  2:68537/150‑1 (MQ=255)
                                                                                                                                                   ctGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaat  <  2:841000/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAAT  >  NC_000913/4275006‑4275302

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