Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,465,1670GT31.1% 67.3 / 17.5 46C29F (TGT→TTT) gspMType II secretion system protein GspM
Reads supporting (aligned to +/- strand):  ref base G (9/22);  new base T (14/0);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.81e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCC  >  NC_000913/3464923‑3465412
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                               
gAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACGCAGGGAAATGATTAAATCATGGTGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCTGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑t                                                                                                                                                                                                                                                         <  2:569851/251‑1 (MQ=255)
 aaCTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑tt                                                                                                                                                                                                                                                        <  1:241718/251‑1 (MQ=255)
       aGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGt                                                                                                                                                                                                                                                   <  2:630323/250‑1 (MQ=255)
       aGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCCCGC‑CC‑GGCG‑TTTTTTTTT‑GGG‑TTgggggg                                                                                                                                                                                                                                                  >  1:122105/1‑249 (MQ=255)
               gCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑cc                                                                                                                                                                                                                                         <  1:367240/251‑1 (MQ=255)
                      tCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGGTTTCTCCC‑CC‑CGCG‑GTTTTTTTT‑TGC‑TTGGGGTGGTAAAA‑CCgatacat                                                                                                                                                                                                                                  >  1:577874/1‑248 (MQ=255)
                                tAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACAt                                                                                                                                                                                                                        <  1:339777/249‑1 (MQ=255)
                                  aaaaTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑CGCG‑TTTTTTTTTTGGT‑TTGGGG‑GGTAAAA‑ACCATAG‑ATA‑CCCACATTGcc                                                                                                                                                                                                                    >  1:15492/1‑250 (MQ=255)
                                          tGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑CGCGGTTTTTTTTT‑GG‑‑TTGGGGGGGTAAAA‑ACGATAG‑ATA‑ACTAccattcagaggatcc                                                                                                                                                                                                            >  1:518904/1‑237 (MQ=255)
                                             tAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑GC‑GGCG‑TTTTTTTTT‑TGG‑TTGGGGGGATAAAA‑CCGCTAA‑ATA‑CCTAcctagtcaggcactagtt                                                                                                                                                                                                         >  2:398994/1‑234 (MQ=255)
                                               agagTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑CC‑GGCG‑GTTTTTTTT‑GGT‑TTGGGGGGGTAAAA‑AACATAG‑ATA‑CCTaaaattctgagaaataatccc                                                                                                                                                                                                       >  2:473526/1‑231 (MQ=255)
                                                    gAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTttgtt‑g                                                                                                                                                                                                                                                             <  1:395616/196‑1 (MQ=255)
                                                               cgatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTaa                                                                                                                                                                                     <  1:89547/251‑1 (MQ=255)
                                                                gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCgata                                                                                                                                                                                                                                     <  1:383129/206‑1 (MQ=255)
                                                                           tgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑GTTTTTTTT‑GGG‑TTGGGGGGGGAAAA‑CACATAG‑ATA‑CCTACCTTGCATAGCAGCCG‑CCCCCCACGC‑CCAAAATTAAGAA‑GGACATCa                                                                                                                                                                         >  2:257551/1‑250 (MQ=255)
                                                                            gAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCaa                                                                                                                                                                        <  2:122105/250‑1 (MQ=255)
                                                                                                         aaaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTTTCCC‑CC‑GGCG‑TTTTTTTTT‑GGG‑TTGGGGGGGGAAAA‑ACGATAA‑AAA‑CCTTCCTTTCAGGGGATCAGGTCCC‑CCCGC‑AAAAAAATAAAAA‑AGACCTAAAAAGGATGCAAGATCAAACGAGG‑Accccatg                                                                                                                                         >  2:15801/1‑247 (MQ=255)
                                                                                                          aaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGt                                                                                                                                        <  2:465769/251‑1 (MQ=255)
                                                                                                             aCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGCCTC‑G‑CG‑TTTTTTTTTTGGT‑TTGGGGGGA‑AAAA‑ACGGTAGGTTA‑CC‑ACATTGCAGAGCAACCG‑TCCCACCCCCCAAAAAATAAAAAA‑GGA‑AACAAATAGAAGCAGGAttcatacgggcccccaggggtt                                                                                                                                      >  2:322054/1‑230 (MQ=255)
                                                                                                                gAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTACC‑TC‑GGCG‑GTTTTTGTT‑GGC‑GTGGGGGGGTAAAA‑ACGCTAA‑ATA‑CCTACATTTTAGGGAATCAG‑TCCCACGCGC‑AAAAAATTAAAAA‑AGAAATCAAAAAGAATAAGGATCAAACGAGC‑ACCCAAGGTTTtctc                                                                                                                                  >  2:147129/1‑251 (MQ=255)
                                                                                                                      gATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑t                                                                                                                           <  2:628925/251‑1 (MQ=255)
                                                                                                                              cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCt                                                                                                                   <  2:565027/251‑1 (MQ=255)
                                                                                                                                                         aGGGAAATGATTAACTCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCCCGC‑TC‑GGCG‑GTTTTTTTT‑GGC‑TTTGGGTGGTAAAAACCAAAAC‑A‑A‑CCCACATTGCCGAGCATCAG‑GCCCCCCCGC‑AAAAAATTAAAAA‑GGGAATCAAATGGCTGCAGGATCAAGAGAGC‑CCCCCTGGTTTACTCCGTCA‑TCCTGCGCTTTCGCAACCCCCCaaaaaaattctttt                                                                                        >  1:63241/1‑242 (MQ=255)
                                                                                                                                                          gggAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTtctc                                                                                       <  1:15801/251‑1 (MQ=255)
                                                                                                                                                               aTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCC‑CC‑GGCG‑TTTTTTTTT‑GGC‑TTGGGGGGGTAAAA‑CCGATAG‑ATA‑ACCACATTTCAGAGCATCAA‑TCCCCCGCGC‑CAAAAATTAAAAA‑GGACATCAAAATGATGCAGGATCCAGCGAGC‑CCCCCTGGTTTTCTCCGTCA‑TCCTTCGCCTTCGCAGCCCCTCAAAAATATCCTTCTCGAggg                                                                                  >  1:612267/1‑250 (MQ=255)
                                                                                                                                                                       aTCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCaaaat                                                                          <  2:524086/251‑2 (MQ=255)
                                                                                                                                                                          aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTg                                                                       <  1:472797/251‑1 (MQ=255)
                                                                                                                                                                               gggCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGaa                                                                   <  2:362581/250‑1 (MQ=255)
                                                                                                                                                                                 gCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTCCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGG‑TTGGGGGGGTTAAA‑CCCCTAG‑ATA‑ACTAAATTTTCCAGGATCAG‑TCCCCCGCGC‑CAAAAATTTACAA‑GGGCAACCAATTGATGCCAGCCCAAGCGAGC‑CACCATGGTTTAACCGGGCACTCC‑GCGCTTTCGCCCCCCATCCAAAATATTCTTCTCGCGGAAGccacaaggcaggattt                                                                >  1:410568/1‑236 (MQ=255)
                                                                                                                                                                                    gAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGc                                                                               <  2:238344/233‑1 (MQ=255)
                                                                                                                                                                                    gAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑GTTTTTTTT‑TGC‑TTGGGGGGGTAAAA‑CCGATAG‑AAA‑CCTACATTTCAGAGCATCAG‑TCCCCCGCGC‑AAAAAAATAAAAA‑GGACATCAAATTGATGCAGGATCAAGCGAGC‑ACCCATGGTTTTATCGGTCA‑TCCTGCGCTTACGCAGCCCCTCCAAAATATCCCTCTCGAGGGAGCagatcgggagagcacacc                                                             >  1:238344/1‑234 (MQ=255)
                                                                                                                                                                                          tCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGC‑TC‑GGCG‑GTTTTTGTT‑TGC‑TTGGGGGGGTAAAA‑ACGATAG‑ATA‑CCTACATTTCAGAGCATCAG‑TCCCCCGCGC‑AAAAAAATAAAAA‑GGACATCAAAAAGATGCAGGGTCAAACGAGC‑ACCCCTGGTTTACTCGGGCA‑TCCTGCGCCTACGCCGCCCCTCAAAAAAATCCTTCTCGCGGAAGCAAAACGTGAGAAATTTGCTATCAc                                                       >  1:501411/1‑251 (MQ=255)
                                                                                                                                                                                           ccACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACg                                                      <  1:453700/251‑1 (MQ=255)
                                                                                                                                                                                               aTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑CC‑GGCG‑TTTTTTTTT‑GGCGTTGGGG‑GGTAAAA‑ACGATAG‑ATA‑ACTACATTTCAGAGCATCAG‑TCCCCCGCGC‑AAAAAATTAAGAA‑AGACATCAAATGGATGCAGGAACAAACGAGG‑ACCCCTTGGTTAATCGGGCC‑TCCTGCGCCTAAGCAGCCCCACAAAAAAATCCTTCTCGAGGGAGCAAAAAGAGAGAATTTGGGTATCACGCTgg                                                  >  1:326203/1‑251 (MQ=255)
                                                                                                                                                                                                tCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑CGCG‑GTTTTTTTT‑GGG‑TTGGGGGGGTAAAA‑ACGGTAG‑ATA‑ACTACATTTCAGAGCATCAA‑TCCCCCGCGC‑AAAAAAATAAAAA‑AGACATCAAAAGGAAGCAGGGACAAACGAGCCACCC‑TGGTTTACTCGGTCA‑TCCTGCCGTTACGCCGCCCCACAAAAAAATTCTTCTCCAGGAAGCAAAACGGGAGAATTTTGCTATTACGccggg                                                 >  1:378949/1‑247 (MQ=255)
                                                                                                                                                                                                  ggAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TCCGGGG‑TTTTTT‑TT‑GGC‑TTGGGGTGGGAAAA‑ACGGTAG‑ATA‑CCTCCATTTTAAAGCATCCA‑TCCCCCGCGC‑AAAAAAATAAGAA‑GGGCAACAAAAAGAAGCAGGGTAAAGCGAGC‑ACCCCTGGGTTACTCGGTCA‑TCCCGCGCCTATCCACACCATCAAAAAAATACTTCTCGAGGGAGCAAAACGCGAGAATTTTGCTaatccgccggaa                                                >  2:25522/1‑239 (MQ=255)
                                                                                                                                                                                                     aaaaCAGATCGGTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑ACCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACg                                            <  2:34994/251‑1 (MQ=255)
                                                                                                                                                                                                      aaaCAGATCGTTGCGGCACTGGCGGTTCTTTTCC‑TC‑GGGG‑TTTTTTGTT‑GGT‑TTTGGGGGATAAAA‑ACGCTAG‑ATAACCT‑CATTGCAGAGCATCAA‑GCCCCCCCGC‑CAAAAATTAAAAAAGGA‑ATTAAAAGGAAGCAGGATCAAGCGAGC‑ACCCCCGGTTTTCTTGGTCC‑TCCTCCGCTTACTCCCCCCCTCAAAAAAACCCTCTTCCAGGAAGAAAAACATGGGAAttttggcaataccccgggaaaacg                                           >  2:438375/1‑230 (MQ=255)
                                                                                                                                                                                                                 tGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAAt                                 <  2:15492/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATAct                              <  2:180142/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggTTGTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTc                    >  2:594009/1‑251 (MQ=255)
                                                                                                                                                                                                                                tctctcGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATc                 >  2:374083/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGc       >  1:866/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gcggtttttGGT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGGGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGGCA‑TCCTGCGCTTACGCATCCCATCAAAAAAATCCTTCTCGAGGAAGCAAAACGTGGGAATTTGGCTATCTCGCTGGGAAACGGCCCGGACAAGACTCTGACAAGTCATCCTGTCACAGCt      <  2:243251/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttttttGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTccc   >  1:529133/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tttttGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTcccc  >  1:292903/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ttttGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTcccc  <  2:612267/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       gtt‑gGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATAct                              >  2:352666/1‑218 (MQ=255)
                                                                                                                                                                                                                                                       gtt‑gGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATAct                              <  1:352666/218‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                               
GAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGC‑TC‑GGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAA‑CCGATAG‑ATA‑CCTACATTGCAGAGCATCAG‑TCCCACGCGC‑AAAAAATTAAGAA‑GGACATCAAATGGATGCAGGATCAAGCGAGC‑ACCCATGGTTTACTCGGTCA‑TCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCC  >  NC_000913/3464923‑3465412

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