Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,3640CG44.0% 152.4 / 14.3 75A71P (GCC→CCC) yjiRputative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (5/37);  new base G (33/0);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 2.79e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGT  >  NC_000913/4571218‑4571510
                                                                                                                                                  |                                                                                                                                                   
gggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCttt                                                                                                                                                  <  2:643444/150‑1 (MQ=255)
   cgcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGt                                                                                                                                              <  2:741009/151‑1 (MQ=255)
    gcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGtt                                                                                                                                             >  2:84336/1‑151 (MQ=255)
         gACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCg                                                                                                                                        <  2:277600/151‑1 (MQ=255)
           cATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGTTGTGCggc                                                                                                                                      >  2:1296639/1‑149 (MQ=255)
               cGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTTTGGGGCaaa                                                                                                                                  >  1:190424/1‑150 (MQ=255)
               cGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGGTGTGCggccaa                                                                                                                                   >  2:217414/1‑145 (MQ=255)
                ggTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGCGGGGGTTTCGTTGTGCGGCaaaa                                                                                                                                 >  1:1246143/1‑150 (MQ=255)
                 gtgtGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATa                                                                                                                               <  2:382509/151‑1 (MQ=255)
                  tgtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGGTTTTCGGTGTGGGGCaaaaaa                                                                                                                               >  1:1220297/1‑148 (MQ=255)
                   gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCCTTGTGGGGCAaaaaaa                                                                                                                              >  2:1314627/1‑150 (MQ=255)
                    tgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTTCGGTGTGCGGCAaaaaaaac                                                                                                                            >  2:11409/1‑149 (MQ=255)
                      cTTTTGCTTAACGGAACAATGGAACTCTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTGTTCGACAAAAAAAaccc                                                                                                                          >  2:165771/1‑150 (MQ=255)
                       ttttGCTTAACGGAACAATGGAACTTTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTGCTGGCCGGCGGGGTTTCAATTGTGTGGGAAAATaaaaccg                                                                                                                         >  2:39251/1‑146 (MQ=255)
                        tttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCg                                                                                                                          <  2:800633/149‑1 (MQ=255)
                            cTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGTTTTGCGACAAAAAAAACCCAAcgcg                                                                                                                    >  2:20820/1‑151 (MQ=255)
                              tAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCCCTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAcgcg                                                                                                                    <  1:463701/149‑1 (MQ=255)
                                aCGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTCGGTTTGCGGCAAAAAAAACCCCAAGCGGcgg                                                                                                                >  2:1312115/1‑149 (MQ=255)
                                 cGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGc                                                                                                               <  1:1315441/151‑1 (MQ=255)
                                  ggAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGCGGGGTTTTCGTTTTGGGGCAAAAAAAACCCAAAGCGGCcgcg                                                                                                              >  2:685258/1‑151 (MQ=255)
                                  ggAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCg                                                                                                              <  2:353391/151‑1 (MQ=255)
                                  ggAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTACACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCg                                                                                                              <  2:1158479/151‑1 (MQ=255)
                                      cAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGt                                                                                                          <  1:836246/151‑1 (MQ=255)
                                        aTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCTGGGGTTTTCGGTTTTGGGCAAAAAAAACCCAAAACGGCTGCggggggg                                                                                                        >  1:1031186/1‑146 (MQ=255)
                                         tGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCTGGGGTTTTCGGTGTGGGGCAAAAAAAACCCAACGCGGCCGCGGAGggat                                                                                                       >  2:1315226/1‑151 (MQ=255)
                                                gTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTCGTTTTGGGACAAAAAAAACCCAACGCGGCCGCGGGGTGGTGAGGTTc                                                                                                >  2:614175/1‑151 (MQ=255)
                                                gTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCtt                                                                                                 <  2:1149242/150‑1 (MQ=255)
                                                  cGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCa                                                                                               <  1:409958/150‑1 (MQ=255)
                                                   gCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATc                                                                                             <  1:1300259/151‑1 (MQ=255)
                                                   gCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATc                                                                                             <  1:446677/151‑1 (MQ=255)
                                                           ccTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGCTTTCGTTGTGGGGCAAAAAAAACCGCACGCGGCCGCGGGGGGGTGAGCTTCATCGTCTcccg                                                                                     >  2:426504/1‑149 (MQ=255)
                                                            cTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTGTGCGGCAAAAAAAACCCAACGCGGCTGCGGGGGGGTGAGCTTTATCCTCTCCCGc                                                                                    >  2:35065/1‑151 (MQ=255)
                                                            cTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTTCGTTTTGGGGCCAAAAAAACCCAACACGGCCGCGGGGGGGTTAGCtttctcctccccccc                                                                                    >  1:200968/1‑138 (MQ=255)
                                                              ccAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGc                                                                                    <  2:1596611/149‑1 (MQ=255)
                                                               cAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTCGTTTTGCGGCAAAAAAACCCCAACCCGGCCGGGGGGGGGTGAGGTTCATCGTCTCCcgcggc                                                                                 >  2:1430703/1‑149 (MQ=255)
                                                               cAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTTCGTTTTTGGGCAAAAAAAACCCCACACGGCCGCGGGGGGGTGAGGTTTATTGTCTcccccggc                                                                                 >  1:1418590/1‑145 (MQ=255)
                                                                     atCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTCGGTTTGCGGCAAAAAAAACCCCAACCGGCTGCGGGGGGGTGAGGTTCATTCTCTCCCGCGGCTGATAc                                                                           >  2:1500050/1‑151 (MQ=255)
                                                                     atCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGGGTTTCGGTGTGCGGCAAAATAACCCCAAAGCGGGCGCGGGGGGATGGGGTTCATCGTCTCCcgcgcccggta                                                                            >  2:1282825/1‑143 (MQ=255)
                                                                       ccAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGc                                                                         <  1:11782/151‑1 (MQ=255)
                                                                         aGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCct                                                                       <  1:834300/151‑1 (MQ=255)
                                                                             ccTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctgctg                                                                   <  2:1012515/151‑1 (MQ=255)
                                                                                 aTCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAATGGCGGCTTGGGTTTTCGGTGTGGGGCAAAAAAAACCCAAACCGGCCGCGGGGGGGTGAGGTTTTTCCTCTCCCGCGGCTGGTAAGCccgccggccg                                                               >  1:568684/1‑145 (MQ=255)
                                                                                   cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGCGGGGTTTTCGTTTGGCGGCAAAAAAAACCCAAAGCCGCCGCGGGGGGGTGAGCGTTATCGTCTCCCCCGTCTGATaaccccccctccccgg                                                             >  1:625897/1‑136 (MQ=255)
                                                                                   cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGCGGGGTTTTCCGTGTGCGGCAAAAAAAACCGCAAGCGGGCGCGGGGGGATGAGCTTCATCGTCTCCCGCGGCTGATAAGCCCGCTGCACggg                                                             >  2:732078/1‑150 (MQ=255)
                                                                                   cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACgg                                                              <  1:1074896/150‑1 (MQ=255)
                                                                                   cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACgg                                                              <  1:102116/150‑1 (MQ=255)
                                                                                   cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGGTTGGCGGCAAAAAAAACCCAAAGCGGCCGCGGGGGGGTGAGGTTCATTGTCTCCCGCGTCTTATACGCCcgctgccaggg                                                             >  2:48224/1‑146 (MQ=255)
                                                                                    ccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTGTGCGGCAAAAAAAACCCAAAGCGGGTGCGGGGGGGTGAGCTTTCTCGTCTCCCACGGCTGATAACCCCGCTGCACgggg                                                            >  2:459876/1‑149 (MQ=255)
                                                                                    ccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGGTTTGCGGCAAAATAACCCGAACGCGGCTGCGGGGGGGTGAGCTTCATCGTCTCCAGCGGCTGATACGCCCGCTGCCCgggg                                                            >  2:998724/1‑149 (MQ=255)
                                                                                          ggTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATa                                                      <  2:307362/151‑1 (MQ=255)
                                                                                            tAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCCGGGGGTTTCCGTGTTGGGCCAAATAAACCCAACGCGGGGGCGGGGGGGTGAGCTTTCTCGCCTCCAGCGGCTGATAAGGCCGCTGCCCGGTGGTGGTACt                                                    >  1:223353/1‑151 (MQ=255)
                                                                                             aaTTTCCACCGGGCGCTGCCCCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTCACGTTGTGCGACAAAATACCCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTg                                                   <  1:1262978/151‑1 (MQ=255)
                                                                                                 tCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACg                                                <  2:1091941/150‑1 (MQ=255)
                                                                                                    aCCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTGTGGGGCAAAAAAAACCCAACGCGGCCGCGGGGGGGTGAGGTTTATCGTCTCCCGCGGCTGATACGCCTGCTGCACGGGGCTGATACTGACGCCgg                                            >  2:435985/1‑150 (MQ=255)
                                                                                                        ggCGCTGCACCGGACGCGTCATCGGCGGTACCTGCGGCTGGGGTTTACGTTGGGCGGCCAAATAAACCGGACCCGGCGGCGGGGGGGTGAGGCTCATCGTCTCCCCCGGCCGCTTCGCCTTCCTCACGGGGCGGGTACTGACGCCGCGCTc                                        >  2:168136/1‑151 (MQ=255)
                                                                                                               cACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACtt                                 <  1:1364379/151‑1 (MQ=255)
                                                                                                                  cGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAg                              <  1:1606523/151‑1 (MQ=255)
                                                                                                                    gACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCt                            <  1:1296426/151‑1 (MQ=255)
                                                                                                                      cgcgTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgc                          <  1:1295609/151‑1 (MQ=255)
                                                                                                                       gcgTCATCGGCGGTACTGGCGGGTGGGGTTTTCGGTGTGCGGCAAAAAAAACCCAACGCGGCTGCGGGGGGGTGAGCTTCATCCTCTCCCGCGTCTGATACGCCCGCCGCACGGGGCTGGTAATGACGCCGTGCCCCTGACTTaaaccggg                         >  1:413072/1‑145 (MQ=255)
                                                                                                                         gTCATCGGCGGTACTGGCGGCTGGGGCTTTCGGTGTGCGGCAAAAAAAACCCAACGCGGCCGCGGGGGGGTGGGGTTCCTCCTCTCCCCCGGCTGATACGCCCGCCGCCCGGGGGTGATACTGACGCCGGGGTCCTGACTTaaaccgggc                        >  1:1171730/1‑143 (MQ=255)
                                                                                                                         gTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTTTGCGGCAAAATAACCCCAAACCGGCTGCGGGGGGGTGAGCTTCATCTTCTCCCGCGGCTGGTACGCCTGCTGCACGGGGCTGGTACTGACGCCGTGCTCCTGACTTAAGCTgcgc                        >  2:1063952/1‑150 (MQ=255)
                                                                                                                           cATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCAc                      <  2:1164362/150‑1 (MQ=255)
                                                                                                                            aTCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCg                    <  2:101771/151‑1 (MQ=255)
                                                                                                                                gcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                >  1:764345/1‑151 (MQ=255)
                                                                                                                                gcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                <  2:1510826/151‑1 (MQ=255)
                                                                                                                                gcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                <  2:873516/151‑1 (MQ=255)
                                                                                                                                   gTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCaa             <  2:750545/151‑1 (MQ=255)
                                                                                                                                       tGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtc         >  1:376622/1‑151 (MQ=255)
                                                                                                                                        ggcggcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtct        <  2:1064569/151‑1 (MQ=255)
                                                                                                                                        ggcggcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtc         <  2:414788/150‑1 (MQ=255)
                                                                                                                                         gcggcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtctc       <  2:22930/151‑1 (MQ=255)
                                                                                                                                         gcggcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGAAACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtctc       >  2:277598/1‑151 (MQ=255)
                                                                                                                                            gcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTcccc    <  2:421746/151‑1 (MQ=255)
                                                                                                                                               gggCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGt  <  2:1210813/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGT  >  NC_000913/4571218‑4571510

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