Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,503,125 G→A *223* (TAG→TAA ydcL → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,503,1250GA98.7% 256.0 / NA 79*223* (TAG→TAAydcLlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/38);  minor base C (1/0);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGATG  >  NC_000913/1502991‑1503257
                                                                                                                                      |                                                                                                                                    
gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                                  >  2:1131500/1‑139 (MQ=255)
gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                                  >  1:2052301/1‑139 (MQ=255)
gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                                  >  2:589073/1‑139 (MQ=255)
        aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:1970205/139‑1 (MQ=255)
        aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:1842355/139‑1 (MQ=255)
        aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:1131500/139‑1 (MQ=255)
        aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:2052301/139‑1 (MQ=255)
             ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           >  2:1168217/1‑133 (MQ=255)
             ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           <  1:1168217/133‑1 (MQ=255)
             ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCg                                                                                                                     >  2:1160893/1‑139 (MQ=255)
                       cGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  2:907290/1‑124 (MQ=255)
                       cGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:907290/124‑1 (MQ=255)
                               aaGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCATCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTAt                                                                                                   >  2:374089/1‑139 (MQ=255)
                                     tGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:247759/110‑1 (MQ=255)
                                     tGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  1:247759/1‑110 (MQ=255)
                                            cGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             <  2:2004707/116‑1 (MQ=255)
                                            cGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             >  1:2004707/1‑116 (MQ=255)
                                                    cGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           <  2:1596417/94‑1 (MQ=255)
                                                    cGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           >  1:1596417/1‑94 (MQ=255)
                                                     gCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCAt                                                                             >  2:477572/1‑139 (MQ=255)
                                                          tGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTc                                                                        >  1:976192/1‑139 (MQ=255)
                                                           ggCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCa                                                                                                                        >  2:441923/1‑90 (MQ=255)
                                                           ggCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCa                                                                                                                        <  1:441923/90‑1 (MQ=255)
                                                              tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCATCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  1:374089/139‑1 (MQ=255)
                                                              tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  2:235648/139‑1 (MQ=255)
                                                              tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  1:477572/139‑1 (MQ=255)
                                                                 cGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTACTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCgg                                                                 >  1:1197799/1‑139 (MQ=255)
                                                                              cAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                                    <  2:976192/139‑1 (MQ=255)
                                                                                 gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATc                                                                                                  >  2:923607/1‑90 (MQ=255)
                                                                                 gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATc                                                                                                  <  1:923607/90‑1 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               <  1:423555/51‑1 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               >  2:423555/1‑51 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               <  1:536504/51‑1 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               >  2:536504/1‑51 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  2:1193166/110‑1 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    >  1:1193166/1‑110 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTg                                       >  2:1681708/1‑139 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTg                                       >  1:808483/1‑139 (MQ=255)
                                                                                           aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGa                                      >  1:596893/1‑139 (MQ=255)
                                                                                            cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               <  1:1814075/50‑1 (MQ=255)
                                                                                            cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               >  2:1814075/1‑50 (MQ=255)
                                                                                            cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           <  2:836593/54‑1 (MQ=255)
                                                                                            cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           >  1:836593/1‑54 (MQ=255)
                                                                                            cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                                  <  1:15671/111‑1 (MQ=255)
                                                                                            cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                                  >  2:15671/1‑111 (MQ=255)
                                                                                             aTGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAc                                     <  1:1160893/139‑1 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               <  1:1548864/47‑1 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               >  2:1548864/1‑47 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           >  1:2728/1‑51 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           <  2:2728/51‑1 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAc                                                                                                                              <  2:1102329/48‑1 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAc                                                                                                                              >  1:1102329/1‑48 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  2:317610/1‑52 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:317610/52‑1 (MQ=255)
                                                                                               ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                                   >  1:1878341/1‑139 (MQ=255)
                                                                                                gCGACCGATGCCACCATGTGTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAg                                  <  2:1878341/139‑1 (MQ=255)
                                                                                                  gACCGATGCCACCATGTTTGACGTGAACAAAAACTAATTGCAGACGCGCCATCGTTACATGACGCTTTTTATCAGGTGCGACGTAAAACCCCATCCTCCAGGGCGGGGATACAGGACGCGGTTTTGCACCTGACCAgga                                >  1:1356844/1‑138 (MQ=255)
                                                                                                   aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  1:1064553/102‑1 (MQ=255)
                                                                                                   aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    >  2:1064553/1‑102 (MQ=255)
                                                                                                   aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGgtg                               >  1:442835/1‑139 (MQ=255)
                                                                                                     cGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGtt                             >  2:89610/1‑139 (MQ=255)
                                                                                                     cGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGtt                             <  2:808483/139‑1 (MQ=255)
                                                                                                       aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTg                           >  1:1127766/1‑139 (MQ=255)
                                                                                                          ccaccaTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtg                        <  2:1197799/139‑1 (MQ=255)
                                                                                                           caccaTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtgt                       <  2:1127766/139‑1 (MQ=255)
                                                                                                                tGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             <  1:1969745/48‑1 (MQ=255)
                                                                                                                tGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             >  2:1969745/1‑48 (MQ=255)
                                                                                                                  tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                                          >  1:353832/1‑113 (MQ=255)
                                                                                                                  tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                                          <  2:353832/113‑1 (MQ=255)
                                                                                                                  tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtg                        <  2:458561/131‑1 (MQ=255)
                                                                                                                  tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtg                        >  1:458561/1‑131 (MQ=255)
                                                                                                                     gACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTg                           <  1:308268/125‑1 (MQ=255)
                                                                                                                     gACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTg                           >  2:308268/1‑125 (MQ=255)
                                                                                                                     gACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGt             >  2:1230272/1‑139 (MQ=255)
                                                                                                                          tAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTgg        <  1:1681708/139‑1 (MQ=255)
                                                                                                                           aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGc       >  2:558411/1‑139 (MQ=255)
                                                                                                                            aCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCg      >  1:23215/1‑139 (MQ=255)
                                                                                                                             cAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGa    <  2:596893/139‑1 (MQ=255)
                                                                                                                                aaaGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGATg  <  2:442835/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGATG  >  NC_000913/1502991‑1503257

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