Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,392,261 (A)8→9 intergenic (+112/‑99) orn → / → glyV oligoribonuclease/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,392,2531.A97.4% 1067.0 / 18.8 265intergenic (+104/‑107)orn/glyVoligoribonuclease/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base . (0/7);  new base A (133/125);  total (133/132)
Fisher's exact test for biased strand distribution p-value = 7.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

Only 100 of 306 total aligned reads displayed.

TTTTCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAA  >  NC_000913/4392153‑4392348
                                                                                                     |                                                                                               
ttttCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGc                                                                                                  >  1:5806142/1‑101 (MQ=255)
ttttCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGc                                                                                                  <  2:5805776/101‑1 (MQ=255)
 tttCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCa                                                                                                >  2:5806145/1‑101 (MQ=255)
   tCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGGGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGGGATAAAAACGCaaa                                                                                              <  1:5805786/101‑1 (MQ=255)
   tCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaa                                                                                              <  1:5805829/101‑1 (MQ=255)
    ccTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaa                                                                                             >  2:5806149/1‑101 (MQ=255)
     cTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaa                                                                                            >  1:5806150/1‑101 (MQ=255)
     cTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaa                                                                                            <  1:5805791/101‑1 (MQ=255)
      tGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaaa                                                                                           <  1:5805581/101‑1 (MQ=255)
      tGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaaa                                                                                           >  1:5806153/1‑101 (MQ=255)
       ggTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaaaa                                                                                          <  2:5805644/101‑1 (MQ=255)
           aCGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaaaaaact                                                                                      <  1:5805820/101‑4 (MQ=255)
             gTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCaaaaaaaaacttt                                                                                    <  2:5805767/101‑6 (MQ=255)
               aaGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAActtttt                                                                                  <  1:5805800/101‑1 (MQ=255)
               aaGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAActtttt                                                                                  <  2:5805808/101‑1 (MQ=255)
               aaGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAActtttt                                                                                  <  2:5805698/101‑1 (MQ=255)
                aGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACtttttt                                                                                  >  2:5806157/1‑101 (MQ=255)
                aGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACtttttt                                                                                  <  1:5805722/101‑1 (MQ=255)
                  cgcCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTgg                                                                                >  1:5806161/1‑101 (MQ=255)
                  cgcCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTgg                                                                                <  1:5805694/101‑1 (MQ=255)
                  cgcCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTgg                                                                                <  1:5805508/101‑1 (MQ=255)
                   gcCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTggg                                                                               >  2:5806162/1‑101 (MQ=255)
                      ggAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTgggggg                                                                            <  2:5805710/101‑1 (MQ=255)
                        aaTCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGtt                                                                          >  2:5806168/1‑101 (MQ=255)
                        aaTCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGtt                                                                          <  1:5805894/101‑1 (MQ=255)
                        aaTCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAAATTTTTTGGGGGGtt                                                                          >  1:5806167/1‑101 (MQ=255)
                            cACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCag                                                                      >  2:5806172/1‑101 (MQ=255)
                             aCTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCaga                                                                     <  1:5805743/101‑1 (MQ=255)
                              cTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCagag                                                                    >  2:5806175/1‑101 (MQ=255)
                                tGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAggg                                                                  <  1:5805648/101‑1 (MQ=255)
                                 gcggcggcgTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGa                                                                 <  1:5805937/101‑1 (MQ=255)
                                  cggcggcgTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGaa                                                                >  1:5806178/1‑101 (MQ=255)
                                    gcggcgTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAg                                                              <  1:5805762/101‑1 (MQ=255)
                                      ggcgTGATTTTGACGCTAAATTGAACGTCTTGCTGATAAGTTCAGCACTTGAGATAAAAACGCAAAAAAAAACtttttt                                                                                  <  2:5805703/79‑1 (MQ=255)
                                         gTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTtc                                                         <  2:5805870/101‑1 (MQ=255)
                                         gTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTtc                                                         >  1:5806182/1‑101 (MQ=255)
                                          tGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTtct                                                        <  2:5805757/101‑1 (MQ=255)
                                           gATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTtctc                                                       >  2:5806184/1‑101 (MQ=255)
                                            aTTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCg                                                      >  1:5806187/1‑101 (MQ=255)
                                              tttGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGta                                                    <  2:5805664/101‑1 (MQ=255)
                                               ttGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGtat                                                   <  1:5805733/101‑1 (MQ=255)
                                                tGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGtata                                                  >  2:5806192/1‑101 (MQ=255)
                                                  aCGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAAt                                                <  1:5805811/101‑1 (MQ=255)
                                                   cGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTtct                                                        <  2:5805904/92‑1 (MQ=255)
                                                    gCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATgc                                              >  2:5806197/1‑101 (MQ=255)
                                                     cTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATgcg                                             >  2:5806200/1‑101 (MQ=255)
                                                       aaaTTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGcc                                           >  1:5806204/1‑101 (MQ=255)
                                                       aaaTTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGcc                                           >  2:5806203/1‑101 (MQ=255)
                                                        aaTTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCt                                          <  2:5805966/101‑1 (MQ=255)
                                                           tGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTccc                                       >  2:5806208/1‑101 (MQ=255)
                                                           tGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTccc                                       <  2:5805741/101‑1 (MQ=255)
                                                           tGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTcc                                        <  2:5805480/100‑1 (MQ=255)
                                                              aCGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTa                                    >  1:5806211/1‑101 (MQ=255)
                                                               cGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTaa                                   >  2:5806215/1‑101 (MQ=255)
                                                               cGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTaa                                   >  1:5806216/1‑101 (MQ=255)
                                                                gTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAac                                  <  1:5805932/101‑1 (MQ=255)
                                                                 tCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacg                                 >  2:5806220/1‑101 (MQ=255)
                                                                 tCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacg                                 >  1:5806222/1‑101 (MQ=255)
                                                                  cTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacga                                <  2:5805997/101‑1 (MQ=255)
                                                                  cTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAac                                  <  2:5805838/99‑1 (MQ=255)
                                                                   ttGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCaga                                                                     <  2:5805816/63‑1 (MQ=255)
                                                                    tGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacgacg                              <  2:5805927/101‑1 (MQ=255)
                                                                    tGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacgacg                              <  2:5805914/101‑1 (MQ=255)
                                                                    tGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacgacg                              >  1:5806227/1‑101 (MQ=255)
                                                                    tGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAacgacg                              <  2:5805892/101‑1 (MQ=255)
                                                                      cTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCa                            >  2:5806229/1‑101 (MQ=255)
                                                                       tGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAg                           >  1:5806233/1‑101 (MQ=255)
                                                                         aTAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGaa                         <  2:5805583/101‑1 (MQ=255)
                                                                         aTAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGaa                         >  2:5806234/1‑101 (MQ=255)
                                                                         aTAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGaa                         <  1:5805906/101‑1 (MQ=255)
                                                                         aTAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGaa                         <  1:5805985/101‑1 (MQ=255)
                                                                          tAATTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGaaa                        >  1:5806238/1‑101 (MQ=255)
                                                                           aaTTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAAt                       >  2:5806241/1‑101 (MQ=255)
                                                                           aaTTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAAt                       <  1:5805878/101‑1 (MQ=255)
                                                                            aTTTCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATg                      >  2:5806244/1‑101 (MQ=255)
                                                                             tttCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGc                     >  2:5806248/1‑101 (MQ=255)
                                                                              ttCAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCg                    >  1:5806249/1‑101 (MQ=255)
                                                                                cAGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGaa                  <  1:5805860/101‑1 (MQ=255)
                                                                                 aGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGaaa                 <  2:5805785/101‑1 (MQ=255)
                                                                                 aGCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGaaa                 >  2:5806254/1‑101 (MQ=255)
                                                                                  gCACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGaaaa                >  1:5806256/1‑101 (MQ=255)
                                                                                   cACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAAt               >  2:5806259/1‑101 (MQ=255)
                                                                                   cACTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAAt               <  1:5805670/101‑1 (MQ=255)
                                                                                    aCTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAAtt              <  1:5805976/101‑1 (MQ=255)
                                                                                     cTTGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTa             >  1:5806262/1‑101 (MQ=255)
                                                                                      ttGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTAc            <  1:5806014/101‑1 (MQ=255)
                                                                                      ttGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTAc            >  1:5806267/1‑101 (MQ=255)
                                                                                      ttGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTAc            >  2:5806265/1‑101 (MQ=255)
                                                                                      ttGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTAc            <  2:5805682/101‑1 (MQ=255)
                                                                                       tGAGATAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGc                                                                        >  2:5806272/1‑40 (MQ=38)
                                                                                        gagaTAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGa          >  1:5806273/1‑101 (MQ=255)
                                                                                         agaTAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGaa         <  2:5805771/101‑1 (MQ=255)
                                                                                          gaTAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGaaa        <  1:5805845/101‑1 (MQ=255)
                                                                                          gaTAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGaaa        >  1:5806276/1‑101 (MQ=255)
                                                                                           aTAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAg       <  1:5805884/101‑1 (MQ=255)
                                                                                            tAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGc      <  2:5805730/101‑1 (MQ=255)
                                                                                            tAAAAACGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGc      >  2:5806281/1‑101 (MQ=255)
                                                                                             aaaaaCGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGAtt                                                           >  1:5806283/1‑47 (MQ=255)
                                                                                             aaaaaCGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCa     <  2:5806040/101‑1 (MQ=255)
                                                                                               aaaCGCAAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCaaa   >  1:5806288/1‑101 (MQ=255)
                                                                                                     |                                                                                               
TTTTCCTGGTCACGTAAGCGCCGGAATCCACTTGCGGCGGCGTGATTTTGACGCTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAA  >  NC_000913/4392153‑4392348

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