Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,605,732 G→A 6.2% A228T (GCC→ACC)  lsrB → autoinducer 2‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,605,7320GA 6.2% 200.2 / 2.7 64A228T (GCC→ACC) lsrBautoinducer 2‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (2/2);  total (37/27)
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.12e-01

CAAAGAGCATCCAGGCTGGGAAATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTG  >  NC_000913/1605587‑1605878
                                                                                                                                                 |                                                                                                                                                  
cAAAGAGCATCCAGGCTGGGATATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCCTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACTCCCTGCCCGCTGCCg                                                                                                                                                 >  3:300300/1‑149 (MQ=255)
              gCTGGGAATTTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGa                                                                                                                                    >  8:272655/1‑148 (MQ=255)
                tGGGAAATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTACCGCACAAGCCGCAGaaaa                                                                                                                                 >  8:34555/1‑149 (MQ=255)
                  ggAAATTGTCTCTACGCAGTTTGGCTATAACGATGCCAGTAAATCGTTACAAAACGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAAc                                                                                                                                <  4:300300/148‑1 (MQ=255)
                  ggAAATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGca                                                                                                                                               <  4:271883/133‑1 (MQ=255)
                  ggAAATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGca                                                                                                                                               >  3:271883/1‑133 (MQ=255)
                  ggAAATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAAc                                                                                                                                <  7:86217/148‑1 (MQ=255)
                      aTTGCCACTTCGCACTTTGGCTATAACGATGGCACTAAATCGTTACAAACCGCTGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATTCCAACGCCCTGCCCGCTGCCGCACAAGCCGGAGAAAACTTGaa                                                                                                                           <  3:279838/149‑1 (MQ=255)
                             cTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGa                                                                                                                     >  7:299684/1‑148 (MQ=255)
                              tACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACa                                                                                                                   <  3:312477/149‑1 (MQ=255)
                               aCGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTACCGCACAAGCCGCAGAAAACTTGAAAAATGACa                                                                                                                   >  5:308878/1‑148 (MQ=255)
                                     tttGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGc                                                                                                            >  8:304328/1‑149 (MQ=255)
                                          cTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCAAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCTATTg                                                                                                       >  1:322892/1‑149 (MQ=255)
                                            ataACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCCACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTc                                                                                                     >  5:340082/1‑149 (MQ=255)
                                               aCGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGa                                                                                                  >  8:201456/1‑149 (MQ=255)
                                                 gATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGAtt                                                                                                <  1:311218/149‑1 (MQ=255)
                                                 gATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGAtt                                                                                                <  2:342333/149‑1 (MQ=255)
                                                 gATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGAtt                                                                                                >  4:108981/1‑149 (MQ=255)
                                                 gATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTACCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGAtt                                                                                                <  7:34555/149‑1 (MQ=255)
                                                        cTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACg                                                                                         >  8:56050/1‑149 (MQ=255)
                                                               gTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAAt                                                                                   <  8:194181/148‑1 (MQ=255)
                                                                ttACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCa                                                                                                                                      <  1:388351/96‑1 (MQ=255)
                                                                ttACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCa                                                                                                                                      >  2:388351/1‑96 (MQ=255)
                                                                   cAAACCGCAGCAGGAATATTAACAGCGTATAGCGATCTCGACGCCATAATCGCCCCCGATGCCAACGCCCTGCCCGCGGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGa                                                                               <  6:224076/148‑1 (MQ=255)
                                                                   cAAACCGCAGAAGGAATATTAAAACCGTATAGCGAACTCGACGCCATTATCGCCCCCGAAGCCAACGCCCTCCCCCCTGCCCCACAAGCCGCAGAAAACTTGAAAAATGACAAAGAAGCGATTCTCGGATACAGAACGCCAAATGTCAt                                                                              >  8:78972/1‑149 (MQ=255)
                                                                    aaaCCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATg                                                                             >  8:268638/1‑149 (MQ=255)
                                                                    aaaCCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATg                                                                             >  7:203304/1‑149 (MQ=255)
                                                                           aGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGt                                                                      <  6:337142/149‑1 (MQ=255)
                                                                            gAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTa                                                                     >  2:401589/1‑149 (MQ=255)
                                                                              aGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATg                                                                   <  8:129637/149‑1 (MQ=255)
                                                                              aGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATg                                                                   <  2:254482/149‑1 (MQ=255)
                                                                               ggAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCAAACTCCCTGCCCGCTGCAGCACAAGCCGCAGAAAACTTGAAAAATGACATAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATg                                                                   >  3:114708/1‑148 (MQ=255)
                                                                                   tatTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCCCAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTagag                                                              >  6:218497/1‑149 (MQ=255)
                                                                                    atTAAAAGCGTATAGCGATCTCGCCGCGATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGTTTGTCGGATGCAGTACGCCAATTGTGATGCGCCCGTATGTAGAgc                                                             >  3:88028/1‑149 (MQ=255)
                                                                                       aaaaGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGcgcg                                                          >  5:53150/1‑149 (MQ=255)
                                                                                        aaaGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGGAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCgg                                                         >  6:221944/1‑149 (MQ=255)
                                                                                                taGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGt                                                   <  8:294186/147‑1 (MQ=255)
                                                                                                taGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGt                                                   >  7:294186/1‑147 (MQ=255)
                                                                                                          gACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAAttt                                         <  7:272655/147‑1 (MQ=255)
                                                                                                          gACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTACCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAAttt                                         <  6:308878/147‑1 (MQ=255)
                                                                                                           aCGGCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACCAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGc                                      <  7:78972/149‑1 (MQ=255)
                                                                                                             gCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCt                                    <  8:203304/149‑1 (MQ=255)
                                                                                                               cATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCtgt                                  >  6:208173/1‑149 (MQ=255)
                                                                                                               cATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCtgt                                  >  5:188095/1‑149 (MQ=255)
                                                                                                               cATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTAGGTAGAGCGCGGCACGGTGAAAGAATTTGGCCtgt                                  >  3:251822/1‑149 (MQ=255)
                                                                                                               cATTATCGCCCCCGATGACAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCATTACGCCAAATGTGATGCGCCCGTATGTAGAGCTCGGCACGGTTAAAGAATTTGGCCtgt                                  <  4:114708/149‑1 (MQ=255)
                                                                                                                   aTCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGa                              >  1:275858/1‑149 (MQ=255)
                                                                                                                    tCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGa                              <  8:299684/148‑1 (MQ=255)
                                                                                                                         cccGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGAGGCAAGGTGAAAGAATTTGGCCTGTGGGATGTGGt                        >  7:394480/1‑149 (MQ=255)
                                                                                                                         cccGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGGAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGGATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGt                        >  1:123979/1‑149 (MQ=255)
                                                                                                                          ccGATGCCACCGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGCCAAAGTCGCGATTGTCGGCTTCAGTACGCCAAATGTGCTGCGCCCGTATGTAGAGCGAGGCCCGGTGATAGCATTTGGCCTGTGGGATGTgggt                       >  7:297136/1‑147 (MQ=255)
                                                                                                                           cGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTc                      <  5:363030/149‑1 (MQ=255)
                                                                                                                              tGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGc                   <  7:201456/149‑1 (MQ=255)
                                                                                                                                 cAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGATAGAATTTGGCCTGTGGGATGTGGTTCAGCAAg                >  3:90831/1‑149 (MQ=255)
                                                                                                                                     gCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAATATGAGAAAGTAGCGAGTGTCGGATTCAGTCCGCCAAGTGTGATGCCCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGGGGGATGTGGTTCAGCAAGGCa             >  3:124865/1‑148 (MQ=255)
                                                                                                                                      cccTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCaa            >  6:220823/1‑148 (MQ=255)
                                                                                                                                       ccTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCaaaa          >  1:270862/1‑149 (MQ=255)
                                                                                                                                          gCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAAttt       >  2:165946/1‑149 (MQ=255)
                                                                                                                                           cccGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTc      <  8:394480/149‑1 (MQ=255)
                                                                                                                                           cccGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATCTGATGCGCCCGTATGTAGAGCCCGGCACGGTCAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAACGCAAAAttt       >  4:168189/1‑148 (MQ=255)
                                                                                                                                            ccGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCa     <  7:268638/149‑1 (MQ=255)
                                                                                                                                              gCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAgt   <  3:108981/149‑1 (MQ=255)
                                                                                                                                              gCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAgt   <  5:220823/149‑1 (MQ=255)
                                                                                                                                               cTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAgtg  >  1:247901/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAAAGAGCATCCAGGCTGGGAAATTGTCACTACGCAGTTTGGCTATAACGATGCCACTAAATCGTTACAAACCGCAGAAGGAATATTAAAAGCGTATAGCGATCTCGACGCCATTATCGCCCCCGATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTG  >  NC_000913/1605587‑1605878

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