Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,3921.G41.5% 42.0 / 137.2 65intergenic (+563/+173)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (12/26);  new base G (4/23);  total (16/49)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-02

CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTT‑GTCGAACCCCGG  >  NC_000913/1097255‑1097486
                                                                                                                                            |                                                                                               
cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTg                                                                                      <  2:1058228/149‑1 (MQ=11)
cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTg                                                                                      <  2:25563/149‑1 (MQ=11)
     cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑Tg‑ggggggggtgcgcacccc                                                                               >  1:810979/1‑140 (MQ=11)
     cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCGCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACt                                                                                   <  1:1357518/147‑1 (MQ=17)
      gTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑Gg‑gggggggggttggaccccc                                                                              >  2:981232/1‑139 (MQ=1)
       ttCACGGTCTGTCCAACTGGCTGTCGCCATTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACtc                                                                                  <  1:562378/146‑1 (MQ=18)
        tCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACtcgctt                                                                              <  2:594949/150‑1 (MQ=9)
         cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTGAGtcccct                                                                              >  2:48806/1‑144 (MQ=9)
          aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGG‑GGGGGGGGTTTGCtcgcttcgc                                                                           >  2:1718930/1‑151 (MQ=11)
           cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTGCGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTtcgctc                                                                         <  2:336188/151‑1 (MQ=18)
           cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGT‑GGT‑GA‑GGGGGGGATTGACTCGCTtcgctc                                                                         >  1:2656/1‑151 (MQ=18)
           cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGAAATATACG‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTtcgctc                                                                         <  1:1222970/150‑1 (MQ=16)
            ggTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTGACTCCCCccgctc                                                                         >  2:1533363/1‑150 (MQ=255)
             gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTtcgctcg                                                                        <  2:1008746/151‑1 (MQ=9)
             gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTGAcacccccccctccc                                                                       >  4:560079/1‑138 (MQ=255)
             gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTtcgctcgc                                                                       <  3:18647/151‑1 (MQ=33)
                 gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGcc                                                                      <  4:432720/149‑1 (MQ=9)
                 gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTg                                                                   <  2:866453/151‑1 (MQ=33)
                  tCCAACTGGCTGTCGCCTGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTtcgctcgc                                                                       <  1:1040810/147‑1 (MQ=9)
                   ccAACTGGCTGTGGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGc                                                                  <  2:115679/151‑1 (MQ=9)
                     aaCTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg                                                                <  3:528189/151‑1 (MQ=9)
                     aaCTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCggg                                                               <  2:1130624/151‑1 (MQ=14)
                       cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑GA‑GGGGGGGGTTGACTCGCTTCGCCCGCCCCGCCGGcc                                                             >  2:979253/1‑150 (MQ=255)
                         ggCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                                             <  2:1409130/149‑1 (MQ=33)
                          gCTGTCGCCAGTTGTCGAACCCCGGACGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTACGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                                             <  2:1221839/148‑1 (MQ=18)
                             gTCGCCAGTTGTCGGGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                                                         <  1:1282882/150‑1 (MQ=9)
                             gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCg                                                        <  2:1017470/150‑1 (MQ=33)
                             gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGc                                                       <  3:272871/151‑1 (MQ=33)
                               cGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                                                         <  1:1699087/148‑1 (MQ=9)
                               cGCCAGTTGCCGAACCCCGGTCGCGGCTTCTCGCCCCCCCGGTGTGTGCAATATTCGAAAAAAAAGCCCTTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc                                                     <  1:473732/151‑1 (MQ=11)
                                 ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc                                                   <  2:669107/151‑1 (MQ=33)
                                   aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                  <  1:1547391/151‑1 (MQ=9)
                                   aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                  <  1:1043576/151‑1 (MQ=9)
                                   aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                  <  2:152099/151‑1 (MQ=9)
                                         cGAGCCCCGTTCGGGGCTTCTCATCCCCCCGGTGTTTGCAATATACGAAAAAAAAGCCCGTTCTTTCGTTCGAGCTCTTCTTTAAATATGGC‑GGT‑TAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc                                            <  3:114423/151‑1 (MQ=11)
                                         cGAGCCCCGCTGGGGGGTTCCGGTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGGCCGTTCTTTCGGACGAGCTCTTCTTTAAAAATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt                                              <  1:701590/149‑1 (MQ=11)
                                         cGAACCCCGGTCGGGGCTTCTCATCCCCCCGTTGTGTGAAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTCCTTTAAATATGGC‑GGT‑AA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCa                                           <  1:467533/151‑1 (MQ=255)
                                         cGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc                                            <  2:350876/150‑1 (MQ=33)
                                           agccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGG‑GGGGGGGGGTGACTCGCTTCCCCCCCCCCGCGGGGAAcccccccccctcgtt                                             >  3:485482/3‑136 (MQ=11)
                                            gccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACg                                         <  2:1533406/150‑1 (MQ=14)
                                             ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCAc                                          <  4:300778/148‑1 (MQ=31)
                                             ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGt                                       <  2:1096375/151‑1 (MQ=31)
                                              cctggTCGGGCCTTTTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACgg                                        <  1:1480167/147‑1 (MQ=255)
                                               ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTc                                      <  3:61414/151‑1 (MQ=14)
                                                cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCt                                     <  3:236899/151‑1 (MQ=14)
                                                   tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGATGAATCGCTTCGCTCGCCCCGCCGGCAGCCCGCGCACTggg                                               >  1:1328917/1‑135 (MQ=11)
                                                    cGGGGCTTCTTATCCCCCCGGTGTCTCTAACATTCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACgg                                        <  4:483088/143‑1 (MQ=11)
                                                    cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGGTTGAATCGCTTCGCTCGCCCCCCGGGCAGCCCCCCCACCCCGTTCACCGGCTGt                                   >  2:1357666/1‑149 (MQ=11)
                                                    cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTcaa                                <  2:518649/151‑3 (MQ=14)
                                                       ggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGGCCGTACTTTCGGACGAGCTCTTCTTTAAAGATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAc                              <  1:7929/151‑1 (MQ=255)
                                                          ttCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt                                              <  4:373007/132‑1 (MQ=17)
                                                           tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGAACCGCCTCGCCCCCCCCCCGCGCGCCCCCCCCCCCGCGTTCACGGTGTGTccca                               >  1:1608104/1‑144 (MQ=11)
                                                           tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTgg                           <  2:1736425/149‑1 (MQ=31)
                                                            ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTGAATCCCCTCGCTCGCCCTGCCGGCAGCCCCCTCACCGCGTTCACCGTCTGTTCCACTgg                           >  2:1448566/1‑148 (MQ=11)
                                                             tcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑GA‑GGGGGGGAGTGAATATCTCCGCTCGCCCCGCCGGCCGCCCCCCCACTGCGGTCACGGCCTGttc                                 >  1:1314335/1‑140 (MQ=11)
                                                              cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCt                         <  2:621095/149‑1 (MQ=14)
                                                                tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTGACTCGCCTTGCTCCCCCTCCGGGCAACCCCCTCTCTGCGGTCACTGCCGGTCCCACCGGCAGTCGc                    >  2:1314191/1‑151 (MQ=11)
                                                                  cccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCa                  <  2:911066/151‑1 (MQ=31)
                                                                   ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTTCTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCa                  <  1:391832/151‑1 (MQ=2)
                                                                   ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTGCTTACGACTTCTTCTTTAAATATGGC‑GGT‑GC‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAg                 <  1:654572/151‑1 (MQ=2)
                                                                   ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTCCGTACGGGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAg                 <  1:25602/151‑1 (MQ=11)
                                                                            gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTT‑GTCGaa        <  2:1215460/151‑1 (MQ=11)
                                                                                 atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGG‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTT‑GTCGAACCCCg   <  1:1645472/151‑1 (MQ=9)
                                                                                 atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTT‑GTCGAACCCCg   <  2:845484/151‑1 (MQ=14)
                                                                                  tataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACGCCCCGCGCTCTCCCCTCGGGCAGCCCCCCCCCCGCGTTCCCGGTCCTTTCAACCGGCTGCGGCCCGGT‑GGCCAAccccc   >  1:1418811/1‑148 (MQ=2)
                                                                                   ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGGTTGACTCGCCTCGCTCGCCCCGCGGGCAGCCCCCCCACTGCGTTTCCGGTCTGTCCAACTGGGTGTCGCCCGTT‑GTCGGACCCCgg  >  2:158656/1‑150 (MQ=11)
                                                                                    taCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGGTTGACTCGCCTTCCCCCCCCCGCGGGCAGCCCCCCCACCGCGTTCAAGGGCTGGTCAACCCGCCGTTGCCCGTTGGTCGAACCCCg   >  2:467607/1‑149 (MQ=2)
                                                                                                                                            |                                                                                               
CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTT‑GTCGAACCCCGG  >  NC_000913/1097255‑1097486

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