Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,152 G→A 10.1% intergenic (+40/‑168) ihfB → / → ycaI integration host factor subunit beta/conserved inner membrane protein YcaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,1520GA10.1% 169.5 / 7.7 60intergenic (+40/‑168)ihfB/ycaIintegration host factor subunit beta/conserved inner membrane protein YcaI
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (2/4);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACT  >  NC_000913/964004‑964297
                                                                                                                                                    |                                                                                                                                                 
tACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                                                                                                  >  1:58039/1‑150 (MQ=255)
       cGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAGAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttt                                                                                                                                           <  1:950097/150‑1 (MQ=255)
         tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                                                                         <  1:783169/150‑2 (MQ=255)
             ccGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAt                                                                                                                                     >  1:684710/1‑150 (MQ=255)
              cGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                                                                    <  1:1117246/150‑1 (MQ=255)
               gAAGACTGGCGTTAAAGTAGAACTGGAAGGAAACTACGTTCCTCTCTTTAAACCTGGTAAATACGTGCGCGATCGCTCCAATATTTACGGTTAAGTTTTTTACTCATACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTTGTctcgt                                                                                                                                   <  2:367965/150‑5 (MQ=255)
               gAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                                                                   <  2:476786/150‑1 (MQ=255)
               gAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                                                                   >  2:590898/1‑150 (MQ=255)
                      ggCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGAGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACTAGAGAAAAGAACCTGTCGGGTGCTTTTTTCATTTCTCTaaa                                                                                                                            >  2:841829/1‑149 (MQ=255)
                        cGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                                                                          <  2:1077540/150‑1 (MQ=255)
                           tAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCTGGTGCTTTTTTCATTTCTCTAATctgga                                                                                                                       <  2:490325/150‑1 (MQ=255)
                                 aGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgga                                                                                                                 <  1:504021/150‑1 (MQ=255)
                                     cTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGACCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAGCTGTCGGGTGCTGTTTTCATTTCTGTAATCTGGAACTGGAAGCt                                                                                                             >  2:1348314/1‑150 (MQ=255)
                                     cTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCt                                                                                                             >  1:551592/1‑150 (MQ=255)
                                     cTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCt                                                                                                             <  2:853420/150‑1 (MQ=255)
                                      tGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTg                                                                                                            >  1:1286079/1‑150 (MQ=255)
                                          aggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTc                                                                                                        >  2:1308322/1‑150 (MQ=255)
                                                aTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCagag                                                                                                  <  1:1007268/150‑1 (MQ=255)
                                                 tACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGt                                                                                                 <  2:92554/150‑1 (MQ=255)
                                                    gTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                              >  2:765762/1‑150 (MQ=255)
                                                    gTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                              <  2:1342292/150‑1 (MQ=255)
                                                        cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                          >  1:141499/1‑150 (MQ=255)
                                                         tCACTTTAAACCTGGTAAAGACCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAACAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACa                                                                                         <  2:479695/150‑1 (MQ=255)
                                                         tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCt                                                                                                         <  1:488815/134‑1 (MQ=255)
                                                         tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCt                                                                                                         >  2:488815/1‑134 (MQ=255)
                                                         tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACa                                                                                         >  2:65483/1‑150 (MQ=255)
                                                                    cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCTGAGTTTTGAACAGTTTTCACCCt                                                                              <  2:1179577/150‑1 (MQ=255)
                                                                       gTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGATCTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                                                                           <  2:1362841/150‑1 (MQ=255)
                                                                        tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg                                                                          <  1:1308322/150‑1 (MQ=255)
                                                                              aCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAt                                                                    >  2:1335757/1‑150 (MQ=255)
                                                                                   gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc                                                               >  2:641914/1‑150 (MQ=255)
                                                                                   gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc                                                               <  1:666526/150‑1 (MQ=255)
                                                                                   gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc                                                               <  1:600377/150‑1 (MQ=255)
                                                                                    cgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct                                                              >  2:247324/1‑150 (MQ=255)
                                                                                     gATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg                                                             <  2:551592/150‑1 (MQ=255)
                                                                                     gATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg                                                             <  2:684710/150‑1 (MQ=255)
                                                                                     gATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCAGTCGGGAGCTTTTTTCATTTCTCTAATCTTGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg                                                             >  2:929038/1‑150 (MQ=255)
                                                                                          cgcCANTATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAatat                                                        >  1:546241/1‑150 (MQ=255)
                                                                                           gcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATg                                                       <  1:765762/150‑1 (MQ=255)
                                                                                           gcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATg                                                       <  1:377470/150‑1 (MQ=255)
                                                                                             cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGcc                                                     <  2:1160205/150‑1 (MQ=255)
                                                                                                   ttACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCgggg                                               <  2:58039/150‑1 (MQ=255)
                                                                                                       ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTTATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACg                                           >  2:391746/1‑150 (MQ=255)
                                                                                                          tAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCaa                                        >  2:25117/1‑150 (MQ=255)
                                                                                                             gTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAAt                                     <  1:247324/150‑1 (MQ=255)
                                                                                                                   tACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccac                               <  1:65483/150‑1 (MQ=255)
                                                                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGTGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCTGGGAACGCAAAATTCccacc                              <  2:1091689/150‑1 (MQ=255)
                                                                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                              <  1:590898/150‑1 (MQ=255)
                                                                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                              >  2:931378/1‑150 (MQ=255)
                                                                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                           >  1:1381069/1‑150 (MQ=255)
                                                                                                                                   tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCt                   <  1:40642/142‑1 (MQ=255)
                                                                                                                                   tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCt                   >  2:40642/5‑146 (MQ=255)
                                                                                                                                   tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac               >  1:308332/5‑150 (MQ=255)
                                                                                                                                   agagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATATGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACAGCTcaaa               >  2:1008154/1‑149 (MQ=255)
                                                                                                                                       aaaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAg           >  2:715389/1‑150 (MQ=255)
                                                                                                                                        aaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTcc                                  <  1:978587/126‑1 (MQ=255)
                                                                                                                                        aaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTcc                                  >  2:978587/1‑126 (MQ=255)
                                                                                                                                        aaaGCACCTGTCGCGTGCTTTTTTCATTTCTCTAATCTCGAACTGCAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACACCCGCTCAACAAAGt          >  2:481109/1‑150 (MQ=255)
                                                                                                                                          aGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTg        <  1:560943/150‑1 (MQ=255)
                                                                                                                                               cTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCacac   <  2:546241/150‑1 (MQ=255)
                                                                                                                                                tGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACt  >  1:340359/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACT  >  NC_000913/964004‑964297

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