Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,531,8290AG45.3% 12.1 / 82.9 64T123A (ACC→GCC) cysZsulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base G (14/15);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGC  >  NC_000913/2531690‑2531964
                                                                                                                                           |                                                                                                                                       
gtgGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATa                                                                                                                                         >  2:425917/1‑140 (MQ=255)
   gCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCg                                                                                                                                      >  1:1162787/1‑140 (MQ=255)
        tGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGTGCTACACCGCCAGATGCCGGGAtt                                                                                                                                 <  1:679438/140‑1 (MQ=255)
        tGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCgg                                                                                                                                     <  2:1132401/136‑1 (MQ=255)
        tGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCgg                                                                                                                                     >  1:1132401/1‑136 (MQ=255)
        tGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGAtt                                                                                                                                 >  1:1206692/1‑140 (MQ=255)
              tCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGt                                                                                                                           <  2:284183/140‑1 (MQ=255)
              tCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGt                                                                                                                           <  1:300067/140‑1 (MQ=255)
              tCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGt                                                                                                                           <  2:98042/140‑1 (MQ=255)
              tCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGt                                                                                                                           <  1:886941/140‑1 (MQ=255)
                aTCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTTGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTAt                                                                                                                         >  1:190579/1‑140 (MQ=255)
                aTCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTAt                                                                                                                         <  2:98780/140‑1 (MQ=255)
                       tgCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGaaa                                                                                                                  >  2:812810/1‑140 (MQ=255)
                            ttAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGAtgt                                                                                                             >  1:1093387/1‑140 (MQ=255)
                             tAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGAtgtg                                                                                                            >  2:757177/1‑140 (MQ=255)
                             tAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGt                                                                                                                           <  2:54833/125‑1 (MQ=255)
                             tAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGt                                                                                                                           >  1:54833/1‑125 (MQ=255)
                                   ttGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGa                                                                                                      <  1:757177/140‑1 (MQ=255)
                                    tGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCcg                                                                                                         >  1:218505/1‑136 (MQ=255)
                                    tGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCcg                                                                                                         <  2:218505/136‑1 (MQ=255)
                                       cTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCa                                                                                                  <  1:248111/140‑1 (MQ=255)
                                                  ccACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGaa                                                                                       <  1:1030153/140‑1 (MQ=255)
                                                    aCGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATg                                                                                     >  1:1030319/1‑140 (MQ=255)
                                                      gATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGc                                                                                   <  2:470212/140‑1 (MQ=255)
                                                      gATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGc                                                                                   >  2:381691/1‑140 (MQ=255)
                                                        ttGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCaa                                                                                 >  1:56083/1‑140 (MQ=255)
                                                        ttGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCaa                                                                                 >  2:961380/1‑140 (MQ=255)
                                                             aaCTGGATTGCCGCTCCGTTTAACGGTTTATTGTCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAAtt                                                                            <  2:17300/140‑1 (MQ=255)
                                                                tGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGc                                                                         <  1:746325/140‑1 (MQ=255)
                                                                tGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGc                                                                         <  2:22086/140‑1 (MQ=255)
                                                                tGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGc                                                                         <  1:369873/140‑1 (MQ=255)
                                                                  gATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCt                                                                       >  1:649902/1‑140 (MQ=255)
                                                                       cctcTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTGGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTAt                                                                  >  2:1000170/4‑140 (MQ=255)
                                                                        cGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAAttt                                                                           >  2:786754/1‑130 (MQ=255)
                                                                        cGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAAttt                                                                           <  1:786754/130‑1 (MQ=255)
                                                                        cGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATc                                                                 >  2:592373/1‑140 (MQ=255)
                                                                          cTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTg                                                               >  1:1071386/1‑140 (MQ=255)
                                                                          cTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTg                                                               >  2:1174138/1‑140 (MQ=255)
                                                                            ccGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGcc                                                             <  2:273884/140‑1 (MQ=255)
                                                                              gTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCcgc                                                           <  1:1174138/140‑1 (MQ=255)
                                                                               tttAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCcgcg                                                          <  2:1181041/140‑1 (MQ=255)
                                                                               tttAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCcgcg                                                          <  2:618885/140‑1 (MQ=255)
                                                                               tttAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCcgcg                                                          <  2:1206692/140‑1 (MQ=255)
                                                                                ttAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCcgcgc                                                         <  1:961380/140‑1 (MQ=255)
                                                                                    cGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCaa                                                     <  2:184133/140‑1 (MQ=255)
                                                                                      gTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAAtt                                                   <  2:755578/140‑1 (MQ=255)
                                                                                      gTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAAtt                                                   <  1:381691/140‑1 (MQ=255)
                                                                                         tATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTa                                                <  2:56083/140‑1 (MQ=255)
                                                                                         tATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTa                                                >  2:1002741/1‑140 (MQ=255)
                                                                                              gCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACt                                           >  2:141707/1‑140 (MQ=255)
                                                                                              gCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACt                                           >  2:812289/1‑140 (MQ=255)
                                                                                               cTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTa                                          <  2:1093387/140‑1 (MQ=255)
                                                                                                 gAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAAt                                        >  1:917032/1‑140 (MQ=255)
                                                                                                  aacaacTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAAtt                                       >  2:576036/1‑140 (MQ=255)
                                                                                                        tGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGt                                                 <  1:1027005/124‑1 (MQ=255)
                                                                                                        tGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGt                                                 >  2:1027005/1‑124 (MQ=255)
                                                                                                           aaGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCTTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTc                              <  2:1162787/140‑1 (MQ=255)
                                                                                                                cGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATccc                         <  2:702922/140‑1 (MQ=255)
                                                                                                                cGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATccc                         <  2:1071386/140‑1 (MQ=255)
                                                                                                                           cgcTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGc              >  1:187182/1‑140 (MQ=255)
                                                                                                                            gcTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCa             <  1:1213752/140‑1 (MQ=255)
                                                                                                                             cTACACCGCCAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCaa            >  2:670438/1‑140 (MQ=255)
                                                                                                                                  ccGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGt       <  1:865175/140‑1 (MQ=255)
                                                                                                                                    gCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGt       >  2:284941/1‑138 (MQ=255)
                                                                                                                                    gCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGt       <  1:284941/138‑1 (MQ=255)
                                                                                                                                     ccAGATGCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTcgc    >  1:870948/1‑140 (MQ=255)
                                                                                                                                           aCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCgcgc  >  1:1163064/1‑136 (MQ=255)
                                                                                                                                           aCCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCgcgc  <  2:1163064/136‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGCTATTTCTTCTCCACGATTGCTAACTGGATTGCCGCTCCGTTTAACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGC  >  NC_000913/2531690‑2531964

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