Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,995 T→G 11.2% intergenic (+16/+23) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9950TG11.2% 277.5 / 28.6 89intergenic (+16/+23)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (40/39);  new base G (5/5);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

ACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATG  >  NC_000913/2364858‑2365143
                                                                                                                                         |                                                                                                                                                    
aCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTc                                                                                                                                           >  2:733985/1‑149 (MQ=255)
               tAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                            <  1:733985/149‑1 (MQ=255)
                 aaCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                                                                          <  4:802916/149‑1 (MQ=255)
                    gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                                                                       >  3:108549/1‑149 (MQ=255)
                    gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                                                                       >  3:147422/1‑149 (MQ=255)
                    gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                                                                       >  3:39083/1‑149 (MQ=255)
                    gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                                                                       >  1:218867/1‑149 (MQ=255)
                    gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTc                                                                                                                                           >  2:457892/1‑129 (MQ=255)
                    gAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTc                                                                                                                                           <  1:457892/129‑1 (MQ=255)
                          tttCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                 <  4:147422/149‑1 (MQ=255)
                          tttCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                 <  4:39083/149‑1 (MQ=255)
                            tCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                                                                    <  2:1164287/144‑1 (MQ=255)
                            tCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                                                                    >  1:1164287/1‑144 (MQ=255)
                                gACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                 >  2:877658/1‑127 (MQ=255)
                                gACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                 <  3:1320133/127‑1 (MQ=255)
                                gACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                 >  4:1320133/1‑127 (MQ=255)
                                gACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                 <  1:877658/127‑1 (MQ=255)
                                   tACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGAGTTGAATGTCACCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTAGTTAACGAATTCCCCGTc                                                                                                        <  4:108549/149‑1 (MQ=255)
                                            gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               >  4:1213625/1‑149 (MQ=255)
                                             tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCAGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAtt                                                                                                                               <  3:394524/116‑1 (MQ=255)
                                             tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCAGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAtt                                                                                                                               >  4:394524/1‑116 (MQ=255)
                                               aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACttt                                                                                            <  2:218867/149‑1 (MQ=255)
                                                aaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAATGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                            <  1:646817/116‑1 (MQ=255)
                                                aaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                                                                            >  2:646817/1‑116 (MQ=255)
                                                   cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                                                                        >  3:686444/1‑149 (MQ=255)
                                                     gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                                                                      >  2:307063/1‑133 (MQ=255)
                                                     gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                                                                      <  1:307063/133‑1 (MQ=255)
                                                        gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                                                                             >  2:155438/1‑123 (MQ=255)
                                                        gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                                                                             <  1:155438/123‑1 (MQ=255)
                                                        gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                                                                   <  1:738985/149‑1 (MQ=255)
                                                        gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                                                                   >  2:738985/1‑149 (MQ=255)
                                                            tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                                                                               >  4:725857/1‑149 (MQ=255)
                                                              gTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAGCGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatg                                                                             >  1:384143/1‑149 (MQ=255)
                                                                       gATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcat                                                                          >  3:1329727/1‑143 (MQ=255)
                                                                       gATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcat                                                                          <  4:1329727/143‑1 (MQ=255)
                                                                        aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                                                                    >  1:1183499/1‑100 (MQ=255)
                                                                        aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                                                                    <  2:1183499/100‑1 (MQ=255)
                                                                                    ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCTCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa                                                            >  1:876052/1‑144 (MQ=255)
                                                                                    ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCTCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa                                                            <  2:876052/144‑1 (MQ=255)
                                                                                      aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCt                                                                                                                                              >  3:1424803/1‑60 (MQ=255)
                                                                                      aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCt                                                                                                                                              <  4:1424803/60‑1 (MQ=255)
                                                                                      aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAGCGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATc                                                     <  2:384143/149‑1 (MQ=255)
                                                                                      aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATc                                                     <  3:1213625/149‑1 (MQ=255)
                                                                                       cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa                                                    <  4:686444/149‑1 (MQ=255)
                                                                                         cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                                                                   <  4:84594/132‑1 (MQ=255)
                                                                                         cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                                                                   >  3:84594/1‑132 (MQ=255)
                                                                                          gAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt                                                 <  3:725857/149‑1 (MQ=255)
                                                                                                 gttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa                                                       <  2:818860/136‑1 (MQ=255)
                                                                                                 gttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa                                                       >  1:818860/1‑136 (MQ=255)
                                                                                                  ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                 >  1:484514/1‑77 (MQ=255)
                                                                                                  ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                 <  2:484514/77‑1 (MQ=255)
                                                                                                            aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                               >  4:635207/1‑149 (MQ=255)
                                                                                                             aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcat                                                                          >  1:1209319/1‑105 (MQ=255)
                                                                                                             aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcat                                                                          <  2:1209319/105‑1 (MQ=255)
                                                                                                               gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcccc                                                                                                           <  3:1191818/70‑1 (MQ=255)
                                                                                                               gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcccc                                                                                                           >  4:1191818/1‑70 (MQ=255)
                                                                                                                  ttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                                                               >  1:1371562/1‑63 (MQ=255)
                                                                                                                  ttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                                                               <  2:1371562/63‑1 (MQ=255)
                                                                                                                     tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               >  3:1056593/1‑76 (MQ=255)
                                                                                                                     tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               <  4:1056593/76‑1 (MQ=255)
                                                                                                                            atgatgACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTAt                           >  3:1052573/1‑137 (MQ=255)
                                                                                                                            atgatgACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTAt                           <  4:1052573/137‑1 (MQ=255)
                                                                                                                             atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                                <  3:1274798/112‑1 (MQ=255)
                                                                                                                             atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                                >  4:1274798/4‑115 (MQ=255)
                                                                                                                              gatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat             >  3:965190/1‑149 (MQ=255)
                                                                                                                               ctgaCAGCCAGTAATGGCTGTCGGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgc                                      <  3:1113081/122‑1 (MQ=255)
                                                                                                                               ctgaCAGCCAGTAATGGCTGTCGGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgc                                      >  4:1113081/2‑123 (MQ=255)
                                                                                                                               ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               <  3:244548/65‑1 (MQ=38)
                                                                                                                               ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               >  4:244548/2‑66 (MQ=38)
                                                                                                                                 gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                               >  1:1299139/1‑128 (MQ=255)
                                                                                                                                 gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                               <  2:1299139/128‑1 (MQ=255)
                                                                                                                                    aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                                                                             <  3:1433508/47‑1 (MQ=255)
                                                                                                                                    aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                                                                             >  4:1433508/1‑47 (MQ=255)
                                                                                                                                    aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAgg                                                  <  2:812925/106‑1 (MQ=255)
                                                                                                                                    aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAgg                                                  >  1:812925/1‑106 (MQ=255)
                                                                                                                                    aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAgg       <  3:635207/149‑1 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt                                                 <  2:769355/105‑1 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt                                                 >  1:769355/1‑105 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgc                                      <  2:960910/116‑1 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgc                                      >  1:960910/1‑116 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                 <  4:456837/137‑1 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                 <  1:632069/137‑1 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                 >  3:456837/1‑137 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                 >  2:632069/1‑137 (MQ=255)
                                                                                                                                      ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCa     >  2:1127022/1‑149 (MQ=255)
                                                                                                                                       cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAg        <  3:1340536/145‑1 (MQ=255)
                                                                                                                                       cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAg        >  4:1340536/1‑145 (MQ=255)
                                                                                                                                        aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGa                                                                        <  1:144365/80‑1 (MQ=255)
                                                                                                                                        aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGa                                                                        >  2:144365/1‑80 (MQ=255)
                                                                                                                                         ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                              <  3:806925/57‑1 (MQ=255)
                                                                                                                                         ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                              >  4:806925/1‑57 (MQ=255)
                                                                                                                                         ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATg  >  2:1202590/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
ACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATG  >  NC_000913/2364858‑2365143

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