Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,359,685 A→G 6.9% D163G (GAT→GGT)  gltD → glutamate synthase, 4Fe‑4S protein, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,359,6850AG 6.9% 171.5 / 3.4 58D163G (GAT→GGT) gltDglutamate synthase, 4Fe‑4S protein, small subunit
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base G (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TTGGCGCGGTGACCATCGGCAACATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAA  >  NC_000913/3359538‑3359833
                                                                                                                                                   |                                                                                                                                                    
ttGGCGCGGTGACCATCGGCAACATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGAt                                                                                                                                                     <  2:857376/149‑1 (MQ=255)
       ggTGACCATCGGCAACATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGa                                                                                                                                              <  1:749489/149‑1 (MQ=255)
             cATCGGCAACATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTa                                                                                                                                        >  1:182662/1‑149 (MQ=255)
                  gCAACATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGc                                                                                                                                   >  1:773981/1‑149 (MQ=255)
                    aaCATTGAGCGCTATATCAACGGTAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGt                                                                                                                                 >  1:2790/1‑149 (MQ=255)
                    aaCATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGt                                                                                                                                 >  1:1181614/1‑149 (MQ=255)
                             cgcTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGt                                                                                                                        <  2:845572/149‑1 (MQ=255)
                             cgcTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGt                                                                                                                        >  2:128934/1‑149 (MQ=255)
                              gcTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGtt                                                                                                                       >  2:1368824/1‑149 (MQ=255)
                              gcTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGtt                                                                                                                       <  2:1206287/149‑1 (MQ=255)
                                 atatCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTc                                                                                                                    <  1:1255325/149‑1 (MQ=255)
                                 atatCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTc                                                                                                                    <  1:1147149/149‑1 (MQ=255)
                                 atatCAACGATAAAGCGCTCGAGATGGGCTGGCGTCCGGATATGGCGGGTGTGCAACAGACCGGTAAAAAAGTGGCGAATATCGGCGCAGGCCCGGCAGGTCTGGCGTGGGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTc                                                                                                                    <  1:199923/149‑1 (MQ=255)
                                    tCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGCAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTc                                                                                                                 >  2:1307112/1‑149 (MQ=255)
                                    tCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGTGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTc                                                                                                                 <  2:944107/149‑1 (MQ=255)
                                      aaCGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGa                                                                                                               <  1:1205780/149‑1 (MQ=255)
                                        cGAAAAAGCGGTCGAGATCGGCTGGCGTCCGGATACGGCTGGTGGTACACAGACCGGTAAAAAAGTGCCGATTATCGGCGCAGGCCCGGCAGGACTGGCGTGTGCGGAGGTCCTGCCGCGTAACGGGGTAAAAGCCGTTGTCTTCGTcc                                                                                                             <  1:818090/149‑1 (MQ=255)
                                         gATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGGTGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCg                                                                                                            >  2:464531/1‑149 (MQ=255)
                                         gATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGGTGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCg                                                                                                            >  2:470518/1‑149 (MQ=255)
                                             aaGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTATAAAAGTGGGGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTGCTGACGCGTATCGGCGTAAAAGCCGTTGTCTTCGACCGTCAt                                                                                                        >  1:240176/1‑149 (MQ=255)
                                              aGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGATGTCTTCGACCGTCATc                                                                                                       >  2:923658/1‑149 (MQ=255)
                                               gCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATcc                                                                                                      >  1:1283179/1‑149 (MQ=255)
                                                 gTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGCCCGTCATCCAg                                                                                                    >  1:465156/1‑149 (MQ=255)
                                                 gTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTAGTCTTCGACCGTCATCCAg                                                                                                    >  1:1269666/1‑149 (MQ=255)
                                                         tGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGc                                                                                            <  2:1169454/149‑1 (MQ=255)
                                                           ggCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGCCCGTCATCCAGAAATTGGCgg                                                                                          >  2:37846/1‑149 (MQ=255)
                                                               ggCGTCCGGATATGTCTGGTGTGAAACAGACCTGTAAAAAAGTGGCGATTATCGGCGCAGGCCCTGCAGGTCTTGCGTTTGCGTATGTCCTGACGCGTAACGGCGTAAAAGCCTTTGTCTTCGACCGTCATCCAGAAATTGGCGGgctg                                                                                      <  1:951423/149‑1 (MQ=255)
                                                                     cGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGAcc                                                                                <  2:563314/149‑1 (MQ=255)
                                                                         tatGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCg                                                                            <  2:182662/149‑1 (MQ=255)
                                                                         tatGTCTGGTGGGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCGGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCg                                                                            <  1:188981/149‑1 (MQ=255)
                                                                          atGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCgg                                                                           <  1:285065/149‑1 (MQ=255)
                                                                                    tgAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTTGCGGGCTGCTGACCTTCGGTATTCCGGcc                                                                 >  2:644274/1‑149 (MQ=255)
                                                                                      aaaCAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGAAACGGCGTAAAAGCCGTTGTCTTCGCCCGACATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGccgt                                                               >  1:1523417/1‑147 (MQ=255)
                                                                                      aaaCAGACCGGTAAAAAAGTGGCGATTATCGCCGCCGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCATAAATTGGCGGGCTGCTGACCTTCGGTATTCAGGCCtt                                                               >  2:711227/1‑149 (MQ=255)
                                                                                           gACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGc                                                          >  2:1307488/1‑149 (MQ=255)
                                                                                            aCCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGGTGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCt                                                         <  1:470518/149‑1 (MQ=255)
                                                                                            aCCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGGTGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCt                                                         <  1:464531/149‑1 (MQ=255)
                                                                                            aCCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCt                                                         <  2:729763/149‑1 (MQ=255)
                                                                                            aCCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCt                                                         >  2:413484/1‑149 (MQ=255)
                                                                                             ccGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTg                                                        <  1:117158/149‑1 (MQ=255)
                                                                                                   aaaaaGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGaaaaa                                                  <  2:1274098/149‑1 (MQ=255)
                                                                                                   aaaaaGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGaaaaa                                                  <  2:508439/149‑1 (MQ=255)
                                                                                                     acaGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCATTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAaga                                                <  1:788232/147‑1 (MQ=255)
                                                                                                     aaaGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAaga                                                <  2:2790/149‑1 (MQ=255)
                                                                                                     aaaGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAaga                                                <  2:1181614/149‑1 (MQ=255)
                                                                                                      aaGTGGCGAGTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAagag                                               <  1:786054/149‑1 (MQ=255)
                                                                                                        gTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCTTTTGCGGATGTCCTGACGCGTAACGGCGTAAATGCCGTTGTGTTCGACCGTCATCCAGAAATTCGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGt                                             <  2:57264/149‑1 (MQ=255)
                                                                                                           gCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGTGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAAt                                          >  1:1542193/1‑149 (MQ=255)
                                                                                                             gATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGCCCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGa                                        >  1:319379/1‑149 (MQ=255)
                                                                                                             gATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGa                                        >  1:647240/1‑149 (MQ=255)
                                                                                                              aTTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGAc                                       <  2:773981/149‑1 (MQ=255)
                                                                                                                         aGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTg                            >  2:728221/1‑149 (MQ=255)
                                                                                                                           gCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGaa                          <  1:518042/149‑1 (MQ=255)
                                                                                                                              cGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATc                       <  1:345831/149‑1 (MQ=255)
                                                                                                                                gCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGCCGCGTCGCCGTGAAATCtt                     >  1:184357/1‑149 (MQ=255)
                                                                                                                                  aGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCa                   <  2:513831/149‑1 (MQ=255)
                                                                                                                                  aGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCa                   <  2:1246385/149‑1 (MQ=255)
                                                                                                                                       tGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCGCCGGc              >  1:1233401/1‑149 (MQ=255)
                                                                                                                                                   aTGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGaa  <  1:1368824/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGGCGCGGTGACCATCGGCAACATTGAGCGCTATATCAACGATAAAGCGTTCGAGATGGGCTGGCGTCCGGATATGTCTGGTGTGAAACAGACCGGTAAAAAAGTGGCGATTATCGGCGCAGGCCCGGCAGGTCTGGCGTGTGCGGATGTCCTGACGCGTAACGGCGTAAAAGCCGTTGTCTTCGACCGTCATCCAGAAATTGGCGGGCTGCTGACCTTCGGTATTCCGGCCTTCAAGCTGGAAAAAGAGGTAATGACGCGTCGCCGTGAAATCTTCACCGGCATGGGTATTGAA  >  NC_000913/3359538‑3359833

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