Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,105,4990AC40.7% 259.0 / 12.9 118Y356S (TAC→TCC) mltCmembrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base A (65/5);  new base C (0/48);  total (65/53)
Fisher's exact test for biased strand distribution p-value = 9.11e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCC  >  NC_000913/3105350‑3105647
                                                                                                                                                       |                                                                                                                                                    
    cGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAccgcc                                                                                                                                                 >  2:1287253/1‑151 (MQ=255)
    cGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAccgcc                                                                                                                                                 <  2:278405/151‑1 (MQ=255)
     gAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACCCCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAccgccg                                                                                                                                                <  3:197433/151‑1 (MQ=255)
      aGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTAccgccgc                                                                                                                                               >  3:1399141/1‑151 (MQ=255)
        ttttttCGAATGATAAGATTCGGGCTGCAAATATTATTAACCCCATGCCGCCGGGCGATGTTTATCAGCCGCTGACGCCCCGCCACCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTA‑AATCCCGCGCAAAAATCCTAccgccgc                                                                                                                                               <  3:1050126/147‑1 (MQ=255)
          ttttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGAt                                                                                                                                           >  1:190919/1‑151 (MQ=255)
           tttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGAt                                                                                                                                           <  1:61909/150‑1 (MQ=255)
            ttCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATa                                                                                                                                          >  2:1210479/1‑150 (MQ=255)
             tCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAAt                                                                                                                                        >  3:204930/1‑151 (MQ=255)
               gAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTc                                                                                                                                      >  2:1307331/1‑151 (MQ=255)
                  tGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAt                                                                                                                                   >  2:555883/1‑151 (MQ=255)
                  tGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAt                                                                                                                                   >  3:851659/1‑151 (MQ=255)
                    aTAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCAtt                                                                                                                                  >  2:1011980/1‑150 (MQ=255)
                     tAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTa                                                                                                                                 >  3:489440/1‑150 (MQ=255)
                       aGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAcc                                                                                                                              >  3:382056/1‑151 (MQ=255)
                       aGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAAcc                                                                                                                              >  3:394115/1‑151 (MQ=255)
                           tCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGccc                                                                                                                          >  2:1039202/1‑151 (MQ=255)
                            cAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGcccc                                                                                                                         >  4:21202/1‑151 (MQ=255)
                              ggCTGCAAATATTATTAACCCCATGACGCCGGGCGATTTTTATCAGACGCTGACGACCCGCCATCCCTCGGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                                       <  1:283988/151‑1 (MQ=255)
                               gCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTg                                                                                                                       >  2:1107839/1‑150 (MQ=255)
                                 tGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACg                                                                                                                    >  4:454915/1‑151 (MQ=255)
                                   ccAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGAt                                                                                                                  >  1:155467/1‑151 (MQ=255)
                                    cAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑ACTACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATg                                                                                                                 >  4:365040/1‑151 (MQ=255)
                                     aaTATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGc                                                                                                                >  2:1137747/1‑151 (MQ=255)
                                      atatTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCt                                                                                                               >  4:36637/1‑151 (MQ=255)
                                        attattAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                              >  2:187600/1‑150 (MQ=255)
                                        attattAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                                                              >  3:423807/1‑150 (MQ=255)
                                         ttattaCCCCCATGCCGCGGGGCGTTTTTTACCAGCCGCTGACGCCCCCCCATCCCTCGGCGAATTCTCG‑CCTTTTTTTTTAAAAAGTA‑AATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCt                                                                                                               <  2:97895/148‑1 (MQ=255)
                                             aacccccAGGACGCCGGGCGATTTTAAACAGACGATGCCGCCCCCCCCCCCCCCTGGGAATTCCCC‑CCTTTATTTTTAAAAAGGA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATCCAATTAACCGCCCCTGACGATGCTCAgggg                                                                                                         <  2:698182/145‑1 (MQ=255)
                                              acccccATGCCGCCGGGCGATTTTTATCAGACGCTGCCGCCCCCCCACCCCTCGGCGGAATCTCG‑CCTTTATTTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                                                                           <  2:1231672/143‑1 (MQ=255)
                                                  ccaggccgcgGGGGGATTTTTACCGGACGCTGACGCCCCCCCACCCCCCGGCGAAATCCCC‑CCGTTATCTTTAAAAAGTA‑AATCCCGCACAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaaa                                                                                                   <  2:18324/145‑1 (MQ=255)
                                                  ccATGCCGCGGGGCGATTTTTATCAGACGCTGACGCCCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTA‑AATACCGCGCAAAAATCCTCCCGCCGCCGATATTTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaaa                                                                                                   <  1:344074/151‑1 (MQ=255)
                                                  ccATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaa                                                                                                    >  2:807926/1‑150 (MQ=255)
                                                  caatgccccggggCGTTTTTTATCAGCGGCGGCCGCCCCCCCCCCCCTCGGGGAATTCCCG‑CCTTTATTTTTATAAGGGG‑AATCCCGCAAAAAAACCCCCCCCCCGCAGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCaaaaa                                                                                                   <  3:823801/141‑1 (MQ=255)
                                                   cATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAAt                                                                                                  >  3:946383/1‑151 (MQ=255)
                                                    agaaCGCGGGGCGATGTTAACCAGCCCTTGAGCCCCCCCCCCCCCTCGGCGGAATCTCC‑CCGTTATCTTTATAAAGTA‑AATCCCGCGCAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                 <  3:486306/148‑1 (MQ=255)
                                                    aTGCCGCGGGGCGTTTTTTACCAGCCGCTGACGCCCCGCCATCCCTCGGGGGACTCCCG‑CCTTTATTTTTATAAAGTG‑AACCCCGCGAAAAAATCCTCCCGCCGCAGAAAATTCCATTACCCGCCCCTGACGATGCTCAGGGGCAAAAAt                                                                                                  <  4:452979/150‑1 (MQ=255)
                                                    aTGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                                                                 >  2:543912/1‑151 (MQ=255)
                                                          cggggggttttttACCAGCCGCTGCCGCCCCCCCATCCCTCTGCGAAATCTCG‑CCGTTATCTTTATAAAGTG‑AATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTa                                                                                              <  4:160467/138‑1 (MQ=255)
                                                          ccGGGCGATTTTTACCAGCCCCTGCCCCCCCCCCACCCCTCTGCGAAATCTCG‑CCTTTATTTTTAAAAAGTG‑AATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATc                                                                                            <  2:292242/150‑1 (MQ=255)
                                                            gggCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                                                                         >  1:34899/1‑151 (MQ=255)
                                                              gCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACAt                                                                                       >  1:461901/1‑151 (MQ=255)
                                                                 aTGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcac                                                                                    >  2:1325659/1‑151 (MQ=255)
                                                                  tttttACCAGCCGCTGACCCCCCCCCACCCCTCTGGGAAATCCCG‑CCTTTATTTTTAAAAGGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATca                                                                                     <  3:1361788/147‑1 (MQ=255)
                                                                   gTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATcaca                                                                                   >  3:1044041/1‑150 (MQ=255)
                                                                    tttACCAGACGTTGACACCCCCCCTCCCCTCGGCGAATTCCCG‑CCGTTTTTTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACaa                                                                                  <  3:250542/150‑1 (MQ=255)
                                                                      taccAGCCGCTGCCGCCCCCCCACCCCCCTGCGAATTCTCC‑CGGTTATTTTTAAAAAGTA‑AATCCCGCCCAAAAACCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                                                                <  2:1215402/147‑1 (MQ=255)
                                                                      tATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                               >  3:364714/1‑151 (MQ=255)
                                                                       accAGCCGTTGACGCCCCCCCACCCCTCGGCGAATTCCCG‑CCGTTATTTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAttt                                                                               <  2:14667/148‑1 (MQ=255)
                                                                       aTCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CTGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCACCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                                              >  4:402082/1‑151 (MQ=255)
                                                                         cAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                            >  3:831129/1‑151 (MQ=255)
                                                                         cAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                            >  3:838435/1‑151 (MQ=255)
                                                                              gcatccaccccGCCACCCCTCTCCGGAATCTCG‑CCGTTATCTTTAAAAAGTA‑AATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTg                                                                        <  1:310941/142‑1 (MQ=255)
                                                                              gCTGACGCCCCGCCATCCCTCTGCGGAATCCCC‑CCGTTATCTTTAAAAAGTG‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTc                                                                              <  2:1362155/144‑1 (MQ=255)
                                                                                tGACGCCCCCCCTCCCCTCGGCGAATTCCCC‑CCTTTTTTTTTAAAAAGTG‑AATCCCGCAAAAAATCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGt                                                                            <  3:1107892/144‑1 (MQ=255)
                                                                                 gacgacCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaa                                                                     >  2:750462/1‑150 (MQ=255)
                                                                                  acgacCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAt                                                                   >  2:1401620/1‑151 (MQ=255)
                                                                                    gacCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                                                                 >  3:1193324/1‑151 (MQ=255)
                                                                                     acCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTgg                                                                >  2:1197090/1‑151 (MQ=255)
                                                                                      cccccctcccccCTGCGAATTCCCG‑CCTTTATTTTTAAAAAGTA‑AATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTg                                                                        <  3:58798/134‑1 (MQ=255)
                                                                                      cccccctccccTCTGGGAAACCCCG‑CCGTTATCTTTAAAAAGTA‑AATCCCGCGAAAAAATCCTCCCCCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTggg                                                               <  2:1309510/143‑1 (MQ=255)
                                                                                      cccGCCACCCTACTGAGCAACCTCT‑CCGTTATCTTTACAAAGGA‑AAGCCCGCGAAAAACTCCACCCGCCCGCGATATTGCTATTACCCGCCCCTGATGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                                                                      <  2:851623/144‑1 (MQ=255)
                                                                                      cccGCCACCCCTCTGGGAAATCTCG‑CCTTTATCTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                                                                           <  2:332705/139‑1 (MQ=255)
                                                                                        cccccccccTCGGGGAAATCCCG‑CGTTTTTTTTTAAAAAGTA‑AACCCCGCAAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGaa                                                             <  3:977549/145‑1 (MQ=255)
                                                                                         gCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                                                            >  1:213790/1‑151 (MQ=255)
                                                                                          cccccccTCGGGGAAATCCCG‑CCTTTATTTTAAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATg                                                           <  3:850215/147‑1 (MQ=255)
                                                                                             cccccCTGCGAATTCCCC‑CTGTTTTTTTTAAAAAGTG‑AATCCCGCAAAAAACCCCTCCCCCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGttt                                                        <  3:1287321/150‑1 (MQ=255)
                                                                                              cccTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                       >  3:533220/1‑151 (MQ=255)
                                                                                              cccTCTG‑GGAATCCCCCCCGTTATCTTTAAAAAGTA‑AATCCCGCGAAAAAATCCTCCCGCCGCCAATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGttt                                                        <  2:561650/150‑1 (MQ=255)
                                                                                              cccTCGGGGGAATCCCC‑CCGTTATCTTTAAAAAGTT‑AAACCCGCGAAAAAACCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg                                                       <  2:423777/151‑1 (MQ=255)
                                                                                               cccctgggaaaTCCCC‑CCTTTATTTTTAAAAAGTG‑AAACCCGCAAAAAACCCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                      <  3:1255451/142‑1 (MQ=255)
                                                                                                ctctGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc                                                      >  2:231781/1‑150 (MQ=255)
                                                                                                    gCGGAATCTCG‑CCGTTATTTTTAAAAAGTG‑AATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAat                                                   <  1:392752/149‑1 (MQ=255)
                                                                                                      ggAATCCCG‑CCTTTTTTTTTAAAAAGTA‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat                                                <  2:394083/150‑1 (MQ=255)
                                                                                                          tctcG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGcc                                           >  3:174471/1‑151 (MQ=255)
                                                                                                                cctttATTTTTAAAAAGTG‑AATCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAg                                       <  4:136719/147‑1 (MQ=255)
                                                                                                                ccGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAgg                                      >  3:1035433/1‑151 (MQ=255)
                                                                                                                ccGTTACCTTTAAAAAGAA‑AATCCCCCACAAAATCCCTCCCCCCCCCGATAATTCCATAAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAg                                       <  2:849535/150‑1 (MQ=255)
                                                                                                                 cGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGt                                     >  3:474370/1‑151 (MQ=255)
                                                                                                                  gTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGt                                     >  2:327660/1‑150 (MQ=255)
                                                                                                                   ttATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTa                                    >  1:388159/1‑150 (MQ=255)
                                                                                                                    tATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                  >  2:1187220/1‑151 (MQ=255)
                                                                                                                    tATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCACCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                  >  2:602168/1‑151 (MQ=255)
                                                                                                                     aTCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAAc                                  >  2:1404871/1‑150 (MQ=255)
                                                                                                                         ttATAAAGTG‑AATACCCCGCAAAAACCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaaa                              <  1:325579/150‑1 (MQ=255)
                                                                                                                         ttATAAAGTA‑AATCCCGCGCAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACaaa                               <  1:377420/149‑1 (MQ=255)
                                                                                                                            taAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                          >  3:320206/1‑151 (MQ=255)
                                                                                                                            aaaaGTGAAATCCCCCGAAAAATCCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAc                           <  2:500906/150‑1 (MQ=255)
                                                                                                                             aaaGTG‑AATCCCGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCa                         <  1:172898/151‑1 (MQ=255)
                                                                                                                              aaGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAAcc                          >  2:482880/1‑149 (MQ=255)
                                                                                                                                   aaaTCCCGCAAAAAAACCCTCCCGCCGCCGATAATCCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcg                    <  2:204901/150‑1 (MQ=255)
                                                                                                                                    aaTCCCGCGAAAAAACCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTcc                     <  2:382020/149‑1 (MQ=255)
                                                                                                                                    aaTACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgc                   >  3:1144101/1‑151 (MQ=255)
                                                                                                                                       aCCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCcgcg                  >  2:391639/1‑149 (MQ=255)
                                                                                                                                        ccGCGAAAAAATCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGa                 <  1:442104/149‑1 (MQ=255)
                                                                                                                                         cgcgCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAg               >  2:1141780/1‑150 (MQ=255)
                                                                                                                                           cgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTg            >  3:535436/1‑151 (MQ=255)
                                                                                                                                              aaaaaTCCTCCCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATa         <  1:325837/151‑1 (MQ=255)
                                                                                                                                               aaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGAt          >  1:317869/1‑149 (MQ=255)
                                                                                                                                                  aTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAAt     >  3:1218633/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GCTGCGAGTCTTTTCGAATGATAAGATTCAGGCTGCCAATATTATTAACACCATGACGCCGGGCGATGTTTATCAGACGCTGACGACCCGCCATCCCTCTGCGGAATCTCG‑CCGTTATCTTTATAAAGTG‑AATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTTGCCACAGGTAACAAAAAACCAGTCCGCGAAGTTGATAGAATCCC  >  NC_000913/3105350‑3105647

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