Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 851,934 Δ1 bp 5.1% coding (664/1584 nt) opgE ← OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913851,9340C. 5.1% 268.2 / 7.8 78coding (664/1584 nt)opgEOPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base . (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGATAGC  >  NC_000913/851788‑852080
                                                                                                                                                  |                                                                                                                                                  
cGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATcacca                                                                                                                                                  <  4:178557/149‑1 (MQ=255)
cGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATcacccc                                                                                                                                                 >  2:266383/1‑147 (MQ=255)
     aGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACcaca                                                                                                                                             >  4:172985/1‑149 (MQ=255)
        cAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAg                                                                                                                                          <  4:285204/149‑1 (MQ=255)
               tttCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGa                                                                                                                                    >  4:9951/1‑148 (MQ=255)
                 tCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATa                                                                                                                                  <  6:133652/148‑1 (MQ=255)
                   tGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTc                                                                                                                                       >  3:337061/1‑141 (MQ=255)
                   tGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTc                                                                                                                                       <  4:337061/141‑1 (MQ=255)
                   tGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACtt                                                                                                                               >  6:58627/1‑149 (MQ=255)
                   tGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACtt                                                                                                                               >  6:117976/1‑149 (MQ=255)
                       tGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGc                                                                                                                           <  5:173463/149‑1 (MQ=255)
                            ccACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAAc                                                                                                                      >  6:207943/1‑149 (MQ=255)
                            ccACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAAc                                                                                                                      <  6:233512/149‑1 (MQ=255)
                                tGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACgg                                                                                                                    >  1:186255/1‑147 (MQ=255)
                                  gcggcAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCCGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCGTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACggacgt                                                                                                                >  6:341376/1‑145 (MQ=255)
                                  gcggcAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCa                                                                                                                  <  5:207943/147‑1 (MQ=255)
                                    ggcAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGa                                                                                                              >  2:9177/1‑149 (MQ=255)
                                       aaTGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGCCCGCCAAAGGCACCGCGGGCATCGCGCCGCGCGCTGTCCCCGATCCCCCCCACATAGGTCTGAGACTTCGGCTTACCGGCAGTGCCCg                                                                                                           >  1:52063/1‑149 (MQ=255)
                                         tGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTc                                                                                                         >  2:386684/1‑149 (MQ=255)
                                            aGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATTGGTCTg                                                                                                                                     >  5:141393/1‑118 (MQ=255)
                                            aGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATTGGTCTg                                                                                                                                     <  6:141393/118‑1 (MQ=255)
                                             gTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAg                                                                                                      <  1:116070/148‑1 (MQ=255)
                                             gTCAGCAAATATCAAACCCTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCCCCGAGGCCATCGCGACGCGAGCTTTCCCCGATCACCACCACATAGGTCTGCGACTTCGGCTTAACGGCAGTGACCGTCCAgg                                                                                                     >  1:104677/1‑149 (MQ=255)
                                              tCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAgg                                                                                                     >  7:143708/1‑148 (MQ=255)
                                                 gCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATc                                                                                                 >  6:238961/1‑149 (MQ=255)
                                                  cAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCTGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCt                                                                                                >  8:217686/1‑149 (MQ=255)
                                                    aaTATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTcc                                                                                                        >  8:343683/1‑139 (MQ=255)
                                                    aaTATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTcc                                                                                                        <  7:343683/139‑1 (MQ=255)
                                                         cAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTc                                                                                         >  2:422023/1‑149 (MQ=255)
                                                                ttAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAg                                                                                  <  4:86946/149‑1 (MQ=255)
                                                                   aCGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGtt                                                                                <  5:294792/148‑1 (MQ=255)
                                                                    cGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGtttt                                                                              >  1:369472/1‑149 (MQ=255)
                                                                    cGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGtttt                                                                              <  1:422023/149‑1 (MQ=255)
                                                                     gctgctGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCCCCGATCACCACCACATAGGTCTGAAACATCGGCTTAACGGCAGTGACCGACCAGGAATCTTTCATTCCGCAAAGTTATg                                                                             >  1:390301/1‑149 (MQ=255)
                                                                          tgGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATc                                                                        >  2:65564/1‑149 (MQ=255)
                                                                            gCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCg                                                                      >  3:47301/1‑149 (MQ=255)
                                                                              aaaCGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTc                                                                   >  6:47166/1‑149 (MQ=255)
                                                                               aaCGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTc                                                                   <  3:9951/149‑1 (MQ=255)
                                                                                aCGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCt                                                                  <  5:58627/149‑1 (MQ=255)
                                                                                 cGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTg                                                                <  5:47166/149‑1 (MQ=255)
                                                                                   gggTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGc                                                               <  5:117976/149‑1 (MQ=255)
                                                                                     gTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGctc                                                             <  3:172985/149‑1 (MQ=255)
                                                                                      tATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCGGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAACGTTTTGCCATCCGTTCTTGctcg                                                            >  2:286807/1‑148 (MQ=255)
                                                                                        ttGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTc                                                          <  2:113573/149‑1 (MQ=255)
                                                                                           tCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAAt                                                       >  3:80913/1‑149 (MQ=255)
                                                                                           tCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGTTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCtt                                                           <  8:147272/144‑1 (MQ=255)
                                                                                           tCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCtt                                                           >  7:147272/1‑144 (MQ=255)
                                                                                             ccAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCt                                                                                                <  3:266280/106‑1 (MQ=255)
                                                                                             ccAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCt                                                                                                >  4:266280/1‑106 (MQ=255)
                                                                                             ccAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCa                                                     >  7:113428/1‑149 (MQ=255)
                                                                                               aGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCAcc                                                   >  7:47799/1‑149 (MQ=255)
                                                                                                   aCCGCCAACGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCg                                                <  5:154213/148‑1 (MQ=255)
                                                                                                   aCCGCCAAAGGCACCGAGGGCTGCGCGACGCGCGCTTTCACCGATCACCACCACATCGGTGTGATACTTCGGCTTATCGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCg                                                <  2:390301/148‑1 (MQ=255)
                                                                                                     cGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGtta                                             <  7:230429/149‑1 (MQ=255)
                                                                                                     cGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGtta                                             <  1:386684/149‑1 (MQ=255)
                                                                                                           aGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTATTATTGa                                      <  5:369015/149‑1 (MQ=255)
                                                                                                             gCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGaaa                                                                                   >  2:189598/1‑103 (MQ=255)
                                                                                                             gCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGaaa                                                                                   <  1:189598/103‑1 (MQ=255)
                                                                                                               aCCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGa                                   <  3:333016/149‑1 (MQ=255)
                                                                                                                ccGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGAt                                  >  7:25069/1‑149 (MQ=255)
                                                                                                                     ggCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACAt                             <  6:271201/149‑1 (MQ=255)
                                                                                                                      gCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATc                            <  7:217686/149‑1 (MQ=255)
                                                                                                                       cATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCg                                                                                                                             <  5:146900/51‑1 (MQ=255)
                                                                                                                       cATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCg                                                                                                                             >  6:146900/1‑51 (MQ=255)
                                                                                                                        aTCGCGACGCGCGCTTTCACCGATCACCCCCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCCCCTCGTTATTATTGACGATGACATCCt                          >  5:179665/1‑149 (MQ=255)
                                                                                                                         tCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCtt                         <  8:156965/149‑1 (MQ=255)
                                                                                                                          cgcgACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTTTTGACGa                                   <  3:81559/138‑1 (MQ=255)
                                                                                                                          cgcgACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGa                                   >  4:81559/1‑138 (MQ=255)
                                                                                                                           gcgACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAg                       <  2:186255/149‑1 (MQ=255)
                                                                                                                              aCGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCa                     >  5:7319/1‑148 (MQ=255)
                                                                                                                                gcgcgcTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGc                      <  2:279482/145‑1 (MQ=255)
                                                                                                                                gcgcgcTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGc                      >  1:279482/1‑145 (MQ=255)
                                                                                                                                  gcgcTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTTCTCTTCAATCACCTCGTTAATATTGACGATGACATCCTTAGCAAAAc                 >  6:271209/1‑148 (MQ=255)
                                                                                                                                  gcgcTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACg                <  5:238961/149‑1 (MQ=255)
                                                                                                                                     cTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGaaa             <  4:80913/149‑1 (MQ=255)
                                                                                                                                     cTTTCACCGATCACCACCACATACGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGaaa             <  2:362222/149‑1 (MQ=255)
                                                                                                                                          aCCGATCCACACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCAGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCg        <  2:52063/149‑1 (MQ=255)
                                                                                                                                          aCCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCg        >  2:127144/1‑149 (MQ=255)
                                                                                                                                          aCCGATCACCACCACATAGGTCAGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCg        >  8:253615/1‑149 (MQ=255)
                                                                                                                                             gATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTACAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGa      >  8:13085/1‑148 (MQ=255)
                                                                                                                                                 accaccacATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGATAGc  <  2:369472/148‑1 (MQ=255)
                                                                                                                                                  ccaccacATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGAt     >  8:144337/1‑144 (MQ=255)
                                                                                                                                                  ccaccacATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGAt     <  7:144337/144‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGGTATTGTCCCAGTGACCGCCAAAGGCACCGAGGGCATCGCGACGCGCGCTTTCACCGATCACCACCACATAGGTCTGATACTTCGGCTTAACGGCAGTGACCGTCCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGATAGC  >  NC_000913/851788‑852080

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