Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,676,209 G→T 5.5% intergenic (+101/+81) yhjE → / ← yhjG putative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,676,2090GT 5.5% 261.0 / 5.1 73intergenic (+101/+81)yhjE/yhjGputative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base T (2/2);  total (37/36)
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.79e-01

TAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGT  >  NC_000913/3676062‑3676356
                                                                                                                                                   |                                                                                                                                                   
tAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                                    <  4:1886629/149‑1 (MQ=255)
                gACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaaa                                                                                                                                    <  3:1617747/149‑1 (MQ=255)
                    cATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcac                                                                                                                                >  4:1088202/1‑149 (MQ=255)
                          aCGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTg                                                                                                                          >  3:1634194/1‑149 (MQ=255)
                              gACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaa                                                                                                                      <  4:1634194/149‑1 (MQ=255)
                                cATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcaca                                                                                                                               <  4:1662334/138‑1 (MQ=255)
                                cATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcaca                                                                                                                               >  3:1662334/1‑138 (MQ=255)
                                     gTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGcc                                                                                                                        >  3:1137850/1‑140 (MQ=255)
                                     gTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGcc                                                                                                                        <  4:1137850/140‑1 (MQ=255)
                                       cGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTaaa                                                                                                             <  3:1088202/149‑1 (MQ=255)
                                         ttGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             >  1:1704240/1‑115 (MQ=37)
                                         ttGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             <  2:1704240/115‑1 (MQ=37)
                                          tGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTg                                                                                                                          <  3:689929/133‑1 (MQ=255)
                                          tGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTg                                                                                                                          >  4:689929/1‑133 (MQ=255)
                                              ataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGc                                                                                                                         <  4:1645351/130‑1 (MQ=255)
                                              ataaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGc                                                                                                                         >  3:1645351/1‑130 (MQ=255)
                                               taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTgtt                                                                                                                        <  4:453659/130‑3 (MQ=255)
                                               taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTgtt                                                                                                                        >  3:453659/1‑128 (MQ=255)
                                               taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaacaa                                                                                                                   >  4:902630/1‑135 (MQ=255)
                                               taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaacaa                                                                                                                   <  3:902630/135‑1 (MQ=255)
                                               taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATc                                                                                                                 <  3:1783171/137‑1 (MQ=255)
                                               taaGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATc                                                                                                                 >  4:1783171/1‑137 (MQ=255)
                                                  gagTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACtt                                                                                                                            <  2:1225462/123‑1 (MQ=255)
                                                  gagTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACtt                                                                                                                            >  1:1225462/1‑123 (MQ=255)
                                                  gagTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaac                                                                                                                     >  2:1310748/1‑130 (MQ=255)
                                                  gagTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaac                                                                                                                     <  1:1310748/130‑1 (MQ=255)
                                                        aaaGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATc                                                                                                           >  2:366174/1‑134 (MQ=255)
                                                        aaaGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATc                                                                                                           <  1:366174/134‑1 (MQ=255)
                                                            ccTGATGCGCTACGCTTCTCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                                                                                                   >  2:627751/1‑106 (MQ=255)
                                                            ccTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                                                                                                   <  1:627751/106‑1 (MQ=255)
                                                                    gcTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcaca                                                                                                                               >  3:312198/1‑102 (MQ=255)
                                                                    gcTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGcaca                                                                                                                               <  4:312198/102‑1 (MQ=255)
                                                                      tACGCTTATCAGGCCTGAATTTCCATACACTATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                                                       <  4:146230/140‑1 (MQ=255)
                                                                      tACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                                                       >  3:146230/1‑140 (MQ=255)
                                                                           ttATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                                       >  4:1313120/1‑119 (MQ=255)
                                                                           ttATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                                       <  3:1313120/119‑1 (MQ=255)
                                                                             aTCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaa                                                                                                                      >  4:1474226/1‑102 (MQ=255)
                                                                             aTCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaa                                                                                                                      <  3:1474226/102‑1 (MQ=255)
                                                                                         tttCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAt                                                                                                   >  2:697168/1‑109 (MQ=255)
                                                                                         tttCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAAt                                                                                                   <  1:697168/109‑1 (MQ=255)
                                                                                              aTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGc                                                                                                <  4:1399488/107‑1 (MQ=255)
                                                                                              aTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGc                                                                                                >  3:1399488/1‑107 (MQ=255)
                                                                                                        ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATc                                                                                 <  1:815198/112‑1 (MQ=255)
                                                                                                        ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATc                                                                                 >  2:815198/1‑112 (MQ=255)
                                                                                                        ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTAtcatca                                                                  <  1:229201/127‑1 (MQ=255)
                                                                                                        ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTAtcatca                                                                  >  2:229201/1‑127 (MQ=255)
                                                                                                                gCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                               >  4:415569/1‑90 (MQ=255)
                                                                                                                gCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                               <  3:415569/90‑1 (MQ=255)
                                                                                                                   aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                                       <  1:900248/79‑1 (MQ=255)
                                                                                                                   aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGtt                                                                                                       >  2:900248/1‑79 (MQ=255)
                                                                                                                      aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                               >  2:305317/1‑84 (MQ=255)
                                                                                                                      aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                               >  2:1688587/1‑84 (MQ=255)
                                                                                                                      aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                               <  1:305317/84‑1 (MQ=255)
                                                                                                                      aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                                                               <  1:1688587/84‑1 (MQ=255)
                                                                                                                      aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCg                                                                               <  2:671769/100‑1 (MQ=255)
                                                                                                                      aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCg                                                                               >  1:671769/1‑100 (MQ=255)
                                                                                                                       tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTa                             >  3:754277/1‑149 (MQ=255)
                                                                                                                                cGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATcc                                                                                <  3:261690/89‑1 (MQ=25)
                                                                                                                                cGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCt                    >  1:1121424/1‑149 (MQ=255)
                                                                                                                                       tttgTTCATGCCGGATGCGGCGTGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTc                  <  3:109689/141‑1 (MQ=255)
                                                                                                                                       tttgTTCATGCCGGATGCGGCGTGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTc                  >  4:109689/4‑144 (MQ=255)
                                                                                                                                       ggCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc                                         <  2:1396336/121‑1 (MQ=255)
                                                                                                                                       ggCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc                                         >  1:1396336/1‑121 (MQ=255)
                                                                                                                                       ggCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAg             >  3:204260/1‑149 (MQ=255)
                                                                                                                                        gCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                                                          <  2:1188469/87‑1 (MQ=255)
                                                                                                                                        gCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                                                          >  1:1188469/1‑87 (MQ=255)
                                                                                                                                        gCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg                 <  2:103633/144‑1 (MQ=255)
                                                                                                                                        gCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg                 >  1:103633/1‑144 (MQ=255)
                                                                                                                                          gTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTc                        <  3:457584/135‑1 (MQ=255)
                                                                                                                                          gTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTc                        >  4:457584/1‑135 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGt                                                                                                        <  2:1874739/54‑1 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGt                                                                                                        >  1:1874739/1‑54 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTga         >  4:16774/1‑149 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTga         <  4:204260/149‑1 (MQ=255)
                                                                                                                                                  cgcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGt  <  4:754277/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAATCGCCCTGCTGCTGACCCATGAGACGCGACATCAGTCGTTGTAATAAGAGTGAAAAGCCTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGT  >  NC_000913/3676062‑3676356

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