Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,747 C→T P26S (CCG→TCG)  metB → O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,7470CT98.6% 224.3 / NA 70P26S (CCG→TCG) metBO‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/35);  minor base A (1/0);  total (35/35)
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

CAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAG  >  NC_000913/4128616‑4128882
                                                                                                                                   |                                                                                                                                       
agaTGTTTACTCTGGTGTCTGACATTTCACCAACACAGCCCAGGGAACCTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCAt                                                                                                                                 <  2:268218/138‑1 (MQ=255)
     tttACTCTTGTGCCTGCCATTTTACCGACAACGCCAAGGTAACTTCATAACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATTGATCCATCTTTc                                                                                                                            <  2:390292/140‑1 (MQ=255)
        aCTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAg                                                                                                                         <  1:532938/140‑1 (MQ=255)
        aCTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAg                                                                                                                         >  1:739871/1‑140 (MQ=255)
             ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCt                                                                                                                    >  1:593381/1‑140 (MQ=255)
             ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCt                                                                                                                    >  1:824157/1‑140 (MQ=255)
             ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCt                                                                                                                    >  2:795836/1‑140 (MQ=255)
             ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGGTGTGCCGTCGATCGATCTTTCCAGCACCt                                                                                                                    >  1:268218/1‑140 (MQ=255)
             ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGGATGGTTGGGTTGTCCCGTCGAGCCATCTTTCCAGCACCt                                                                                                                    >  2:219707/1‑140 (MQ=255)
              gTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGTCGAACAGTTTGGCTGCGTTGTCCCAACGATCCATCTTTCCAGCACCta                                                                                                                   >  2:361543/1‑140 (MQ=255)
              gTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAg                                                                                                                         >  1:445713/1‑134 (MQ=255)
              gTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAg                                                                                                                         <  2:445713/134‑1 (MQ=255)
              gTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCta                                                                                                                   >  2:624576/1‑140 (MQ=255)
              gTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTGGTCGCATCGATCCATCTTTCCAGCACCta                                                                                                                   >  1:330387/1‑140 (MQ=255)
               tGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCtat                                                                                                                  >  1:18790/1‑140 (MQ=255)
               tGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCtat                                                                                                                  >  1:156493/1‑140 (MQ=255)
               tGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCtat                                                                                                                  <  1:13764/140‑1 (MQ=255)
                   tGACACTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACt                                                                                                              >  1:631596/1‑140 (MQ=255)
                          tCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCgg                                                                                                       <  2:18790/140‑1 (MQ=255)
                                   aaGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGa                                                                                              >  1:143868/1‑140 (MQ=255)
                                         aGGGAACTTCATCACCTGACGCATACACACGCCACCACCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACcgcg                                                                                        <  1:249070/140‑1 (MQ=255)
                                         aGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACcgcg                                                                                        <  1:795836/140‑1 (MQ=255)
                                           ggAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACcgcgcg                                                                                      <  2:81992/140‑1 (MQ=255)
                                           ggAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACcgcgcg                                                                                      >  2:805272/1‑140 (MQ=255)
                                            gAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACcgcgcgc                                                                                     <  1:551707/140‑1 (MQ=255)
                                            gAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACcgcgcgc                                                                                     <  1:689270/140‑1 (MQ=255)
                                                ttCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGCTCCAACTTTCCAGTACCCATAGCTTTACCGAATTTCATGAACCGCGGGAGCac                                                                                 >  2:288524/1‑139 (MQ=255)
                                                ttCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCAGCTTTCCAGGACCAATAGCTTTACCGGATTTGATGAACCGCGCGCGCAt                                                                                 >  2:239858/1‑140 (MQ=255)
                                                     cacaTGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTa                                                                            >  1:183683/1‑140 (MQ=255)
                                                          gACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCCCCTATAACTTTACCGGAGTTAATGAACCGCGCGCGCATGCTTACtcgc                                                                       >  2:418555/1‑140 (MQ=255)
                                                          gACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtcgc                                                                       >  1:58526/1‑140 (MQ=255)
                                                              cgTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtcgcgtcg                                                                   <  2:183683/140‑1 (MQ=255)
                                                                 aaaCAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGAATTGAATGAACCGCGCGCGCAAGATTACTCGCGTCGCgg                                                                >  2:219325/1‑140 (MQ=255)
                                                                       gCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAAccc                                                          >  2:595816/1‑140 (MQ=255)
                                                                        ccaccaTCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCa                                                         >  1:62828/1‑140 (MQ=255)
                                                                             aTCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAAcgcg                                                    <  2:593381/140‑1 (MQ=255)
                                                                             aTCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAAcgcg                                                    <  1:478976/140‑1 (MQ=255)
                                                                               cGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAAcgcgcg                                                  <  2:156493/140‑1 (MQ=255)
                                                                                 cAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGAt                                                <  2:824157/140‑1 (MQ=255)
                                                                                 cAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGAt                                                <  2:549436/140‑1 (MQ=255)
                                                                                   gTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGAtgt                                              >  1:574854/1‑140 (MQ=255)
                                                                                    tGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCCTGATTACGCGCTTCTCGGCAACCCAACGCGCGAtgtg                                             >  1:374265/1‑140 (MQ=255)
                                                                                      cGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGAGTCCCGGCGTCGCTGCAATACAACGCGTGATGTGGt                                           >  2:417887/1‑140 (MQ=255)
                                                                                        tAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAAcgcgcg                                                  >  1:785094/1‑131 (MQ=255)
                                                                                        tAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAAcgcgcg                                                  <  2:785094/131‑1 (MQ=255)
                                                                                            gggTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCg                                     >  1:509871/1‑140 (MQ=255)
                                                                                              gTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCTTCGAGGCAACCCTACGCACGATGTGGTTCAGCGTg                                   >  2:349221/1‑140 (MQ=255)
                                                                                              gTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCCCGTCGCGCCAACCCAACGCGCAATGTGGTTCAGCGTg                                   >  1:334321/1‑140 (MQ=255)
                                                                                                tAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTgc                                  <  2:477136/139‑1 (MQ=255)
                                                                                                tAAATGACGACGAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTgc                                  >  1:477136/1‑139 (MQ=255)
                                                                                                      acgacgAACAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCa                           <  2:143868/140‑1 (MQ=255)
                                                                                                            aaCAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTgg                             >  2:634483/1‑132 (MQ=255)
                                                                                                            aaCAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTgg                             <  1:634483/132‑1 (MQ=255)
                                                                                                            aaCAGTATGGTTGCGATGTCCCTTCTATCGATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTg                     <  1:239858/140‑1 (MQ=255)
                                                                                                              cAGTATGGTTGCGTTGTTCCATCGATCCATCTTTCTAGCGCCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCTATGTGGTTCAGCGTGCGCTGGCAGAACTGGa                   <  2:334321/140‑1 (MQ=255)
                                                                                                              cAGTATGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGa                   <  2:574854/140‑1 (MQ=255)
                                                                                                                  aTGGTTGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAggt               <  2:631596/140‑1 (MQ=255)
                                                                                                                      ttGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGgtggtg           >  2:534173/1‑140 (MQ=255)
                                                                                                                      ttGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGAGGTTCAGCGTGCACTGGGAGAACTGGAAGgtggtg           >  2:332550/1‑140 (MQ=255)
                                                                                                                       tGCGTTGTCCCATCGATTCATCTTTCCAGCACCTACAACATTACCGGACTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgc          <  2:374265/140‑1 (MQ=255)
                                                                                                                       tGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgc          <  2:624541/140‑1 (MQ=255)
                                                                                                                       tGCGTTGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgc          <  1:534173/140‑1 (MQ=255)
                                                                                                                         cGTTGTCCCATCGATTCATCTTTCCTGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctg        <  1:349221/140‑1 (MQ=255)
                                                                                                                           ttGTCCCATCGATCCGTCTTTCCAGTACCCATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggt      <  1:361543/140‑1 (MQ=255)
                                                                                                                           ttGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctg        >  1:762055/1‑138 (MQ=255)
                                                                                                                           ttGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctg        <  2:762055/138‑1 (MQ=255)
                                                                                                                            tGTCTCATCGATCCATCCTTCCTGCACCTATTACTTTACCGGATTTAAGGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg     <  2:330387/140‑1 (MQ=255)
                                                                                                                            tGTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg     <  1:805272/140‑1 (MQ=255)
                                                                                                                            tGTCCCATCGATCCAGCTTTCCGGCACCTGTAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg     <  1:288524/140‑1 (MQ=255)
                                                                                                                            tGTCCCATCGATCCACCTTTCCCGCACCTATAACTTTTCCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg     <  1:219707/140‑1 (MQ=255)
                                                                                                                             gTCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtgc    >  2:564038/1‑140 (MQ=255)
                                                                                                                              tCCCATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCa   <  2:62828/140‑1 (MQ=255)
                                                                                                                               cccATCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAg  >  2:3966/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAG  >  NC_000913/4128616‑4128882

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: