Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9940GT35.0% 204.6 / 36.4 191intergenic (+392/+631)yhaC/garKpentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base G (16/90);  new base T (4/53);  total (20/171)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

GAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAG  >  NC_000913/3269845‑3270142
                                                                                                                                                     |                                                                                                                                                     
 aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                                                     <  5:312060/151‑1 (MQ=255)
   aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                                                                   <  7:846145/151‑1 (MQ=255)
    cccGCCGGACCGGGTTATTACGTACATCAGGTTAAACTGATCGATAAGCCTCTTTCTTTTGGGTATAGTGTCGTGTACAGTCATTCATCTTTCTGCCCCTGCATAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                                                  <  8:104509/151‑1 (MQ=11)
    cacgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                                                  <  5:1319848/149‑1 (MQ=255)
     ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTg                                                                                                                                                 <  7:1089439/151‑1 (MQ=255)
      cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                                                                                <  6:97858/151‑1 (MQ=255)
      cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                                                                                <  6:1108763/151‑1 (MQ=255)
      cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                                                                                <  8:1447197/151‑1 (MQ=255)
       gccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                                                               <  8:681263/151‑1 (MQ=255)
        ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                                              <  7:43914/151‑1 (MQ=255)
         cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAACATCTTTATTATCGTAATATAGTCATAGACAGTAATTAATATTTAAGAACAACCAAAACAAAAAAACATAAAAAAAAAGTTTTTAAGTAAAACAGGTGAAACTGACCg                                                                                                                                             >  8:1214527/1‑151 (MQ=255)
          ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCAATTTCTTTTTGTTATATAATCGTGTACACTAATTAATATTTATGCCCCTACAAAATAAAAAAAAAAAAAAAAAGGGTTTTTACGTAAATCAGGTGAAACTGACCgc                                                                                                                                            >  2:1142374/1‑150 (MQ=9)
          ggAACGGTTTATTACGTAAATCAGGTAAAACTGACCGATAATCCGATTACTTTTGGGTATAGAATCGTGGAAAGAAATTAATAATTAATCCCATCAAAAACAAAAAAAAAAACAAAACCGGTTTTTACGTAAAACAGGTGAAAATGACCg                                                                                                                                             >  2:1062402/1‑150 (MQ=255)
            aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                                                          <  7:211789/151‑1 (MQ=255)
            aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                                                          <  1:1275217/151‑1 (MQ=255)
             aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                                                         <  1:465378/151‑1 (MQ=255)
             aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                                                         <  4:1686416/151‑1 (MQ=255)
                gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                                                                                                                                      <  7:947628/151‑1 (MQ=255)
                   tattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                                                                                                                                    <  6:569646/150‑1 (MQ=255)
                    attatttACATTAGGTAAAACTGTTCGATTAGCCGATTTTTTTTGGTTATTTTGTCGTGGACTGTCTTTCTTATTTATGGCACTCCAAAATTTTAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGct                                                                                                                                   <  6:1021755/144‑1 (MQ=2)
                    attaCGGACATCAGGGAACACTGACCGATCAGCCGCGGTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctt                                                                                                                                  <  3:1350409/151‑1 (MQ=14)
                      taCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTGGTGGCGAGTCATTCATCGTTATGCGACTCCAAAACGAAAAAACCCCAGAATCGGGTTTGTACGTATAAAAGCGTCAACTGAGCGATAATCCGctttc                                                                                                                                >  3:1426878/1‑151 (MQ=9)
                      taCGTACATCAGGTAAAAATGACCGATAAGCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcttt                                                                                                                                 <  3:662345/150‑1 (MQ=255)
                       aCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                                               <  4:584506/151‑1 (MQ=40)
                       aCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCCGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                                               <  5:410412/151‑1 (MQ=255)
                       aCGTACACCAGGTATATCTCAACTATAAGCTTCTTTCTTTTGGTTATAGTCGCTATTACTTTCAGTTATATTTATGCCACTCAAAAATCTTAAACCCGCGGTATCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGctttct                                                                                                                               <  6:355879/151‑1 (MQ=2)
                           acaacaCGTAACACTGACCGATAAGCCGCTTGCTTTTGCCGATCGTGTCGTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGTATAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                                                           <  6:1680746/151‑1 (MQ=11)
                            cATCACTTATAACTTAGCGATAATCCGCTTTCTTTTTTGTATAGTGTCGTGGAAACTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                                          <  8:383631/151‑1 (MQ=11)
                               cagtaaataCTGACCGATAAGGCGCTTTCTTTTGGGTCTAGTGTCGTGGCCAGGGATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                                                                       <  4:1369533/143‑1 (MQ=11)
                                    caaaCTGTCCTATACGCCGGATTGTTTATCGTATCATGTCGGGGACTGCCATACCTATTTCTGCCACTCAAAAATTATAAACCCGCCGTATCGGTTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt                                                                                                                   <  4:165774/149‑1 (MQ=2)
                                    aaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt                                                                                                                   <  5:1448186/150‑1 (MQ=25)
                                          gACCGTTATGCGGCTTTCTTTTGTTTTAAGTGTCGTGGACTGTGTTTCCTCTTTCTGGCACTCAAAAATTATAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                                                             <  6:818584/150‑1 (MQ=11)
                                          gACCGATAAGCCGCATTCATTTGGCTAAAGTGTCGCGGACTGGCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                                              <  4:191582/149‑1 (MQ=11)
                                            ccGATAAGCCGCTTGCTTTTGCGGACGGTGTCGTGGAAAGTCTTTCATCTTTATCGCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAACACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAc                                                                                                          <  8:1284235/151‑1 (MQ=11)
                                                 aaggtgCTTTCTTTTGGGTCTAGTGTCGGGGACAGTCCTTCATATATCTGCCCCTCAAAAAGTTCAAACCCGCCGTAGCGGTTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtca                                                                                                     <  4:221949/146‑1 (MQ=11)
                                                 aaGCCGTTTTCTTTTTAGTATAGTGTCGGGGACAGTCATTCTTCTTTATGGCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                                                                      <  3:1146829/150‑1 (MQ=11)
                                                    gcGCATTCTTTTGGTTATAGTGTCGGGGACAGTCATTCATCTTTCTGGCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                                                                  <  3:1666215/150‑1 (MQ=255)
                                                    ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                                                                  <  3:773416/151‑1 (MQ=25)
                                                          tctttTTGTTTTTGTCTCGTGGAATGTGTTTCATATTTCTGGCACTCCAAAAGTATAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCtt                                                                                             <  1:1239001/150‑1 (MQ=11)
                                                             tttGGGTATAGTGTCGTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                            <  3:302813/148‑1 (MQ=25)
                                                             tgagtcttttgtgtCGTGGAATGTGTTTGATCTTTCTGGCACTCCAAAAGTATAAACCCGCCGAATCGGGTTTTTACGTAAAAGAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCt                                                                                          <  1:264937/140‑1 (MQ=11)
                                                              ttGCGTATAGTGTCGTGGACAGGCACTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTg                                                                                         <  2:1255832/150‑1 (MQ=255)
                                                                ttttattGTGTCGTGGAATGTGATTGATCTTTATGCCACTCAAAAAGTTAAAACCCGCCGTAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcc                                                                                       <  2:1699409/147‑1 (MQ=11)
                                                                         gtcgtcGAAAGTCATTCTTCTTTCTGGCCCTCAAAAAGTAAAAACCCTCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                                             <  2:269401/151‑1 (MQ=11)
                                                                         gtCGAGGACAGAGTTTAATCTTTATGGCAGTCCAAAAGCTCAAACCCGCCGTAGCGGTTTATTACGTAAAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                              <  7:1243597/150‑1 (MQ=11)
                                                                            gcggaaactgtattttCTTTATGGCCCTCAAAAATTTTAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa                                                                          <  6:1578652/136‑1 (MQ=11)
                                                                                 atgcgtctgctatatctcCCACTCAACAAGCAAAAACCCGCCGCAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                     <  5:268769/137‑1 (MQ=11)
                                                                                 atcggtatactatttCTCCCACTCAAAAAGTTAAAACCCGCCGCACCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                     <  3:1011754/139‑1 (MQ=11)
                                                                                  agggttttttattaatGCCACTCAAAAATTTAAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                    <  7:280677/136‑1 (MQ=11)
                                                                                   gtgaatctacTTTCTGCCACTCAAAAATTTAAAACCCGCCGTACCGGGTTTTTACGTAAATCAGGTCATACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                    <  4:377646/141‑1 (MQ=11)
                                                                                   gggtttcttcTTTCTCCCACTCAACATTTATAAACCCGCAGTAGCGGTTTTTTACGTATATCAGGTGTCACTGACCGATAAGCCCCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATCTTTCTGCCCCTGCAAAAGCAAAATccc                                                                    <  6:345366/146‑1 (MQ=2)
                                                                                    ggcgttgtCTTTCTGCCCCGCCACAAGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                                  <  2:867382/144‑1 (MQ=11)
                                                                                      ggtcatCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg                                                                <  7:559006/149‑1 (MQ=25)
                                                                                      attcatCTTTCTGCCCCTCCAAAATTAAAAAAACGCCGAAGAGGTTTTTTACATAAAACAGGTTAAACTAACAGATAATCCTCTTTCTTTTGTTTAAAGAATCGTGGGCAGTGATTGATCTTTGTGGCCCTGCaaaaagataggaaggggg                                                                >  8:1637154/1‑137 (MQ=2)
                                                                                        attaatttCTCCCACTCCAAAAGTATAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                               <  6:466176/145‑1 (MQ=11)
                                                                                          ttatttCTGGCACTCCAAAAGTTAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                             <  5:1054698/147‑1 (MQ=255)
                                                                                              ttCTGGCCCTCAAAAATATGAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑Cggg                                                        <  5:1661362/150‑1 (MQ=11)
                                                                                              ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑Cggg                                                        <  3:247387/150‑1 (MQ=25)
                                                                                              ttCTCCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCCCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑c                                                           <  2:645046/147‑1 (MQ=11)
                                                                                               gcgcccccTCCACAAGCACAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGtt                                                      <  3:1263280/147‑1 (MQ=255)
                                                                                                cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTAGAGTGTCGTCGAAAGTCAATCAAAATTTAACACCCACAAAAACAAAAACCAGCAGAAgaacgggg                                                      >  2:402272/1‑144 (MQ=2)
                                                                                                  gtcagtcaaaaaTTATAAACCCTCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGtttt                                                    <  2:340644/142‑1 (MQ=11)
                                                                                                     ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGTTTCTTTTGGGTATAGTGTCGTGGGCCGTCAGTGCTGTTTGGGCGCCTCCAAAAGCCAAAACCGCCCCCCG‑CGGGttttttcc                                                >  7:836294/1‑148 (MQ=11)
                                                                                                        ccACAAGCACAAACCCGCCGCAGCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTa                                              <  5:563448/150‑1 (MQ=255)
                                                                                                        ccAAAAGTTAAAACCCGCCGGATCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCACTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGt                                               <  1:847672/149‑1 (MQ=255)
                                                                                                        ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACACGTGAAACTGACCGATAAGCCGCTGTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTaa                                             <  3:1501148/151‑1 (MQ=9)
                                                                                                        caaaaaTTATAAACCCGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTATGTaa                                             <  2:1550818/149‑1 (MQ=11)
                                                                                                        caaaaaGCAAAAACCCGCCGAATCGTGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTaa                                             <  2:1463142/149‑1 (MQ=255)
                                                                                                          aaaaGCAAAAACCCGCCGCACCGGGTTATTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAAc                                           <  4:1245747/151‑1 (MQ=255)
                                                                                                           aaatttaaaaCCCGCCGCATCGGTTTTTTAAGTAAAACAGGTTAAACTTACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAAcc                                          <  8:366597/145‑1 (MQ=11)
                                                                                                           aaaTTAAAAACCCGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAAcc                                          <  2:1113775/151‑1 (MQ=255)
                                                                                                            aattataaaCCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCa                                         <  8:43242/145‑1 (MQ=11)
                                                                                                                ttaaaCCGGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGt                                      <  1:192766/148‑1 (MQ=11)
                                                                                                                 taaaCCCGCCGCATCGGTTTTTTACGTAAATCAGTTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTg                                     <  8:872835/149‑1 (MQ=11)
                                                                                                                  aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaa                                   <  3:936289/151‑1 (MQ=35)
                                                                                                                  aaaCCAGCCGTTTCGGTTTTTTACGTATATCTGGTTAAACTGTGCGATAAGCTGCTTTCTTTTGTGTATTGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGa                                    <  1:1074216/150‑1 (MQ=11)
                                                                                                                   aaCCCGCCGTAGCGGTTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaaa                                  <  3:792150/151‑1 (MQ=255)
                                                                                                                   aaCCCGCCGGATCGGTTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaaa                                  <  2:784411/151‑1 (MQ=255)
                                                                                                                   aaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaaa                                  <  5:618278/151‑1 (MQ=25)
                                                                                                                   aaCCCGCCGAATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaaa                                  <  1:752045/151‑1 (MQ=255)
                                                                                                                     cccGTCGTACCGGTTTTTTACGTAAATCAGGTTATACTGACCGATAAGCCGCTTTCATTTGGCTATAGTGTCGTGGACAGTCTTTCATCTTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACt                                <  6:233715/151‑1 (MQ=11)
                                                                                                                       cgccgAAGCGGTTTTTTACGTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGa                              <  7:1599625/151‑1 (MQ=255)
                                                                                                                       cgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGa                              <  3:518559/151‑1 (MQ=35)
                                                                                                                         ccgtatcGGTTTTTTACTTAAAACATGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                            <  1:383721/145‑1 (MQ=11)
                                                                                                                         ccgtaccGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                            <  1:1590962/145‑1 (MQ=255)
                                                                                                                         ccgtaacGGTTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                            <  4:565415/145‑1 (MQ=255)
                                                                                                                         ccgCAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                            <  1:10122/151‑1 (MQ=25)
                                                                                                                         ccgAAGCGGTTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                            <  1:1663707/151‑1 (MQ=25)
                                                                                                                          cgcaGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCg                           <  5:816963/148‑1 (MQ=16)
                                                                                                                          cgcaGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCg                           <  3:410201/148‑1 (MQ=18)
                                                                                                                              ccGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATaa                       <  7:692009/150‑1 (MQ=25)
                                                                                                                               cGGTTTTTTACGTAAAAAATTTTAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAT‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                      <  1:963681/151‑1 (MQ=11)
                                                                                                                                gggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAGAGTCAATCAAATTATATCGCCAACAAAAACAAAAACCCGCCGCAGTCCGGTTTTTACTTAATCAAGGCCAAATTG‑CCGATaatc                     >  3:808360/1‑149 (MQ=2)
                                                                                                                                  ttttttACTTAAAACATTTTAAACTTACCGATAAGCCGCTTTCTTTTGTGTATTGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAT‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                    <  2:901907/149‑1 (MQ=11)
                                                                                                                                  tttattACGTAAAACAGGTAAAAATTATCGATAAGCCTCTTTCTTTTGTGTATTGTGTTGTGTTCTGTGTTTTATCTTTCTGCCCCTCAAAAAGCATAAACCCGCCGTAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                    <  1:1521637/146‑1 (MQ=11)
                                                                                                                                  gTTATTACGTAAATGAGGTGAACCTGACCGATCAGCGGCTTTCTTTTTGGTATCGTGTCGGCGACTGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                   <  4:826714/151‑1 (MQ=11)
                                                                                                                                   ttattACGTAGAACAGGTGACACTGACCGATAAGCTGCTTGCTTTTGGGTATAGTGTCGTCGCCAGTGATTCCTCTTTCTGCCCCTCCAAAAGCAGAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                   <  7:658104/147‑1 (MQ=11)
                                                                                                                                             aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGGGGCGAGGCATTCAAGAGTTGGCGCCGACAAACGCAAAAACCCGCCGAAG‑CGGGGTGTTTCGTTTACCAGCGGCAAATTGCCGATAAGCCgggtgctgtttg        >  2:1681823/1‑140 (MQ=2)
                                                                                                                                             aaaCACTTGAAACTTACCTATCCGTCGACTTCTTTTTTGTTTTCTGTCGGGGAATGTCTTTTTTATTTATGCCACTCAAAAATCATAAACCCGCCGAAG‑CGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGggttctgtcgt        <  2:527438/151‑6 (MQ=9)
                                                                                                                                                   gTGACACTGCCCGAGACGCCGCTGGCTTTAGGCTATCGTGTCGTGGAAAGTCAGTCATCTTTCTGCCACTCCAAAAGCAAAAACCCGCCGCAG‑CGGGTTTTTACGTAAACCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg  <  4:1569060/151‑1 (MQ=11)
                                                                                                                                                     |                                                                                                                                                     
GAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAG  >  NC_000913/3269845‑3270142

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