Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5080GT45.5% 51.7 / 30.5 56Q335K (CAA→AAA) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (10/20);  new base T (25/0);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 7.01e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTA  >  NC_000913/366363‑366654
                                                                                                                                                  |                                                                                                                                                      
gAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTgtgtt                                                                                                                                                    >  1:538075/1‑149 (MQ=255)
     aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGgc                                                                                                                                             <  1:1019073/151‑1 (MQ=255)
     aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGgc                                                                                                                                             <  1:97127/151‑1 (MQ=255)
       aaaGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGgcg                                                                                                                                            <  1:768731/150‑1 (MQ=255)
        aaGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGcc                                                                                                                                          <  2:169994/151‑1 (MQ=255)
            gtAAAGCCTGGGGCGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGTTC‑TTGGGCGcctg                                                                                                                                        >  2:767240/1‑147 (MQ=255)
              aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGTA‑TTGGGGGCCAggggg                                                                                                                                    >  2:710184/1‑149 (MQ=255)
              aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGggtg                                                                                                                                    <  1:51546/151‑1 (MQ=255)
                   cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGTA‑TTTGGCGCCAGGGGGGgttt                                                                                                                               >  1:155358/1‑151 (MQ=255)
                     ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTACAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttc                                                                                                                             <  2:415572/151‑1 (MQ=255)
                       ggTGCCTAATGAGTGAGCTAACTCACATTGATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTTGGCGCCCggggggggtttttt                                                                                                                           >  2:936422/1‑143 (MQ=255)
                         tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTGGGCGCCCGGGGGGGttttttttt                                                                                                                         >  2:697953/1‑151 (MQ=255)
                          gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                         <  2:30992/150‑1 (MQ=255)
                          gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                         <  2:27743/150‑1 (MQ=255)
                           ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGCATTTTGGGGCCGGGGGGGTtttttttt                                                                                                                          >  2:259185/1‑149 (MQ=255)
                              aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                    <  1:896561/151‑1 (MQ=255)
                               aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                     >  1:296935/1‑149 (MQ=255)
                                           aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCaa                                                                                                      <  2:341218/151‑1 (MQ=255)
                                             tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGGA‑TTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGA‑CGGGCaa                                                                                                      >  2:859052/1‑149 (MQ=255)
                                                  ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCGTGAGA‑CGGGCAACa                                                                                                    >  2:104448/1‑146 (MQ=255)
                                                  ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGA‑CGGGCAAc                                                                                                     >  2:759837/1‑145 (MQ=255)
                                                      ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGGAACAGCTGAtt                                                                                             >  2:353014/1‑149 (MQ=255)
                                                          gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGGA‑TTGGGGGGCCGGGGGGGTTTTTTTTTTTCCAGGGAGA‑AGGGGGACAGCTGATT‑G‑tcc                                                                                       >  1:748805/1‑149 (MQ=255)
                                                           ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGA‑CGGGGAACAACTGATT‑G‑cccct                                                                                     >  1:567202/1‑148 (MQ=255)
                                                                  cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGc                                                                               <  2:885320/149‑1 (MQ=255)
                                                                   aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGGGGT‑TTTGGCGCCCGGGGGGGTTTTCTTTTTCCCCGTGGGA‑CGGGCAACAGCTGATT‑G‑CCCCTCACcccc                                                                              >  1:195724/1‑149 (MQ=255)
                                                                     tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTG‑TTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGG‑AGGGGAACAGCTGGTT‑G‑CCCTTCACcccc                                                                              >  1:764185/1‑147 (MQ=255)
                                                                      gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑CTTTCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTCACCCAGGAGA‑CGGGGCACAGCCGGTTGG‑CCCTTCccccccgg                                                                            >  1:486255/1‑147 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGTG‑TTGGGGGGCAGGGGGGGTTTTTTTTTTCCCCGGGAGA‑CGGGCAAAACCCGATT‑G‑CCCCTCACCCCCCGGccc                                                                        >  2:781833/1‑150 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGA‑CGGGCAACAGCTGATT‑T‑CCCTTCACCCCCCTGccc                                                                        >  2:224325/1‑150 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGCAAAAGCTGATT‑G‑CCCTTCACCCCCTGGccc                                                                        >  2:221220/1‑150 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCt                                                                       <  2:106325/151‑1 (MQ=255)
                                                                          gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGGA‑TTGGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGGGACAGCTTGTT‑G‑CCCCTCCCCCCCCGGccc                                                                        >  1:915757/1‑148 (MQ=255)
                                                                                    cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGG‑CGGGCAACAACTGATT‑G‑CCCTTCACCGCCCGGCCCCGagagaggt                                                              >  2:878026/1‑146 (MQ=255)
                                                                                    cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGCCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGA‑CGGGGCACAGCTGATT‑G‑CCCTTCACCCCCCGGCCCTGGGAGAGGTg                                                             >  2:713345/1‑149 (MQ=255)
                                                                                     gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGA‑AGGGGAACAACTGATT‑G‑CCCTTCACCCCCCGGCCCCGagagaggt                                                              >  2:244512/1‑145 (MQ=255)
                                                                                        aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGAGA‑CGGGCAAAAGCTGATT‑T‑CCCTTCACCCCCTGGCCCTGAGAGAGGTgcagca                                                        >  2:1116175/1‑150 (MQ=255)
                                                                                         aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGCAAAAGCGGATT‑G‑CCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCaa                                                       >  2:123952/1‑150 (MQ=255)
                                                                                          aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTTACCCGTGAGA‑CGGGGAACAGCTGATT‑G‑CCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCaa                                                       >  2:258877/1‑149 (MQ=255)
                                                                                             tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTTCTAT‑TTTGGCGCCGGGGGGGTTTTT‑TTTTTACCCGTGAGA‑CGTGCACCAGCACCTT‑GACCCTCCACC‑CCTGTTCCTCagagaggtttatcaatccg                                                   >  1:209794/1‑137 (MQ=255)
                                                                                             tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGA‑AGGGCAACAGCTGATT‑G‑GCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCggg                                                  >  1:1025718/1‑150 (MQ=255)
                                                                                              gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                 <  1:218461/151‑1 (MQ=255)
                                                                                                cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGCT‑TTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCCTGGGA‑CGGGGCACAGCCTAAT‑G‑TCCTTCACCCCCCCGCCCCGTGAGAGtgggcgccagccggcc                                                >  1:280624/1‑135 (MQ=255)
                                                                                                cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                               <  1:136403/151‑1 (MQ=255)
                                                                                                     cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCCTGAGA‑CGGGGAACAGCTGATT‑T‑CCCCTCCCCCCCCGGCCCTGAGAGAGGTGCAGCCAACGGTCCAc                                              >  1:563343/1‑147 (MQ=255)
                                                                                                           cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                    <  1:598537/151‑1 (MQ=255)
                                                                                                            aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGGG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGGGA‑CAGGCAAAAACCGATT‑G‑CCCTTCCCCGCCAGGCCCTCGGAGAGGTGCAGCAAGAGGGCCACGCCGGgtt                                      >  2:760125/1‑148 (MQ=255)
                                                                                                            aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGGAAAAGCTGATT‑T‑GCCTTCACCGCCTGGGCCTGAGAGAGTTGCAGCAAACGGGCCACGCCGGGTTTcc                                   >  2:740777/1‑151 (MQ=255)
                                                                                                              tAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTCCCCGTGAGA‑CGGGCAACAACCGATT‑G‑CCCCTTACCGCCCGGCCCTGAGAGAGGTGGAGCAAACGGGCCACGCCGGGTTTcccc                                 >  2:1015672/1‑151 (MQ=255)
                                                                                                                aTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                               <  2:395138/151‑1 (MQ=255)
                                                                                                                  gAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                             <  1:416480/151‑1 (MQ=255)
                                                                                                                   aaTCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                            <  1:1152805/151‑1 (MQ=255)
                                                                                                                       ggCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                        <  1:1119601/151‑1 (MQ=255)
                                                                                                                       ggCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                        <  2:995191/151‑1 (MQ=255)
                                                                                                                         ccAACGCGCGGGGAGAGGCGG‑TTTGGCTC‑ATAGGGGCCCGGGGGGTTTTTCTTTTCATCAGTGAGACCGGGCAA‑AGCGGACT‑G‑GCCTTTTCCGCCTTGCCCCGCGGGAGGTGCAGCAAGCCGGTCACCCTGGGGTGCCCCAGCAGGcaca                       >  2:611104/1‑147 (MQ=255)
                                                                                                                                            gg‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt   >  2:560599/1‑151 (MQ=255)
                                                                                                                                               tttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa  >  1:1105291/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTA  >  NC_000913/366363‑366654

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: ‑