Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,358,1200GA21.1% 75.4 / 15.5 38L1132L (CTG→CTA) gltBglutamate synthase subunit GltB
Reads supporting (aligned to +/- strand):  ref base G (19/11);  new base A (4/4);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTG  >  NC_000913/3357885‑3358349
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
agGCGCTGGTTGCTAACGGCTTGCGTCTCAGGATCCGTTTGCAGGTCGATGGCGGCCTGAAGACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAGCGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACg                                                                                                                                                                                                                        <  2:221359/251‑1 (MQ=255)
       ggTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTa                                                                                                                                                                                                                   <  1:196501/249‑1 (MQ=255)
              aaCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGGGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGtt                                                                                                                                                                                                          >  1:42930/1‑251 (MQ=255)
                            caAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGaaa                                                                                                                                                                                            <  2:214294/251‑1 (MQ=255)
                              aGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAAcc                                                                                                                                                                                          >  2:43433/1‑251 (MQ=255)
                                  ccGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCcgcg                                                                                                                                                                                      <  1:21695/251‑1 (MQ=255)
                                     tttGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGc                                                                                                                                                                                   >  1:175501/1‑251 (MQ=255)
                                      gggCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATAGCATCAAGGCGGCGATTCGCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGGCATCTGAACAACTGCTCAACGGGTGTATCAACGCATTATGATAAAGTGCGTAAGAGCTACTATCACGGCCTGCCATTCAAGGTGATGAAGTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGc                                                                                                                                                                                   <  1:162545/248‑1 (MQ=255)
                                          aGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATg                                                                                                                                                                              <  1:49475/251‑1 (MQ=255)
                                            gTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGc                                                                                                                                                                            <  1:237573/251‑1 (MQ=255)
                                              cGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGcac                                                                                                                                                                          <  2:81262/251‑1 (MQ=255)
                                               gATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTATGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGcac                                                                                                                                                                          >  1:71247/1‑250 (MQ=255)
                                                           aaaaCGGGTGTCGATATCATCACGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGCCAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTaa                                                                                                                                                              >  2:95894/1‑250 (MQ=255)
                                                                          atcatcAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCt                                                                                                                                              <  2:175501/251‑1 (MQ=255)
                                                                               caAGGCGGCGATTCTCGGCGCAGACAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACCACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACGATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGAtt                                                                                                                                          <  1:148593/250‑1 (MQ=255)
                                                                                     ggcgATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAATTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGc                                                                                                                                    >  2:186005/1‑250 (MQ=255)
                                                                                       cgATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGTATCTGATTGTTCGCAcc                                                                                                                                 <  2:149594/251‑1 (MQ=255)
                                                                                              cGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACctgc                                                                                                                          >  1:193166/1‑251 (MQ=255)
                                                                                                 cgcAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACctgctg                                                                                                                        >  2:39117/1‑250 (MQ=255)
                                                                                                                   cggcACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCa                                                                                                     <  1:178655/251‑1 (MQ=255)
                                                                                                                                aTGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCACCCTGAACAACTGCGAACCGAGTGCAGCCACTCAGGATGACCACCTGCGTAAGAATCCCCATCCCAGAATGACATTCAAGGCGACGAATTAATTTGAGTTTATCGCCCGTGAAACCCGCCAGCTGATGGCAACGCTTGGCGCAACACGTCTGGCGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTCTCACCGCCaaaaacc                                                                                         >  1:258934/1‑246 (MQ=255)
                                                                                                                                     ggCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGTACGGTTTCACCGCCAAACAGCAGAAg                                                                                    <  1:101310/250‑1 (MQ=255)
                                                                                                                                         cTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGgc                                                                               >  2:41498/1‑251 (MQ=255)
                                                                                                                                            ggCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCt                                                                            >  1:62180/1‑251 (MQ=255)
                                                                                                                                                gTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTc                                                                         <  1:211971/250‑1 (MQ=255)
                                                                                                                                                      atCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAACCTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAgctg                                                                  >  1:250854/1‑251 (MQ=255)
                                                                                                                                                      atCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAgctg                                                                  >  2:55678/1‑251 (MQ=255)
                                                                                                                                                           cGTATTTGCCATCTGAACAACTGCGCAACGGGTGCAGCAACTCAGGATGACAAACTGAGTACGAAACACTATACCGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGCCCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGAACAAGCTGGCGCTGTCGCAGCTGCTgg                                                              >  1:122078/1‑250 (MQ=255)
                                                                                                                                                              aTTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGCCAAACTGCGTCAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAACCAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGAc                                                          >  1:179356/1‑251 (MQ=255)
                                                                                                                                                                    cATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCg                                                     >  2:257606/1‑250 (MQ=255)
                                                                                                                                                                         gAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGc                                               >  1:213681/1‑251 (MQ=255)
                                                                                                                                                                                     aaCGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAgg                                   >  1:50102/1‑251 (MQ=255)
                                                                                                                                                                                              aGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTAc                           <  1:68034/250‑1 (MQ=255)
                                                                                                                                                                                              aGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCACGTGAAACCCGCGCGCTGATGGAACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAACGAGCTGGACGGTTTCCCCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATACAGGTAAGGCACTCTACt                          >  1:269231/1‑251 (MQ=255)
                                                                                                                                                                                                        gATGACAAACTGCGTAAGAACCACTATCACGGCCTACCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaaa                 >  2:175300/1‑250 (MQ=255)
                                                                                                                                                                                                        gATGAAAACCTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTCCTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGaaa                 >  2:32126/1‑250 (MQ=255)
                                                                                                                                                                                                             cAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAAcc           >  1:196910/1‑251 (MQ=255)
                                                                                                                                                                                                                       aaGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTg  >  1:150222/1‑250 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
AGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGATTGGTCGCACCGACCTGCTGAAAGAGCTGGACGGTTTCACCGCCAAACAGCAGAAGCTGGCGCTGTCGAAGCTGCTGGAGACTGCCGAACCGCATCCAGGTAAGGCACTCTACTGCACCGAAAACAACCCGCCGTTTG  >  NC_000913/3357885‑3358349

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