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

GAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTA  >  NC_000913/366363‑366654
                                                                                                                                                  |                                                                                                                                                      
gAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTgtgtt                                                                                                                                                    >  1:538075/1‑149 (MQ=255)
     aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgc                                                                                                                                             <  1:1019073/151‑1 (MQ=255)
     aTAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgc                                                                                                                                             <  1:97127/151‑1 (MQ=255)
       aaaGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgcg                                                                                                                                            <  1:768731/150‑1 (MQ=255)
        aaGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGcc                                                                                                                                          <  2:169994/151‑1 (MQ=255)
            gtAAAGCCTGGGGCGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGTTCTTGGGCGcctg                                                                                                                                        >  2:767240/1‑147 (MQ=255)
              aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTATTGGGGGCCAggggg                                                                                                                                    >  2:710184/1‑149 (MQ=255)
              aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                    <  1:51546/151‑1 (MQ=255)
                   cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGCGCCAGGGGGGgttt                                                                                                                               >  1:155358/1‑151 (MQ=255)
                     ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTACAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                             <  2:415572/151‑1 (MQ=255)
                       ggTGCCTAATGAGTGAGCTAACTCACATTGATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCggggggggtttttt                                                                                                                           >  2:936422/1‑143 (MQ=255)
                         tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTGGGCGCCCGGGGGGGttttttttt                                                                                                                         >  2:697953/1‑151 (MQ=255)
                          gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                         <  2:30992/150‑1 (MQ=255)
                          gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                         <  2:27743/150‑1 (MQ=255)
                           ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGCATTTTGGGGCCGGGGGGGTtttttttt                                                                                                                          >  2:259185/1‑149 (MQ=255)
                              aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                    <  1:896561/151‑1 (MQ=255)
                               aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                     >  1:296935/1‑149 (MQ=255)
                                           aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                      <  2:341218/151‑1 (MQ=255)
                                             tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGGATTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCaa                                                                                                      >  2:859052/1‑149 (MQ=255)
                                                  ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGGGGCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACa                                                                                                    >  2:104448/1‑146 (MQ=255)
                                                  ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAAc                                                                                                     >  2:759837/1‑145 (MQ=255)
                                                      ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGAtt                                                                                             >  2:353014/1‑149 (MQ=255)
                                                          gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTGGGGGGCCGGGGGGGTTTTTTTTTTTCCAGGGAGAAGGGGGACAGCTGATTGtcc                                                                                       >  1:748805/1‑149 (MQ=255)
                                                           ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTGTTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAACAACTGATTGcccct                                                                                     >  1:567202/1‑148 (MQ=255)
                                                                  cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                               <  2:885320/149‑1 (MQ=255)
                                                                   aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGGGGTTTTGGCGCCCGGGGGGGTTTTCTTTTTCCCCGTGGGACGGGCAACAGCTGATTGCCCCTCACcccc                                                                              >  1:195724/1‑149 (MQ=255)
                                                                     tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTGTTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGGAGGGGAACAGCTGGTTGCCCTTCACcccc                                                                              >  1:764185/1‑147 (MQ=255)
                                                                      gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGTATTGGGGGCCCGGGGGGGTTTTTTTTTCACCCAGGAGACGGGGCACAGCCGGTTGGCCCTTCccccccgg                                                                            >  1:486255/1‑147 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTGTTGGGGGGCAGGGGGGGTTTTTTTTTTCCCCGGGAGACGGGCAAAACCCGATTGCCCCTCACCCCCCGGccc                                                                        >  2:781833/1‑150 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCCCCCTGccc                                                                        >  2:224325/1‑150 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCCCCTGGccc                                                                        >  2:221220/1‑150 (MQ=255)
                                                                        ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                       <  2:106325/151‑1 (MQ=255)
                                                                          gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTGGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGGACAGCTTGTTGCCCCTCCCCCCCCGGccc                                                                        >  1:915757/1‑148 (MQ=255)
                                                                                    cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTGGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGGCGGGCAACAACTGATTGCCCTTCACCGCCCGGCCCCGagagaggt                                                              >  2:878026/1‑146 (MQ=255)
                                                                                    cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGCCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGCACAGCTGATTGCCCTTCACCCCCCGGCCCTGGGAGAGGTg                                                             >  2:713345/1‑149 (MQ=255)
                                                                                     gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAACAACTGATTGCCCTTCACCCCCCGGCCCCGagagaggt                                                              >  2:244512/1‑145 (MQ=255)
                                                                                        aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTTCCCTTCACCCCCTGGCCCTGAGAGAGGTgcagca                                                        >  2:1116175/1‑150 (MQ=255)
                                                                                         aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAAAAGCGGATTGCCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCaa                                                       >  2:123952/1‑150 (MQ=255)
                                                                                          aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTCTTTTTACCCGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCaa                                                       >  2:258877/1‑149 (MQ=255)
                                                                                             tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTTCTATTTTGGCGCCGGGGGGGTTTTTTTTTTACCCGTGAGACGTGCACCAGCACCTTGACCCTCCACCCCTGTTCCTCagagaggtttatcaatccg                                                   >  1:209794/1‑137 (MQ=255)
                                                                                             tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGAAGGGCAACAGCTGATTGGCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCggg                                                  >  1:1025718/1‑150 (MQ=255)
                                                                                              gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                 <  1:218461/151‑1 (MQ=255)
                                                                                                cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGCTTTGGGGGCCGGGGGGGTTTTTTTTTTACCCCTGGGACGGGGCACAGCCTAATGTCCTTCACCCCCCCGCCCCGTGAGAGtgggcgccagccggcc                                                >  1:280624/1‑135 (MQ=255)
                                                                                                cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                               <  1:136403/151‑1 (MQ=255)
                                                                                                     cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTACCCCTGAGACGGGGAACAGCTGATTTCCCCTCCCCCCCCGGCCCTGAGAGAGGTGCAGCCAACGGTCCAc                                              >  1:563343/1‑147 (MQ=255)
                                                                                                           cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                    <  1:598537/151‑1 (MQ=255)
                                                                                                            aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGGTTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGGGACAGGCAAAAACCGATTGCCCTTCCCCGCCAGGCCCTCGGAGAGGTGCAGCAAGAGGGCCACGCCGGgtt                                      >  2:760125/1‑148 (MQ=255)
                                                                                                            aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAAAAGCTGATTTGCCTTCACCGCCTGGGCCTGAGAGAGTTGCAGCAAACGGGCCACGCCGGGTTTcc                                   >  2:740777/1‑151 (MQ=255)
                                                                                                              tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTCCCCGTGAGACGGGCAACAACCGATTGCCCCTTACCGCCCGGCCCTGAGAGAGGTGGAGCAAACGGGCCACGCCGGGTTTcccc                                 >  2:1015672/1‑151 (MQ=255)
                                                                                                                aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                               <  2:395138/151‑1 (MQ=255)
                                                                                                                  gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                             <  1:416480/151‑1 (MQ=255)
                                                                                                                   aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                            <  1:1152805/151‑1 (MQ=255)
                                                                                                                       ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                        <  1:1119601/151‑1 (MQ=255)
                                                                                                                       ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                        <  2:995191/151‑1 (MQ=255)
                                                                                                                         ccAACGCGCGGGGAGAGGCGGTTTGGCTCATAGGGGCCCGGGGGGTTTTTCTTTTCATCAGTGAGACCGGGCAAAGCGGACTGGCCTTTTCCGCCTTGCCCCGCGGGAGGTGCAGCAAGCCGGTCACCCTGGGGTGCCCCAGCAGGcaca                       >  2:611104/1‑147 (MQ=255)
                                                                                                                                            ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt   >  2:560599/1‑151 (MQ=255)
                                                                                                                                               tttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa  >  1:1105291/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTA  >  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: