Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,995 T→G 6.2% intergenic (+30/+19) ompX → / ← opgE outer membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913850,9950TG 6.2% 260.9 / 7.1 77intergenic (+30/+19)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (30/42);  new base G (2/3);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCT  >  NC_000913/850850‑851141
                                                                                                                                                 |                                                                                                                                                  
cGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCtctc                                                                                                                                                   <  2:449796/147‑1 (MQ=255)
cGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCtctc                                                                                                                                                   >  1:449796/1‑147 (MQ=255)
   gTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                              <  4:529464/149‑1 (MQ=255)
        cAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtt                                                                                                                                         <  4:828891/149‑1 (MQ=255)
         aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttt                                                                                                                                        <  3:187739/149‑1 (MQ=255)
         aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtt                                                                                                                                         <  1:601135/148‑1 (MQ=255)
          gTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                              >  3:390980/1‑142 (MQ=255)
          gTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                              <  4:390980/142‑1 (MQ=255)
           ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATGAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      >  4:579627/1‑149 (MQ=255)
           ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttt                                                                                                                                      >  4:122761/1‑149 (MQ=255)
           ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttt                                                                                                                                      <  3:122761/149‑1 (MQ=255)
             cAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTg                                                                                                                                     <  2:507807/148‑1 (MQ=255)
              aaCCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGt                                                                                                                                    <  1:644078/148‑1 (MQ=255)
                 ccGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtttt                                                                                                                                 >  1:192184/1‑148 (MQ=255)
                 ccGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtttt                                                                                                                                 >  2:774618/1‑148 (MQ=255)
                   gATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                              <  3:678192/149‑1 (MQ=255)
                         aaaCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTc                                                                                                                        <  3:579627/149‑1 (MQ=255)
                          aaCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTc                                                                                                                        <  2:116722/148‑1 (MQ=255)
                           aCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTa                                                                                                                               >  1:646641/1‑140 (MQ=255)
                           aCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTa                                                                                                                               <  2:646641/140‑1 (MQ=255)
                             gTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCggg                                                                                                                     <  1:103095/148‑1 (MQ=255)
                                  tctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaa                                                                                                               >  4:111385/1‑149 (MQ=255)
                                   ctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaa                                                                                                               >  2:361484/1‑148 (MQ=255)
                                   ctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaa                                                                                                               >  2:175320/1‑148 (MQ=255)
                                     ggACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAata                                                                                                            <  4:583455/149‑1 (MQ=255)
                                     ggACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAata                                                                                                            >  4:864057/1‑149 (MQ=255)
                                     ggACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAat                                                                                                             <  1:548159/148‑1 (MQ=255)
                                        cTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                                                                          <  1:175320/148‑1 (MQ=255)
                                         ttCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  4:887018/119‑1 (MQ=255)
                                         ttCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      >  3:887018/1‑119 (MQ=255)
                                         ttCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCgg                                                                                                        >  3:97747/1‑149 (MQ=255)
                                            tctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTc                                                                                                     >  4:278436/1‑149 (MQ=255)
                                             ctTACGAGCAGAGCCGTACTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCCCTTTGGTGATATAAAAAATCCGCCGCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                                    <  4:608504/149‑1 (MQ=255)
                                               tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                                                                   >  2:428064/1‑148 (MQ=255)
                                                   aGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGAt                                                                                               >  2:108914/1‑148 (MQ=255)
                                                     cAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCGCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATgg                                                                                             >  2:172729/1‑148 (MQ=255)
                                                      agagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                           <  4:539714/149‑1 (MQ=255)
                                                      agagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                           <  3:728293/149‑1 (MQ=255)
                                                      agagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGt                                                                                            <  2:192184/148‑1 (MQ=255)
                                                        agCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtttt                                                                                                                                       <  1:611358/103‑1 (MQ=255)
                                                        agCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGt                                                                                         >  3:749429/1‑149 (MQ=255)
                                                        agCCGTATTCGTAACGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtttt                                                                                                                                       >  2:611358/1‑103 (MQ=255)
                                                            gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                           >  3:138106/1‑143 (MQ=255)
                                                            gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                           <  4:138106/143‑1 (MQ=255)
                                                            gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATc                                                                                     >  3:141551/1‑149 (MQ=255)
                                                              aTTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                                                                          <  1:158165/126‑1 (MQ=255)
                                                              aTTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                                                                          >  2:158165/1‑126 (MQ=255)
                                                                     gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCt                                                                             <  1:428064/148‑1 (MQ=255)
                                                                         tGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTaa                                                                         >  2:409725/1‑148 (MQ=255)
                                                                              tAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGa                                                                    <  1:409725/148‑1 (MQ=255)
                                                                                  cACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg                                                               <  4:637919/149‑1 (MQ=255)
                                                                                   aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAc                                                                                  <  3:650693/129‑1 (MQ=255)
                                                                                   aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGATTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAc                                                                                  >  4:650693/1‑129 (MQ=255)
                                                                                    ccTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAAATGGAAGTTgg                                                              <  1:677108/148‑1 (MQ=255)
                                                                                      tGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGa                                                           <  3:864057/149‑1 (MQ=255)
                                                                                         aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt                                                        <  3:837273/149‑1 (MQ=255)
                                                                                              cGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTatat                                                    <  1:698643/148‑1 (MQ=255)
                                                                                                 tgtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGc                                                <  4:141551/149‑1 (MQ=255)
                                                                                                        tACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCtt                                         <  3:111385/149‑1 (MQ=255)
                                                                                                          ccGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                       <  3:278436/149‑1 (MQ=255)
                                                                                                              ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                       <  3:354724/145‑1 (MQ=255)
                                                                                                              ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                       >  4:354724/1‑145 (MQ=255)
                                                                                                                  aaTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAg                                <  1:886724/148‑1 (MQ=255)
                                                                                                                   aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGa                                                           >  2:640383/1‑120 (MQ=255)
                                                                                                                   aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGa                                                           <  1:640383/120‑1 (MQ=255)
                                                                                                                   aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGat                              >  3:717597/1‑149 (MQ=255)
                                                                                                                    tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                       >  4:904902/1‑139 (MQ=255)
                                                                                                                    tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                       <  3:904902/139‑1 (MQ=255)
                                                                                                                     cACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGata                             >  2:189915/1‑148 (MQ=255)
                                                                                                                        tttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATaaa                         <  4:717597/149‑1 (MQ=255)
                                                                                                                           ggTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACgg                      <  4:749429/149‑1 (MQ=255)
                                                                                                                           ggTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACg                       <  1:361484/148‑1 (MQ=255)
                                                                                                                              gATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTa                   <  4:756872/149‑1 (MQ=255)
                                                                                                                                 ataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTAtt                 <  1:311755/148‑1 (MQ=255)
                                                                                                                                   aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTg               >  1:758682/1‑148 (MQ=255)
                                                                                                                                   aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTCCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGa              <  4:141689/149‑1 (MQ=255)
                                                                                                                                     aaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGAt             >  1:397513/1‑148 (MQ=255)
                                                                                                                                          ccGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTaa        >  1:434634/1‑148 (MQ=255)
                                                                                                                                             ccTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAttt     >  2:451076/1‑148 (MQ=255)
                                                                                                                                               tctcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCt  >  4:116492/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCT  >  NC_000913/850850‑851141

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: