Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,398 T→C 5.8% intergenic (+22/+21) corA → / ← yigF magnesium/nickel/cobalt transporter/DUF2628 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,3980TC 5.8% 411.0 / 11.3 138intergenic (+22/+21)corA/yigFmagnesium/nickel/cobalt transporter/DUF2628 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (60/70);  new base C (4/4);  total (64/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

Only 100 of 139 total aligned reads displayed.

GCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAG  >  NC_000913/4002249‑4002546
                                                                                                                                                     |                                                                                                                                                     
gttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGt                                                                                                                                                       >  2:1686711/5‑150 (MQ=255)
  tATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt                                                                                                                                                     >  2:1203233/1‑150 (MQ=255)
  tATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt                                                                                                                                                     >  2:1389683/1‑150 (MQ=255)
  tATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt                                                                                                                                                     >  1:1588439/1‑150 (MQ=255)
   ctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGCACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg                                                                                                                                                    >  1:326921/2‑150 (MQ=255)
     ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc                                                                                                                                                  <  2:380460/150‑1 (MQ=255)
     ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc                                                                                                                                                  >  2:246682/1‑150 (MQ=255)
           aaCTTGGCGTGTCTGCCAGAACTGAAGGGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCCGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGcc                                                                                                                                           <  2:332231/150‑1 (MQ=255)
              tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCAcc                                                                                                                                        <  1:1484164/150‑1 (MQ=255)
              tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACATTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGAC‑TAAGCCAcc                                                                                                                                        <  1:113699/150‑1 (MQ=255)
               ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTAGCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGAC‑TAAGCCAccg                                                                                                                                       <  2:274351/150‑2 (MQ=255)
                tGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTg                                                                                                                                      <  2:803634/150‑1 (MQ=255)
                  aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGtt                                                                                                                                    >  1:782436/1‑150 (MQ=255)
                  aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑GAAGCCACCTGtt                                                                                                                                    >  1:1713595/1‑150 (MQ=255)
                    tttATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTAt                                                                                                                                  <  2:1358538/150‑1 (MQ=255)
                      tATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTc                                                                                                                                >  1:1628206/1‑150 (MQ=255)
                      tATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCATAGTCTGGC‑TAAGCCACCTGTTATTc                                                                                                                                >  2:805077/1‑150 (MQ=255)
                         gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCaaa                                                                                                                             >  1:1352466/1‑150 (MQ=255)
                         gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCaaa                                                                                                                             <  1:333926/150‑1 (MQ=255)
                         gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCaaa                                                                                                                             >  1:1112278/1‑150 (MQ=255)
                         gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCaaa                                                                                                                             >  2:1419871/1‑150 (MQ=255)
                          ccAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTTCTTTAAGCGGAAGAACTGGTTGTAAAAAACGCGAGCGGTGGCTTTGTCTGGC‑TAAGCCCCCTGTTATTCAAAg                                                                                                                            >  2:1101686/1‑150 (MQ=255)
                            aGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGc                                                                                                                          >  1:1016030/1‑150 (MQ=255)
                             gAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCt                                                                                                                         >  1:414806/1‑150 (MQ=255)
                              aaCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGGGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTc                                                                                                                        <  1:1171672/150‑1 (MQ=255)
                               aCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTcc                                                                                                                       >  2:653050/1‑150 (MQ=255)
                               aCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTAACATTAAGATCCTCGCGGGCCAGGCACCCTATCTGGACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTTTTATTCAAAAGCTcc                                                                                                                       >  1:1670619/1‑150 (MQ=255)
                                cttaaGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCTCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCa                                                                                                                      <  1:1042269/147‑1 (MQ=255)
                                 tGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAg                                                                                                                     >  1:400302/1‑150 (MQ=255)
                                    aGTGGAGCTTCTGCTACCCTGGGGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTGGTAAAAAACGAGAGCGGGGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTa                                                                                                                  <  1:670320/150‑1 (MQ=255)
                                    aGTGGAGCTTCGGCTCCCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAACCGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAggga                                                                                                                  >  1:82457/1‑148 (MQ=255)
                                    aGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTa                                                                                                                  <  2:576459/150‑1 (MQ=255)
                                       ggAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTAttt                                                                                                               <  2:657862/150‑1 (MQ=255)
                                       ggAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTAttt                                                                                                               <  2:1282009/150‑1 (MQ=255)
                                         aGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTTAGCGGAAGAACTGGTTTTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTaa                                                                                                             >  2:488653/1‑150 (MQ=255)
                                          gcctCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAc                                                                                                            <  1:660343/147‑1 (MQ=255)
                                            ttCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAccc                                                                                                          <  2:596345/150‑1 (MQ=255)
                                            ggcGGCTACCCTGGCGCGATTATCTTTATCATCCTCGCGGGCCTGGCACCGTATCTGTACGTTAAGCGGAAGAACTGGTTGTAAAACACGAGAGCGGTGGCTTAGTCGGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAccc                                                                                                          <  2:272781/148‑1 (MQ=255)
                                                gCTACCCTGGCGCGATTATCTTTTTGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑GAAGCCCCCTGTTATTCAAAGGCTCCAGGTATTTAACCCgttt                                                                                                      >  2:1423575/1‑150 (MQ=255)
                                                gCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCtttt                                                                                                      <  2:504679/150‑1 (MQ=255)
                                                       tGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCtt                                                                                               <  1:1043782/150‑1 (MQ=255)
                                                         gcgcTATTATCTTTATGATCCTCGCGTGCCTGGCACCGTATCTGTACTTTAAGAGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGTCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTtc                                                                                             <  1:227131/150‑1 (MQ=255)
                                                            cgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCa                                                                                          <  1:1592243/150‑1 (MQ=255)
                                                            cgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCa                                                                                          <  1:969881/150‑1 (MQ=255)
                                                              aTTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATa                                                                                        <  2:1365840/150‑1 (MQ=255)
                                                                   cTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAGGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAcc                                                                                   >  2:1531731/1‑150 (MQ=255)
                                                                   cTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCAGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAcc                                                                                   >  2:1517677/1‑150 (MQ=255)
                                                                      tATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAtt                                                                                <  2:469645/150‑1 (MQ=255)
                                                                       aTGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAttt                                                                               <  1:1368936/150‑1 (MQ=255)
                                                                          aTCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGATCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCGGTTATTCAAAGGCTCCAGGTAATTAACCCTTTTACCTCTTTCTCATAGAACCATTTGtt                                                                            <  1:698700/150‑1 (MQ=255)
                                                                            ccTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCg                                                                          <  2:939258/150‑1 (MQ=255)
                                                                             cTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCgt                                                                         >  1:958156/1‑150 (MQ=255)
                                                                               cgcgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCgtg                                                                        <  2:642433/150‑1 (MQ=255)
                                                                                gcgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGtt                                                                      <  1:28970/150‑1 (MQ=255)
                                                                                 cgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTa                                                                     <  1:940700/150‑1 (MQ=255)
                                                                                   ggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAAc                                                                   >  1:535352/1‑150 (MQ=255)
                                                                                      cTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGc                                                                >  1:286519/1‑150 (MQ=255)
                                                                                       tGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCa                                                               >  1:1639974/1‑150 (MQ=255)
                                                                                         gCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAAc                                                             <  1:1267339/150‑1 (MQ=255)
                                                                                         gCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAAc                                                             <  2:889175/150‑1 (MQ=255)
                                                                                             cGTATCAGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTATTTCTAATAGAACCATTTGTTCGTGTTAACAGCAAAATAg                                                         >  1:877810/1‑150 (MQ=255)
                                                                                               tATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgc                                                       <  1:949837/150‑1 (MQ=255)
                                                                                                 tCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctg                                                     >  1:147282/1‑150 (MQ=255)
                                                                                                 tCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctg                                                     >  1:1434770/1‑150 (MQ=255)
                                                                                                  cTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgc                                                    >  2:182695/1‑150 (MQ=255)
                                                                                                   tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                   <  2:1373356/150‑1 (MQ=255)
                                                                                                   tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACAATTTGTTCGTGATAACATCAACATAGgctgct                                                   >  2:1311434/1‑150 (MQ=255)
                                                                                                      aCTTTAAGCGGAAGATCTGGTTGTAAAAAACGAGATCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACg                                                <  2:324902/150‑1 (MQ=255)
                                                                                                       cTTTAAGCGGAAGAACTGGTTGTAATAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTCACAGCAACATAGGCTGCTACgg                                               <  2:677572/150‑1 (MQ=255)
                                                                                                       cTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACgg                                               >  1:1470809/1‑150 (MQ=255)
                                                                                                       cTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACgg                                               >  2:506435/1‑150 (MQ=255)
                                                                                                       cTTTAAGCGGAAGAACAGGTTGTAAAAAACGAGAGCTGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACgg                                               <  2:1124824/150‑1 (MQ=255)
                                                                                                          tAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAACATTTTACATCTTTCACATAGACCCATTAGTTCGTCTTAACAGCAACATAGGCTGCTACGGCaa                                            >  2:1124040/1‑150 (MQ=255)
                                                                                                           aaGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAAt                                           >  2:229489/1‑150 (MQ=255)
                                                                                                             gCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACTGCAATTc                                         >  1:13621/1‑150 (MQ=255)
                                                                                                             gCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCGTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTAAACCCTTTAACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTc                                         <  2:264694/150‑1 (MQ=255)
                                                                                                               ggAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTccc                                       <  2:1588439/150‑1 (MQ=255)
                                                                                                               ggAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTccc                                       >  2:1515369/1‑150 (MQ=255)
                                                                                                                gaagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCATGTTAACAGCAACATAGGCTGCTACGGCAATTCCCa                                      <  1:1201323/150‑1 (MQ=255)
                                                                                                                gaagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACACTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCa                                      <  2:810634/150‑1 (MQ=255)
                                                                                                                 aagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCGTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAg                                     <  2:110722/150‑1 (MQ=255)
                                                                                                                   gaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGaa                                   <  2:206508/150‑1 (MQ=255)
                                                                                                                         gTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCTTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTaa                             >  1:684560/1‑150 (MQ=255)
                                                                                                                            gTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCACTTCCCAGCATGTGTACGc                          >  1:311662/1‑150 (MQ=255)
                                                                                                                             tAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGcc                         >  1:517538/1‑150 (MQ=255)
                                                                                                                             tAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGcc                         >  1:511034/1‑150 (MQ=255)
                                                                                                                              aaaaaaCGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCa                        >  1:1415556/1‑150 (MQ=255)
                                                                                                                                aaaaCGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAt                      <  2:414806/150‑1 (MQ=255)
                                                                                                                                      agagCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGc                <  2:1484300/150‑1 (MQ=255)
                                                                                                                                       gagCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACATCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCa               <  2:1628206/150‑1 (MQ=255)
                                                                                                                                        acaggTGGCTTAGCCAGAC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCAc              >  2:1709056/4‑150 (MQ=255)
                                                                                                                                          aggTGGCTTAGCCAGAC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACca            >  1:598527/2‑150 (MQ=255)
                                                                                                                                            gTGGCTTAGCCAGAC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTTACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaac          >  1:1708339/1‑150 (MQ=255)
                                                                                                                                             tGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaaca         >  1:638539/1‑150 (MQ=255)
                                                                                                                                              ggCTTAGCCAGAC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaacaa        >  2:1665674/1‑150 (MQ=255)
                                                                                                                                                cTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAAt      >  2:1242059/1‑150 (MQ=255)
                                                                                                                                                  tAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATag    <  2:638539/150‑1 (MQ=255)
                                                                                                                                                  tAGTCTGGC‑TAAGCCACCGGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATag    <  2:326921/150‑1 (MQ=255)
                                                                                                                                                   aGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATaga   >  2:533028/1‑150 (MQ=255)
                                                                                                                                                    gTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGGACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATagag  <  2:400302/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGC‑TAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAG  >  NC_000913/4002249‑4002546

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