Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 397,915 T→C V15A (GTT→GCT)  yaiW → surface‑exposed outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913397,9150TC98.3% 210.8 / NA 60V15A (GTT→GCT) yaiWsurface‑exposed outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/30);  minor base G (1/0);  total (30/30)
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

TCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGAT  >  NC_000913/397779‑398048
                                                                                                                                        |                                                                                                                                     
tctatctaCAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGctct                                                                                                                                    <  1:915347/140‑1 (MQ=255)
        cAAACGTTTGCGTAGCTATGAACTTAAGCTGGATGGTGACAAATTTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAg                                                                                                                            <  1:627394/140‑1 (MQ=255)
        cAAACGTCTGCGCAGCTTTGAACATGAGCTGGATTGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAg                                                                                                                            >  2:955518/1‑140 (MQ=255)
        cAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAg                                                                                                                            <  2:586216/140‑1 (MQ=255)
        cAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAg                                                                                                                            <  1:543569/140‑1 (MQ=255)
          aaCGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCg                                                                                                                          <  1:564584/140‑1 (MQ=255)
            cGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCggg                                                                                                                        <  2:798823/140‑1 (MQ=255)
            cGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCggg                                                                                                                        <  2:1013017/140‑1 (MQ=255)
            cGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCggg                                                                                                                        >  1:659518/1‑140 (MQ=255)
              tCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCggggg                                                                                                                      >  1:961818/1‑139 (MQ=255)
              tCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGgt                                                                                                                      >  2:727080/1‑140 (MQ=255)
              tCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCCTCGCTGTAATTGTTGGCGGCTATGGTTTTAGCGGGgt                                                                                                                      >  1:404977/1‑140 (MQ=255)
               cTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGgtg                                                                                                                     <  2:192583/140‑1 (MQ=255)
                tGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGgtgt                                                                                                                    <  2:136517/140‑1 (MQ=255)
                tGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGgtgt                                                                                                                    >  2:731896/1‑140 (MQ=255)
                tGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCAATGTAGGCGGCTCTGGTTTTAGCGGGgtgt                                                                                                                    >  1:270903/1‑140 (MQ=255)
                       cTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTa                                                                                                                             >  1:750604/1‑124 (MQ=255)
                       cTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTa                                                                                                                             <  2:750604/124‑1 (MQ=255)
                        tttGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCa                                                                                                            <  2:82757/140‑1 (MQ=255)
                          tGAACATGAGCTGGATGTTGACAAAATTCAGTAAGTAACTCATACCTTTAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAgg                                                                                                          <  2:404977/140‑1 (MQ=255)
                          tGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAgg                                                                                                          >  1:111788/1‑140 (MQ=255)
                                 gAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTGATTGTTGGTGGCGCTGCTTTTAGCGGGGTGAAGTAGTCAGGCTCCGCa                                                                                                   >  2:382707/1‑140 (MQ=255)
                                 gAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCa                                                                                                   >  1:548279/1‑140 (MQ=255)
                                 gAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCa                                                                                                   <  2:725789/140‑1 (MQ=255)
                                  aGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCaa                                                                                                  <  2:1002519/140‑1 (MQ=255)
                                   gCTGGATGGTGACAATATTCAGTAACTAACCCATACCATGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAAc                                                                                                 <  1:264939/140‑1 (MQ=255)
                                    cTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATGGTTGGCGGGTCTGGTTGTAGCGGGGAGTAGAAGTCAGGCTCCGCAAcc                                                                                                >  2:284076/1‑140 (MQ=255)
                                               cAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAgg                                                                                     >  2:921276/1‑140 (MQ=255)
                                                 aaaTTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGc                                                                                   >  2:1000780/1‑140 (MQ=255)
                                                    ttCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGaa                                                                                >  1:589294/1‑140 (MQ=255)
                                                      cAGGAAGTAACCCATACCTTGAGCTAAAAAGGGGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGaaaa                                                                              >  1:579003/1‑140 (MQ=255)
                                                       aGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGaaaaa                                                                             >  2:797296/1‑140 (MQ=255)
                                                        ggAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAg                                                                            >  1:867384/1‑140 (MQ=255)
                                                           aGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGa                                                                         <  2:659518/140‑1 (MQ=255)
                                                           aGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGa                                                                         <  1:727080/140‑1 (MQ=255)
                                                           aGTAACCCATACCTTCAGCTAACAAAGAGCAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGa                                                                         <  1:264908/140‑1 (MQ=255)
                                                              aaCCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAg                                                                      >  2:691440/1‑140 (MQ=255)
                                                                cccATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGAt                                                                    <  2:867384/140‑1 (MQ=255)
                                                                          gAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTagcgagcg                                                          >  1:540499/1‑140 (MQ=255)
                                                                          gAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGCCTCCGCAACCGTTGAAGAAGGGCGAAGAAGCCATAGATGCagcgagcg                                                          >  2:457768/1‑140 (MQ=255)
                                                                           aGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGcgt                                                         >  2:237255/1‑140 (MQ=255)
                                                                           aGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGcgt                                                         >  1:889508/1‑140 (MQ=255)
                                                                                aaaaGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGc                                                    >  2:637794/1‑140 (MQ=255)
                                                                                aaaaGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGc                                                    >  2:65838/1‑140 (MQ=255)
                                                                                      agaAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGCGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCaga                                              <  2:589294/140‑1 (MQ=255)
                                                                                        aaGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCagaag                                            <  1:797296/140‑1 (MQ=255)
                                                                                        aaGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCagaag                                            <  1:65838/140‑1 (MQ=255)
                                                                                        aaGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCagaag                                            <  1:637794/140‑1 (MQ=255)
                                                                                                tCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGc                                    <  1:731896/140‑1 (MQ=255)
                                                                                                  gcgcgTCAATCCTCTATCTTCGGTGACATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGa                                  <  1:457768/140‑1 (MQ=255)
                                                                                                  gcgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGa                                  <  2:540499/140‑1 (MQ=255)
                                                                                                  gcgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGa                                  >  1:553470/1‑140 (MQ=255)
                                                                                                   cgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCATGCGGCAGAAGATGCCTGCGAg                                 >  1:761563/1‑140 (MQ=255)
                                                                                                        cAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCagaag                                            >  1:237547/1‑124 (MQ=255)
                                                                                                        cAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCagaag                                            <  2:237547/124‑1 (MQ=255)
                                                                                                        cAATCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGa                            <  2:548279/140‑1 (MQ=255)
                                                                                                           tCCTCTATCTTCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAg                         <  1:134819/140‑1 (MQ=255)
                                                                                                                     tCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGc               <  1:218449/140‑1 (MQ=255)
                                                                                                                     tCGCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGc               >  2:218449/1‑140 (MQ=255)
                                                                                                                       gCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCt             >  1:192098/1‑140 (MQ=255)
                                                                                                                       gCTGTCATTGTTGGCGGCTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTAGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCt             >  1:827181/1‑140 (MQ=255)
                                                                                                                                  tgGCGGCTCTGGTTTCAGCGGGCTGTAGTTGTCAGGCACCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGAt  <  1:382707/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
TCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGAT  >  NC_000913/397779‑398048

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