Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,410,949 G→C 7.1% intergenic (+65/+64) dbpA → / ← ttcA ATP‑dependent RNA helicase DbpA/tRNA cytosine(32) 2‑sulfurtransferase TtcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,410,9490GC 7.1% 378.8 / 15.5 112intergenic (+65/+64)dbpA/ttcAATP‑dependent RNA helicase DbpA/tRNA cytosine(32) 2‑sulfurtransferase TtcA
Reads supporting (aligned to +/- strand):  ref base G (46/58);  new base C (6/2);  total (52/60)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 113 total aligned reads displayed.

GGCTGTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGG  >  NC_000913/1410803‑1411097
                                                                                                                                                  |                                                                                                                                                     
ggCTGTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCcg                                                                                                                                                    >  1:1526224/1‑150 (MQ=255)
  cTGTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCcgcg                                                                                                                                                  <  1:173444/150‑1 (MQ=255)
  cTGTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCAGATGCGACGCTGGCcgcg                                                                                                                                                  <  1:1055646/150‑1 (MQ=255)
  cTGTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGGTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCcgcg                                                                                                                                                  <  1:979169/150‑1 (MQ=255)
    gTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTc                                                                                                                                                <  1:520982/150‑1 (MQ=255)
    gTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTc                                                                                                                                                <  1:1564736/150‑1 (MQ=255)
    gTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTc                                                                                                                                                >  1:1455925/1‑150 (MQ=255)
    gTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTc                                                                                                                                                <  2:293332/150‑1 (MQ=255)
    gTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTc                                                                                                                                                <  1:1149479/150‑1 (MQ=255)
     ttGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCt                                                                                                                                               <  1:359129/150‑1 (MQ=255)
     ttGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCt                                                                                                                                               <  2:1438416/150‑1 (MQ=255)
      tGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATAAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCtt                                                                                                                                              <  2:1478398/150‑1 (MQ=255)
        cTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTAt                                                                                                                                            >  1:729875/1‑150 (MQ=255)
           aTAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCg                                                                                                                                         >  1:754362/1‑150 (MQ=255)
            tAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAACCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATAAGACGCGGCCAGCGTCGCATCCgg                                                                                                                                        <  2:991750/150‑1 (MQ=255)
               aGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCt                                                                                                                                    <  2:132996/150‑1 (MQ=255)
                gCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCtg                                                                                                                                   <  1:161186/150‑2 (MQ=255)
                gCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTc                                                                                                                                   >  1:2114/1‑150 (MQ=255)
                gCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTc                                                                                                                                   <  1:892874/150‑1 (MQ=255)
                 cATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTcc                                                                                                                                  >  1:896205/1‑150 (MQ=255)
                     gAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCTGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATaa                                                                                                                              >  2:361231/1‑150 (MQ=255)
                     gAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCCGGCCGCGTCTTATCCGGCCTCCATaa                                                                                                                              >  1:25514/1‑150 (MQ=255)
                      aaaCAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAag                                                                                                                             >  1:286430/1‑150 (MQ=255)
                      aaaCAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAag                                                                                                                             <  1:689117/150‑1 (MQ=255)
                       aaCAGTTACAGGGCTGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAaga                                                                                                                            <  1:604072/150‑1 (MQ=255)
                          aGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTa                                                                                                                         <  1:1346330/150‑1 (MQ=255)
                          aGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTa                                                                                                                         <  2:729875/150‑1 (MQ=255)
                              aCAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGccc                                                                                                                     >  2:1348814/1‑150 (MQ=255)
                                 gggCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGAt                                                                                                                  >  2:1260455/1‑150 (MQ=255)
                                  ggCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATa                                                                                                                 <  1:412120/150‑1 (MQ=255)
                                   gCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATAc                                                                                                                >  1:657256/1‑150 (MQ=255)
                                      ggAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCt                                                                                                             >  2:1579789/1‑150 (MQ=255)
                                          gATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTgcg                                                                                                         <  1:817058/150‑1 (MQ=255)
                                          gATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTgcg                                                                                                         <  2:1000990/150‑1 (MQ=255)
                                          gATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTgcg                                                                                                         <  2:286430/150‑1 (MQ=255)
                                           aTTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTgcgc                                                                                                        >  1:1066926/1‑150 (MQ=255)
                                            ttAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCa                                                                                                       >  2:873889/1‑150 (MQ=255)
                                            ttAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCa                                                                                                       >  2:310334/1‑150 (MQ=255)
                                              aaaGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGGAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATc                                                                                                     >  1:1641745/1‑150 (MQ=255)
                                              aaaGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATc                                                                                                     >  1:1647348/1‑150 (MQ=255)
                                                aGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCgg                                                                                                   <  2:133371/150‑1 (MQ=255)
                                                 ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCggg                                                                                                  >  1:1323667/1‑150 (MQ=255)
                                                 ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCggg                                                                                                  <  1:181764/150‑1 (MQ=255)
                                                 ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCAGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                  >  1:1498777/1‑150 (MQ=255)
                                                  gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCGGgc                                                                                                 <  1:349374/150‑1 (MQ=255)
                                                  gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCGGgc                                                                                                 <  1:930824/150‑1 (MQ=255)
                                                   aaaaaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgcg                                                                                                >  1:321365/1‑150 (MQ=255)
                                                    aaaaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgcgc                                                                                               >  1:816650/1‑150 (MQ=255)
                                                    aaaaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgcgc                                                                                               <  2:349741/150‑1 (MQ=255)
                                                       aCGTGCCGGGTGCGGTTTTTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTAGTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGTTACGCTTGCGCATCGGGCGCTag                                                                                            >  2:731497/1‑149 (MQ=255)
                                                                  gCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAttt                                                                                 >  1:1520394/1‑150 (MQ=255)
                                                                    ggTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCa                                                                               >  1:1265651/1‑150 (MQ=255)
                                                                       tattaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACtt                                                                            <  1:261257/150‑1 (MQ=255)
                                                                        attaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTc                                                                           <  1:170390/150‑1 (MQ=255)
                                                                         ttaAAATAATGAAATGTTGAATTGCCGGGTGCAGGAGTAAACATGTTATTCGGGATTGCCGGATGCGACGCTGGCAGCGTCTTATCCGGCCTGCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTAACTTCa                                                                          >  1:837419/1‑150 (MQ=255)
                                                                         ttaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCa                                                                          >  2:836189/1‑150 (MQ=255)
                                                                           aaaaTAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAAc                                                                        <  2:649778/150‑1 (MQ=255)
                                                                              ataatGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCtt                                                                                                                                              >  1:733725/1‑78 (MQ=255)
                                                                              ataatGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCtt                                                                                                                                              <  2:733725/78‑1 (MQ=255)
                                                                                aatGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACAt                                                                   >  2:1208866/1‑150 (MQ=255)
                                                                                           gAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAAc                                                        <  2:1607591/150‑1 (MQ=255)
                                                                                            aaTTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACt                                                       <  1:821194/150‑1 (MQ=255)
                                                                                            aaTTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACt                                                       <  1:411823/150‑1 (MQ=255)
                                                                                                gCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATc                                                   >  1:1160925/1‑150 (MQ=255)
                                                                                                gCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATc                                                   >  2:1369362/1‑150 (MQ=255)
                                                                                                 ccGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATcc                                                  <  2:944739/150‑1 (MQ=255)
                                                                                                 ccGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATcc                                                  <  2:657256/150‑1 (MQ=255)
                                                                                                       gCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGa                                            <  2:321365/150‑1 (MQ=255)
                                                                                                          agagTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTACATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGAttt                                         >  1:1329409/1‑150 (MQ=255)
                                                                                                             gTAAACATCTTATTCGGTCTTGCCTGATGCGAAGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCa                                      <  1:714772/150‑1 (MQ=255)
                                                                                                              tAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAt                                     >  1:1664222/1‑150 (MQ=255)
                                                                                                                aaCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtct                                   <  2:1251367/150‑1 (MQ=255)
                                                                                                                aaCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtct                                   <  2:611385/150‑1 (MQ=255)
                                                                                                                   aTCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttct                                <  2:1498777/150‑1 (MQ=255)
                                                                                                                      ttATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttcttct                             <  2:754362/150‑1 (MQ=255)
                                                                                                                         ttCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGc                          <  1:1348814/150‑1 (MQ=255)
                                                                                                                           cGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCTCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                        <  2:304230/150‑1 (MQ=255)
                                                                                                                           cGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                        <  2:1160925/150‑1 (MQ=255)
                                                                                                                             ggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGcc                      <  2:816650/150‑1 (MQ=255)
                                                                                                                             ggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGcc                      <  2:1066926/150‑1 (MQ=255)
                                                                                                                               aTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAg                    >  1:728121/1‑150 (MQ=255)
                                                                                                                                aggCCGGATAAGACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGc                   >  1:559196/3‑150 (MQ=255)
                                                                                                                                aggCCGGATAAGACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGc                   >  2:1213280/3‑150 (MQ=255)
                                                                                                                                    cGGATAAGACGCGGCCAGCGTCGCATCCGGCAATCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACg                >  2:188843/1‑150 (MQ=38)
                                                                                                                                    cGGATAAGACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGc               >  1:257977/1‑150 (MQ=255)
                                                                                                                                        taagACGCGGCCAGCGTCGCATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGt           >  1:1627778/4‑150 (MQ=255)
                                                                                                                                         aagACGCGGCCAGCGTCACATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGtt          >  2:335653/3‑150 (MQ=255)
                                                                                                                                          cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCCGCACGCCGGTTg         >  2:480454/1‑150 (MQ=255)
                                                                                                                                          cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg         >  1:514516/1‑150 (MQ=255)
                                                                                                                                          cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTAAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg         >  2:1357957/1‑150 (MQ=255)
                                                                                                                                           gACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGt        <  1:1260455/150‑1 (MQ=255)
                                                                                                                                              gCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                        <  1:73263/131‑1 (MQ=255)
                                                                                                                                              gCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                        >  2:73263/1‑131 (MQ=255)
                                                                                                                                              gCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGt     <  2:896205/150‑1 (MQ=255)
                                                                                                                                               gtGGCCGCGTCTTATCCGGCGTCCATAAGAGTAGGCATATACGTTTGCGCTTCGGGCGGTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGGCGGTTGTAGTg    <  2:1576264/149‑1 (MQ=255)
                                                                                                                                               cTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTg    <  2:1526224/150‑1 (MQ=255)
                                                                                                                                               cTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTg    <  1:568968/150‑1 (MQ=255)
                                                                                                                                                tGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCTTAACTCATTCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTgg   <  2:1334503/150‑1 (MQ=255)
                                                                                                                                                tGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTgg   >  2:1573992/1‑150 (MQ=255)
                                                                                                                                                 ggCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTggg  >  1:1602554/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCTGTTGCTCATAAAGCATGGAAACAGTTACAGGGCGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGG  >  NC_000913/1410803‑1411097

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