Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,4960AG41.4% 81.8 / 27.0 58pseudogene (1083/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (4/30);  new base G (24/0);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.60e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTC  >  NC_000913/4573350‑4573643
                                                                                                                                                   |                                                                                                                                                      
tgtgAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAggg                                                                                                                                                     <  2:571223/150‑1 (MQ=255)
tgtgAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAggg                                                                                                                                                     <  1:48509/150‑1 (MQ=255)
tgtgAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAggg                                                                                                                                                     <  2:364033/150‑1 (MQ=255)
    aaGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGggtg                                                                                                                                                  >  1:681711/1‑149 (MQ=255)
         ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGAttgtt                                                                                                                                            >  2:149798/1‑150 (MQ=255)
         ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGGGAtt                                                                                                                                               >  2:219441/1‑147 (MQ=255)
          gTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTgg                                                                                                                                          <  2:674541/151‑1 (MQ=255)
            gAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGGTTGTTtggg                                                                                                                                        >  1:810365/1‑150 (MQ=255)
            gAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGATTGTTgggg                                                                                                                                        >  2:200092/1‑150 (MQ=255)
             aGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAg                                                                                                                                       <  2:245471/151‑1 (MQ=255)
              ggAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGAGGGGGGGTTGTTtgggg                                                                                                                                       >  1:984891/1‑148 (MQ=255)
                  gTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTa                                                                                                                                 <  2:112526/151‑1 (MQ=255)
                  gTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGGTTTTTGGGGGCAtt                                                                                                                                  >  1:297161/1‑150 (MQ=255)
                      ctGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTaa                                                                                                                                <  1:265300/148‑1 (MQ=255)
                        ggCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTtg                                                                                                                           <  1:347035/151‑1 (MQ=255)
                        ggCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGATGATTGTTGGGGGGATTAAACtt                                                                                                                            >  2:658455/1‑150 (MQ=255)
                         gCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTtga                                                                                                                          <  2:754625/151‑1 (MQ=255)
                         gCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGATTGGGGGGGGATTGTTGGGGGGCTTTAACtt                                                                                                                            >  2:649938/1‑149 (MQ=255)
                              gggCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACa                                                                                                                   <  1:258583/151‑1 (MQ=255)
                                gCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGt                                                                                                                 <  2:146003/151‑1 (MQ=255)
                                  cAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGATTTTTGGGGGCATTAAAAtt                                                                                                                            >  2:824904/1‑140 (MQ=255)
                                        aTGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTg                                                                                                          <  2:386800/150‑1 (MQ=255)
                                          gCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGcc                                                                                                        <  1:45386/150‑1 (MQ=255)
                                              ctctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCAtt                                                                                                     <  2:882901/149‑1 (MQ=255)
                                               tctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGGTTGTTGGGGGGATTAAAATTGATGACAGTGAAAATGCCAtt                                                                                                     >  2:548445/1‑148 (MQ=255)
                                                ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTgg                                                                                                   <  1:308059/149‑1 (MQ=255)
                                                  ggAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGt                                                                                                  <  1:935301/148‑1 (MQ=255)
                                                   gAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAgg                                                                                              <  1:541630/151‑1 (MQ=255)
                                                    aaaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGGGGTTGTTGGGGGCATTAAACTTTATGACAGTGAAAATGCCCTTGGGGAgg                                                                                              >  1:347611/1‑150 (MQ=255)
                                                    aaaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGt                                                                                             <  2:6685/151‑1 (MQ=255)
                                                        cagacagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCg                                                                                         <  1:545197/151‑1 (MQ=255)
                                                           acagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTa                                                                                      <  1:31744/151‑1 (MQ=255)
                                                             agGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGCGCTGGGGGGGTGGTTGTTGGGGGCATTAAAATTGATGGCAGGGAAATTGCCATTGGTGGGGGGCCggtta                                                                                     >  1:316810/1‑146 (MQ=255)
                                                              ggCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGGTTTTTGGGGGGATTAAAATTTATGGCAGTGAAAATGCCATTTGGGAGGGGGCGTTAAcc                                                                                   >  1:283578/1‑151 (MQ=255)
                                                                   cAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGATTGTTTGGGGGATTAAAATTTATGACAGGGGAAATGCCCTTGGTGGGGGGGCGTTAAACtt                                                                                 >  2:980985/1‑148 (MQ=255)
                                                                     gTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGGTAGTGGCGCGTGGTGGGGGGGGGGTTGTTGGGGGCATTTAAACTTATGTAAAGTGAAGTTCCAATGGGGGGGGGGCCGTTAACCtcatt                                                                              >  2:465521/1‑146 (MQ=255)
                                                                       tCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGATTGTTGGGGGCATTTAAATTTGTGACAGGGAAATTGCCATTTGTGGGGGGCCGGTAAACTTCAt                                                                              >  1:839718/1‑147 (MQ=255)
                                                                        cGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGtt                                                                          <  2:406891/150‑1 (MQ=255)
                                                                          aTTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGAt                                                                       <  1:401100/151‑1 (MQ=255)
                                                                           ttACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTCGCTGTGGGGGGATGATTGTTGGGAGCCTTTAAATTTTTGAAAGTGAAAATTCCATTGGGGGGGTGCGGCTAAACTTCATCGtttgt                                                                       >  1:81383/1‑147 (MQ=255)
                                                                            tACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATg                                                                      <  1:239581/150‑1 (MQ=255)
                                                                            tACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGGGGCTCGTGGCGGGGGGGGGGTTTTTTGGGGGATTTAAATTTGTTGAAAGGGAAATTTCCTTTGGGGGGGGGCCGTTAACCTTCTCGtt                                                                          >  1:81391/1‑146 (MQ=255)
                                                                             aCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                                     <  2:231981/150‑1 (MQ=255)
                                                                             aCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGGTTTGTGGGGGGATTTAAATTTGTGACAGGGGAATTGTCCTTTGGGGGGGGGCGTTAACCTTTATCGtt                                                                          >  2:163938/1‑145 (MQ=255)
                                                                                     gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCGGGGGGGGAGATTGTTGGGGGCATTAAACTTTATGACAGTGAAATTGCCATTGGGGGGGGGCCGGTAACCTTCATCGTTGATGGAGATCaa                                                              >  2:795508/1‑149 (MQ=255)
                                                                                           tGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGAATGTTGTTGGGGGCGTTTAAATTTTTGACGAGTGAAATGCCATTGGTTGGGGGCCCGTAACCTTTAATCTTGGTGGAGATCAATGGGTg                                                        >  2:608426/1‑149 (MQ=255)
                                                                                            gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTCAACCTTCATCGTTGATGCAGATCAATGGCTGTTa                                                     <  1:520694/151‑1 (MQ=255)
                                                                                                             aTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGATTGTTGGGGGCATTAAAATTGGTGACAGTGGAAATGCCATTGGGGAGGGGGCGTTAAACTTTATCGTTGGTGCAGATCAATGGCTGATACAGGGACAGGcttctt                                     >  2:658161/1‑150 (MQ=255)
                                                                                                               ccATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGGTTTTTGGGAGCATTAAAATTGATGACAGGGAAATTGCCATTGGGGAGGGGCCGTTAAACTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAAGcttctt                                     >  1:841466/1‑148 (MQ=255)
                                                                                                               ccATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCa                                  <  2:763080/151‑1 (MQ=255)
                                                                                                                cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGGTTGTTGGGGGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGGGGGGCCGTTTACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGAAAAGcttctt                                     >  1:500888/1‑147 (MQ=255)
                                                                                                                     gCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGAAGTGTGc                            <  2:479870/151‑1 (MQ=255)
                                                                                                                         gTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAag                        <  2:984788/151‑1 (MQ=255)
                                                                                                                            tGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagca                      <  2:156751/150‑1 (MQ=255)
                                                                                                                                 ttAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTTAAATTGGTGACAGTGAAATTGCCATTGGTGGGGGGCCGTTAACCTTCATCGTTGATGCAGATCAATGGGTGTTACAGGGACAGGCTTCTTGCAGTGTGCGGAGCAGCg                   >  2:846261/1‑148 (MQ=255)
                                                                                                                                          cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtg       <  1:133835/151‑1 (MQ=255)
                                                                                                                                           gTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCTGCGATGTTCTGATtgt      >  1:621492/1‑151 (MQ=255)
                                                                                                                                                tGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTc  >  2:334751/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTC  >  NC_000913/4573350‑4573643

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