Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,3680AT43.8% 111.9 / 30.2 74R69R (CGT→CGAyjiRputative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (5/36);  new base T (32/0);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 9.74e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATAC  >  NC_000913/4571221‑4571518
                                                                                                                                                   |                                                                                                                                                       
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:446677/151‑1 (MQ=255)
                                                gCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATc                                                                                                     <  1:1300259/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:873516/151‑1 (MQ=255)
                                                                                                                             gcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                        <  2:1510826/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)
                                                                                                                                                   aCGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATAc  >  2:573039/1‑151 (MQ=255)
                                                                                                                                                   aCGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATAc  >  1:1545908/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
CGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATAC  >  NC_000913/4571221‑4571518

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: