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

CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG  >  NC_000913/1097255‑1097486
                                                                                                                                            |                                                                                               
cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTg                                                                                      <  2:1058228/149‑1 (MQ=11)
cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTg                                                                                      <  2:25563/149‑1 (MQ=11)
     cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTgggggggggtgcgcacccc                                                                               >  1:810979/1‑140 (MQ=11)
     cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCGCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACt                                                                                   <  1:1357518/147‑1 (MQ=17)
      gTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGggggggggggttggaccccc                                                                              >  2:981232/1‑139 (MQ=1)
       ttCACGGTCTGTCCAACTGGCTGTCGCCATTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtc                                                                                  <  1:562378/146‑1 (MQ=18)
        tCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcgctt                                                                              <  2:594949/150‑1 (MQ=9)
         cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAGtcccct                                                                              >  2:48806/1‑144 (MQ=9)
          aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTTTGCtcgcttcgc                                                                           >  2:1718930/1‑151 (MQ=11)
           cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTGCGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctc                                                                         <  2:336188/151‑1 (MQ=18)
           cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGTGGTGAGGGGGGGATTGACTCGCTtcgctc                                                                         >  1:2656/1‑151 (MQ=18)
           cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGAAATATACGAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctc                                                                         <  1:1222970/150‑1 (MQ=16)
            ggTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCCCCccgctc                                                                         >  2:1533363/1‑150 (MQ=255)
             gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTtcgctcg                                                                        <  2:1008746/151‑1 (MQ=9)
             gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAcacccccccctccc                                                                       >  4:560079/1‑138 (MQ=255)
             gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctcgc                                                                       <  3:18647/151‑1 (MQ=33)
                 gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGcc                                                                      <  4:432720/149‑1 (MQ=9)
                 gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTg                                                                   <  2:866453/151‑1 (MQ=33)
                  tCCAACTGGCTGTCGCCTGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTtcgctcgc                                                                       <  1:1040810/147‑1 (MQ=9)
                   ccAACTGGCTGTGGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGc                                                                  <  2:115679/151‑1 (MQ=9)
                     aaCTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg                                                                <  3:528189/151‑1 (MQ=9)
                     aaCTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCggg                                                               <  2:1130624/151‑1 (MQ=14)
                       cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTTGACTCGCTTCGCCCGCCCCGCCGGcc                                                             >  2:979253/1‑150 (MQ=255)
                         ggCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                                             <  2:1409130/149‑1 (MQ=33)
                          gCTGTCGCCAGTTGTCGAACCCCGGACGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTACGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                                             <  2:1221839/148‑1 (MQ=18)
                             gTCGCCAGTTGTCGGGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                                                         <  1:1282882/150‑1 (MQ=9)
                             gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCg                                                        <  2:1017470/150‑1 (MQ=33)
                             gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGc                                                       <  3:272871/151‑1 (MQ=33)
                               cGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                                                         <  1:1699087/148‑1 (MQ=9)
                               cGCCAGTTGCCGAACCCCGGTCGCGGCTTCTCGCCCCCCCGGTGTGTGCAATATTCGAAAAAAAAGCCCTTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc                                                     <  1:473732/151‑1 (MQ=11)
                                 ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc                                                   <  2:669107/151‑1 (MQ=33)
                                   aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                  <  1:1547391/151‑1 (MQ=9)
                                   aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                  <  1:1043576/151‑1 (MQ=9)
                                   aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                  <  2:152099/151‑1 (MQ=9)
                                         cGAGCCCCGTTCGGGGCTTCTCATCCCCCCGGTGTTTGCAATATACGAAAAAAAAGCCCGTTCTTTCGTTCGAGCTCTTCTTTAAATATGGCGGTTAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc                                            <  3:114423/151‑1 (MQ=11)
                                         cGAGCCCCGCTGGGGGGTTCCGGTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGGCCGTTCTTTCGGACGAGCTCTTCTTTAAAAATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt                                              <  1:701590/149‑1 (MQ=11)
                                         cGAACCCCGGTCGGGGCTTCTCATCCCCCCGTTGTGTGAAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTCCTTTAAATATGGCGGTAAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCa                                           <  1:467533/151‑1 (MQ=255)
                                         cGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc                                            <  2:350876/150‑1 (MQ=33)
                                           agccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGGTGACTCGCTTCCCCCCCCCCGCGGGGAAcccccccccctcgtt                                             >  3:485482/3‑136 (MQ=11)
                                            gccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACg                                         <  2:1533406/150‑1 (MQ=14)
                                             ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCAc                                          <  4:300778/148‑1 (MQ=31)
                                             ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGt                                       <  2:1096375/151‑1 (MQ=31)
                                              cctggTCGGGCCTTTTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACgg                                        <  1:1480167/147‑1 (MQ=255)
                                               ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTc                                      <  3:61414/151‑1 (MQ=14)
                                                cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCt                                     <  3:236899/151‑1 (MQ=14)
                                                   tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATGAATCGCTTCGCTCGCCCCGCCGGCAGCCCGCGCACTggg                                               >  1:1328917/1‑135 (MQ=11)
                                                    cGGGGCTTCTTATCCCCCCGGTGTCTCTAACATTCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACgg                                        <  4:483088/143‑1 (MQ=11)
                                                    cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAATCGCTTCGCTCGCCCCCCGGGCAGCCCCCCCACCCCGTTCACCGGCTGt                                   >  2:1357666/1‑149 (MQ=11)
                                                    cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTcaa                                <  2:518649/151‑3 (MQ=14)
                                                       ggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGGCCGTACTTTCGGACGAGCTCTTCTTTAAAGATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAc                              <  1:7929/151‑1 (MQ=255)
                                                          ttCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt                                              <  4:373007/132‑1 (MQ=17)
                                                           tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGAACCGCCTCGCCCCCCCCCCGCGCGCCCCCCCCCCCGCGTTCACGGTGTGTccca                               >  1:1608104/1‑144 (MQ=11)
                                                           tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTgg                           <  2:1736425/149‑1 (MQ=31)
                                                            ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAATCCCCTCGCTCGCCCTGCCGGCAGCCCCCTCACCGCGTTCACCGTCTGTTCCACTgg                           >  2:1448566/1‑148 (MQ=11)
                                                             tcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGAGGGGGGGAGTGAATATCTCCGCTCGCCCCGCCGGCCGCCCCCCCACTGCGGTCACGGCCTGttc                                 >  1:1314335/1‑140 (MQ=11)
                                                              cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCt                         <  2:621095/149‑1 (MQ=14)
                                                                tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTGACTCGCCTTGCTCCCCCTCCGGGCAACCCCCTCTCTGCGGTCACTGCCGGTCCCACCGGCAGTCGc                    >  2:1314191/1‑151 (MQ=11)
                                                                  cccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCa                  <  2:911066/151‑1 (MQ=31)
                                                                   ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTTCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCa                  <  1:391832/151‑1 (MQ=2)
                                                                   ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTGCTTACGACTTCTTCTTTAAATATGGCGGTGCGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAg                 <  1:654572/151‑1 (MQ=2)
                                                                   ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTCCGTACGGGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAg                 <  1:25602/151‑1 (MQ=11)
                                                                            gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa        <  2:1215460/151‑1 (MQ=11)
                                                                                 atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCg   <  1:1645472/151‑1 (MQ=9)
                                                                                 atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCg   <  2:845484/151‑1 (MQ=14)
                                                                                  tataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACGCCCCGCGCTCTCCCCTCGGGCAGCCCCCCCCCCGCGTTCCCGGTCCTTTCAACCGGCTGCGGCCCGGTGGCCAAcccc  >  1:1418811/1‑148 (MQ=2)
                                                                                   ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGACTCGCCTCGCTCGCCCCGCGGGCAGCCCCCCCACTGCGTTTCCGGTCTGTCCAACTGGGTGTCGCCCGTTGTCGGACCCCgg  >  2:158656/1‑150 (MQ=11)
                                                                                    taCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGACTCGCCTTCCCCCCCCCGCGGGCAGCCCCCCCACCGCGTTCAAGGGCTGGTCAACCCGCCGTTGCCCGTTGGTCGAACCCCg   >  2:467607/1‑149 (MQ=2)
                                                                                                                                            |                                                                                               
CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG  >  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: