Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,276 G→C 6.0% A160G (GCT→GGT)  obgE ← GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,331,2760GC 6.0% 223.1 / 4.7 67A160G (GCT→GGT) obgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base C (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AGCTCTTTTCGTTGTCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGT  >  NC_000913/3331130‑3331423
                                                                                                                                                  |                                                                                                                                                   
aGCTCTTTTCGTTGTCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTcagc                                                                                                                                                    >  8:246362/1‑148 (MQ=255)
     ttttCGTTGTCCATTCGTACCACACCCAGACTTGGCATCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTcag                                                                                                                                                     <  4:92990/142‑1 (MQ=255)
     ttttCGTTGTCCATTCGTACCACACCCAGACTTGGCATCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTcag                                                                                                                                                     >  3:92990/1‑142 (MQ=255)
     ttttCGTTGTCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTcagccagc                                                                                                                                                <  3:317182/147‑1 (MQ=255)
     ttttCGTTGTCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCATCcagcag                                                                                                                                              >  3:300296/1‑149 (MQ=255)
          gTTGTCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCa                                                                                                                                         <  7:180091/149‑1 (MQ=255)
              tCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTc                                                                                                                                     >  7:186409/1‑149 (MQ=255)
              tCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTc                                                                                                                                     >  2:338110/1‑149 (MQ=255)
                cATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGTACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCca                                                                                                                                   <  7:105000/149‑1 (MQ=255)
                  ttCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGAAAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAgcagca                                                                                                                                            <  8:379741/138‑1 (MQ=255)
                  ttCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGAAAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAgcagca                                                                                                                                            >  7:379741/1‑138 (MQ=255)
                     gTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCcagcagc                                                                                                                              <  7:199576/149‑1 (MQ=255)
                     gTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCcagcagc                                                                                                                              <  7:48519/149‑1 (MQ=255)
                        ccACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCcagcagcag                                                                                                                            <  8:374465/148‑1 (MQ=255)
                        ccACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCagcagcagc                                                                                                                           >  8:243157/1‑149 (MQ=255)
                        ccACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTACGGATTAGCCGCCGATACCGCACGAATAAAGGTCGATTAACACGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCAACAGGTCCagcagcagc                                                                                                                           >  4:301961/1‑149 (MQ=255)
                         cacacacaGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACTTCAGCCAGCAGCATCAGCTCCAGCAGCAGCt                                                                                                                          >  1:248733/1‑149 (MQ=255)
                            acCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCCCTTTCCGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCACCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTcgc                                                                                                                       >  1:181343/1‑149 (MQ=255)
                              ccAGACTTTGCACCAGAGTGTTAAACGGATAAACCGCCACTTTCGGTTTAGCCGCAGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACATCAGCCAGCAGCATCAGCTCAAGCAGCAGCTCgagc                                                                                                                     >  1:290032/1‑149 (MQ=255)
                                 gACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTa                                                                                                                  <  8:172960/149‑1 (MQ=255)
                                  aCTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTAt                                                                                                                 >  2:51721/1‑149 (MQ=255)
                                             gagTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTg                                                                                                      >  2:354982/1‑149 (MQ=255)
                                             gagTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTg                                                                                                      <  4:167681/149‑1 (MQ=255)
                                                tGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCg                                                                                                   >  4:269512/1‑149 (MQ=255)
                                                      aCGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGGCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGtt                                                                                             <  2:149728/149‑1 (MQ=255)
                                                        ggATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGtttt                                                                                           <  7:246362/149‑1 (MQ=255)
                                                        ggATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGtttt                                                                                           <  8:331130/149‑1 (MQ=255)
                                                        ggATAATCCGCCACTTTCGGGTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGGCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGtttt                                                                                           <  5:988/149‑1 (MQ=255)
                                                         gATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTc                                                                                          <  3:203225/149‑1 (MQ=255)
                                                         gATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTc                                                                                          <  8:186409/149‑1 (MQ=255)
                                                                     cattCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCACCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCg                                                                                                           >  1:201041/3‑120 (MQ=255)
                                                                     cTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCACCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCg                                                                                                           >  1:201042/1‑120 (MQ=255)
                                                                     cTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCACCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCg                                                                                                           <  2:201042/120‑1 (MQ=255)
                                                                     cTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCACCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCg                                                                                                           <  2:201041/120‑1 (MQ=255)
                                                                      tttCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTAc                                                                             >  8:121324/1‑149 (MQ=255)
                                                                           gTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                                         <  3:185084/148‑1 (MQ=255)
                                                                           gTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                                         >  4:185084/1‑148 (MQ=255)
                                                                           gTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                                         >  2:168555/1‑148 (MQ=255)
                                                                            tttAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCCCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                                                       >  3:3180/1‑149 (MQ=255)
                                                                                gccgccGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGa                                                                   <  6:111694/149‑1 (MQ=255)
                                                                                 ccgccgATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGGTTTCTGCCGCGGTGTACGGTTAACGGAc                                                                  >  6:141992/1‑149 (MQ=255)
                                                                                    ccgATACCGCACGAATAAAGGTCGATTTACCCGGGTTTGGCATCCCCAACATACCGCCGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGCCGTGCCGGTGGTTTGCTGCCGCGGTGCACGGTTATCGGACGAt                                                               <  6:265274/149‑1 (MQ=255)
                                                                                     cgATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGAtt                                                              <  1:354982/149‑1 (MQ=255)
                                                                                      gATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGAttt                                                             >  5:102684/1‑149 (MQ=255)
                                                                                          ccGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaaa                                                         <  6:133580/149‑1 (MQ=255)
                                                                                             cACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACgg                                                      <  7:121324/149‑1 (MQ=255)
                                                                                             cACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACgg                                                      >  5:246623/1‑149 (MQ=255)
                                                                                             cACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACgg                                                      <  5:195734/149‑1 (MQ=255)
                                                                                              aCGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACggg                                                     <  1:51721/149‑1 (MQ=255)
                                                                                               cGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTTGCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGt                                                    <  2:181343/149‑1 (MQ=255)
                                                                                                 aaTAAAGGTCGATTTACCCGCGTTTGGCATCCCGAACATACCGACGTCAGCCAGCAGCATCAGGTGCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCTTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTAt                                                  <  3:301961/149‑1 (MQ=255)
                                                                                                     aaGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGcc                                              <  1:338110/149‑1 (MQ=255)
                                                                                                         tCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGa                                          <  3:269512/149‑1 (MQ=255)
                                                                                                             tttACCCGCGTTTGACATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGt                                      >  6:107798/1‑149 (MQ=255)
                                                                                                               tACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTg                                     <  7:243157/148‑1 (MQ=255)
                                                                                                                aCCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAAc                                                                      >  6:50827/1‑114 (MQ=255)
                                                                                                                aCCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAAc                                                                      <  5:50827/114‑1 (MQ=255)
                                                                                                                          ggCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCt                          >  1:36743/1‑148 (MQ=255)
                                                                                                                           gCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGcc                                              >  6:61804/1‑127 (MQ=255)
                                                                                                                           gCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTa                        <  1:168555/149‑1 (MQ=255)
                                                                                                                            cATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACg                                                                     >  5:216753/1‑103 (MQ=255)
                                                                                                                            cATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACg                                                                     <  6:216753/103‑1 (MQ=255)
                                                                                                                                cccAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGAGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAAc                   <  2:248733/149‑1 (MQ=255)
                                                                                                                                 ccAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAAc                   >  3:163236/1‑148 (MQ=255)
                                                                                                                                   aaCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACcag                >  3:137114/1‑149 (MQ=255)
                                                                                                                                    aCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc                            <  7:135266/136‑1 (MQ=255)
                                                                                                                                    aCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc                            >  8:135266/1‑136 (MQ=255)
                                                                                                                                     cATACAGACGTCTGCCAGCAGCATCAGCTACAGCATCAGCTCGCTGTTATCGCCCGGCGTGCAGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTTTTGCCCAGACCGTGCAAGCCGCCCTTAGCAACcagca              <  2:290032/149‑1 (MQ=255)
                                                                                                                                        aCCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGAc           >  1:115989/1‑149 (MQ=255)
                                                                                                                                         ccGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGAc           <  5:107798/148‑1 (MQ=255)
                                                                                                                                          cGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCACGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTCCCAGCCGCCCTTAGCAACCAGCAGACGc         >  5:69235/1‑149 (MQ=255)
                                                                                                                                                 agccagcAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGCCGCTGACCgt  >  4:175985/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCTCTTTTCGTTGTCCATTCGTACCACACCCAGACTTGGCACCAGAGTGGTAAACGGATAATCCGCCACTTTCGGTTTAGCCGCCGATACCGCACGAATAAAGGTCGATTTACCCGCGTTTGGCATCCCCAACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGT  >  NC_000913/3331130‑3331423

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