Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,561 Δ1 bp 31.1% intergenic (‑217/+73) yhfZ ← / ← trpS putative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,5600T.31.1% 207.4 / 103.6 90intergenic (‑216/+74)yhfZ/trpSputative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (35/27);  new base . (13/15);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATC  >  NC_000913/3512461‑3512658
                                                                                                   |                                                                                                  
aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt                                                                                                    >  2:2151453/1‑100 (MQ=255)
 tAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt                                                                                                   >  2:769250/1‑100 (MQ=255)
 tAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt                                                                                                   >  2:748116/1‑100 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                                 >  1:644528/1‑97 (MQ=255)
     aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACtt                                                                                               >  1:588095/1‑100 (MQ=255)
     aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACttt                                                                                              <  2:1852021/100‑1 (MQ=255)
           gggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGc                                                                                        >  1:1035198/1‑100 (MQ=255)
           gggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGc                                                                                        >  2:155931/1‑100 (MQ=255)
            ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGcc                                                                                       >  1:2146497/1‑100 (MQ=255)
              tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGcc                                                                                       >  1:901930/1‑99 (MQ=255)
              tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCgg                                                                                     <  1:507902/100‑1 (MQ=255)
                tCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCCATATGTCCATTACTTTGTAGGcc                                                                                       >  2:1781979/1‑97 (MQ=255)
                tCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGa                                                                                    <  1:1070599/100‑1 (MQ=255)
                 ccccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGAt                                                                                   >  1:269694/1‑100 (MQ=255)
                 ccccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATa                                                                                  >  1:376534/1‑100 (MQ=255)
                  cccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATa                                                                                  <  2:1189573/100‑1 (MQ=255)
                  cccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATaa                                                                                 >  1:879767/1‑100 (MQ=255)
                  cccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATa                                                                                  >  1:2105678/1‑99 (MQ=255)
                   ccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATaa                                                                                 >  2:2041515/1‑100 (MQ=255)
                     ggTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAgg                                                                               <  1:1602524/100‑1 (MQ=255)
                        ttCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGt                                                                            >  2:1330586/1‑100 (MQ=255)
                        ttCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGt                                                                            >  2:1984092/1‑100 (MQ=255)
                         tCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTc                                                                          <  1:2120928/100‑1 (MQ=255)
                         tCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTc                                                                          <  1:846646/100‑1 (MQ=255)
                          cAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTc                                                                          >  2:525093/1‑100 (MQ=255)
                           agagTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCa                                                                         <  1:213996/100‑1 (MQ=255)
                            gagTGATGATAAAAGCAAAATTGCCTTATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCAc                                                                        >  1:800057/1‑100 (MQ=255)
                            gagTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcg                                                                       >  1:2190014/1‑100 (MQ=255)
                                gatgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                    <  2:1649167/100‑1 (MQ=255)
                                 atgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                  <  2:644528/100‑1 (MQ=255)
                                    atAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCAcgccgaa                                                                  >  1:2237869/1‑96 (MQ=255)
                                      aaaaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTCCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                               <  2:2063620/99‑1 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                           <  2:801827/100‑1 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                           >  2:1527855/1‑100 (MQ=255)
                                         aGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                           <  2:282332/100‑1 (MQ=255)
                                          gCAAAATTGCCTTATGCGTTACGCTTATCAGGCCTACATTTCCTTGCAAGATGTGCATTACTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                           <  1:1761532/99‑2 (MQ=21)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                          >  2:87946/1‑100 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                         >  1:710423/1‑100 (MQ=255)
                                           cAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                        <  1:1499442/100‑1 (MQ=255)
                                            aaaaTTGCCGGATGCGCTACTCTTATCAGGCCTACATTACATTGCAATATGTGCATTACTTTGTAGGCCAGATAAGGCGATCACGCCGCATCCGGc                                                            >  2:1369230/1‑96 (MQ=37)
                                             aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                       <  1:1722481/100‑1 (MQ=255)
                                                ttGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaaa                                                   <  1:1783276/100‑1 (MQ=255)
                                                  gCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                  >  2:193960/1‑100 (MQ=255)
                                                  gCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCCGCATGAACAAAg                                                  <  1:1472493/100‑1 (MQ=255)
                                                     tGATGCGCTACGCTTATCAGGTCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgcgc                                               <  1:679358/100‑1 (MQ=255)
                                                      gATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCa                                              >  2:1518524/1‑100 (MQ=255)
                                                      gATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCa                                              >  1:1743089/1‑100 (MQ=255)
                                                       aTGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCaa                                             >  1:716751/1‑100 (MQ=255)
                                                         gcgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAtt                                           <  2:1049265/100‑1 (MQ=255)
                                                           gcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTg                                         <  2:1642706/100‑1 (MQ=255)
                                                           gcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTg                                         <  1:922091/100‑1 (MQ=255)
                                                           gcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAttt                                          >  1:1844736/1‑98 (MQ=255)
                                                           gcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGc                                        >  2:1444044/1‑100 (MQ=255)
                                                             tACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCa                                      <  2:1915891/100‑1 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                     <  2:1199458/100‑1 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                     <  1:2151468/100‑1 (MQ=255)
                                                                 cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                    <  1:1961319/99‑1 (MQ=255)
                                                                 cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCa                                   >  2:1257041/1‑100 (MQ=255)
                                                                 cTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCaa                                  <  1:1115619/100‑1 (MQ=255)
                                                                   tATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCTAGCAAt                                 >  1:412188/1‑100 (MQ=255)
                                                                   tATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                 <  1:882104/100‑1 (MQ=255)
                                                                    aTCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGAAATTTGCCAGTAATa                                >  1:138700/1‑100 (MQ=255)
                                                                     tCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAg                               <  1:2062213/100‑1 (MQ=255)
                                                                      cAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGt                              >  2:619663/1‑100 (MQ=255)
                                                                       aGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTg                             <  1:1110433/100‑1 (MQ=255)
                                                                        ggCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGa                            >  2:1586080/1‑100 (MQ=255)
                                                                        ggCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGa                            >  2:1352104/1‑100 (MQ=255)
                                                                        ggCCTACATTTCCTTACAATATGTGCATTACTTTGTAGGCCGGATAAGGCGATCACGCCACATACGTCATGCAGAAAGTGCAATTTGGCAGCAATAg                               >  2:854136/1‑97 (MQ=21)
                                                                          ccTACATTTCCTTGCAATATGTGCATACTTTGTAAGCCGGATAAGGCGTTCACGCCCCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGatgg                         >  2:577818/1‑97 (MQ=25)
                                                                           cTACATTTCCTTGCAATATGTGCATTACTTTGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGaa                           >  2:169686/1‑98 (MQ=255)
                                                                            tACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTa                        >  1:592568/1‑100 (MQ=255)
                                                                            tACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTa                        >  2:1581528/1‑100 (MQ=255)
                                                                            tACATTTCCTTGCAATAGGTGCATATTTTGTAGGTCGGATAATGCGTTCACGCTTCATCCGGCATTAACAAATTGTAATTTGCCAGCAATAGTGAATTa                        <  2:854838/99‑1 (MQ=21)
                                                                             aCATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCTATTTGCCAGCAATAGTGAATTAc                       <  2:708/100‑1 (MQ=255)
                                                                             aCATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTAc                       >  2:197795/1‑100 (MQ=255)
                                                                             aCATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTAc                       <  1:619663/100‑1 (MQ=255)
                                                                              cATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCCATAATGAATTAc                       >  1:1304061/1‑99 (MQ=255)
                                                                                tttCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCt                   <  2:1035196/100‑1 (MQ=255)
                                                                                  tCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCtt                  >  1:2104927/1‑100 (MQ=255)
                                                                                   ccTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg                <  1:883628/100‑1 (MQ=255)
                                                                                    cTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg                >  2:1448323/1‑100 (MQ=255)
                                                                                     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg                <  1:1787283/99‑1 (MQ=255)
                                                                                     ttGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc               <  2:2146482/99‑1 (MQ=255)
                                                                                      tGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGcc              <  2:1325500/100‑1 (MQ=255)
                                                                                       gCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCca             >  1:242386/1‑100 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca           >  1:1834705/1‑100 (MQ=255)
                                                                                            atGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAc       <  2:1900250/100‑1 (MQ=255)
                                                                                             tgtgCATTACTTTGTAGGCCTGATAAGGCGTTCACGCCGTATCCGGCATGAACAAAGCGCAATTTGCTAGCAATAGTGAATTACGGCTTCGCCACAAAAc       <  2:363944/100‑1 (MQ=255)
                                                                                             tgtgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAc       >  1:1481709/1‑100 (MQ=255)
                                                                                             tgtgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAc       >  1:1052554/1‑100 (MQ=255)
                                                                                                gCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCaa    >  2:248085/1‑100 (MQ=255)
                                                                                                gCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCaa    >  2:851131/1‑100 (MQ=255)
                                                                                                  aTTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGCACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATc  <  1:660372/100‑1 (MQ=255)
                                                                                                  aTTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATc  <  2:1663343/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATC  >  NC_000913/3512461‑3512658

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