Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,206,675 T→G 5.1% C71W (TGT→TGG) ‡ ttdA → L‑tartrate dehydratase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,6750TG 5.1% 271.8 / 6.3 78C71W (TGT→TGG) ‡ttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (42/32);  new base G (2/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

GCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTGACGAA  >  NC_000913/3206529‑3206822
                                                                                                                                                  |                                                                                                                                                   
gCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCa                                                                                                                                                   >  4:54765/1‑149 (MQ=255)
     gATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacg                                                                                                                                              >  4:574589/1‑148 (MQ=255)
     gATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacg                                                                                                                                              >  4:216902/1‑148 (MQ=255)
      aTTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacga                                                                                                                                             <  3:574589/149‑3 (MQ=255)
      aTTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacga                                                                                                                                             <  4:616119/149‑3 (MQ=255)
           tACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCgggg                                                                                                                                        >  1:25260/1‑149 (MQ=255)
              ccGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGgaga                                                                                                                                     >  1:171959/1‑149 (MQ=255)
              ccGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGgaga                                                                                                                                     <  1:562929/149‑1 (MQ=255)
                  aTGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTAt                                                                                                                                 <  4:608230/149‑1 (MQ=255)
                  aTGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTAt                                                                                                                                 >  4:135323/1‑149 (MQ=255)
                     ccTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGtt                                                                                                                              <  2:268337/149‑1 (MQ=255)
                        gatgaCGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGttctt                                                                                                                           <  2:607087/149‑1 (MQ=255)
                               tggtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTaaa                                                                                                                    >  1:168450/1‑149 (MQ=255)
                                   ggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCg                                                                                                                >  2:484115/1‑149 (MQ=255)
                                   ggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCg                                                                                                                <  1:484115/149‑1 (MQ=255)
                                    gATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCgg                                                                                                               <  3:335615/149‑1 (MQ=255)
                                     aTAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGt                                                                                                              >  2:13751/1‑149 (MQ=255)
                                      tAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTAAATCGTCCAGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGtt                                                                                                             >  4:504947/1‑149 (MQ=255)
                                         aCTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTccc                                                                                                          >  1:173027/1‑149 (MQ=255)
                                         aCTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTccc                                                                                                          >  2:587373/1‑149 (MQ=255)
                                          cTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCg                                                                                                         >  1:184609/1‑149 (MQ=255)
                                           tAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGc                                                                                                        >  4:421965/1‑149 (MQ=255)
                                             aaaCAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAGAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCtg                                                                                                      >  2:446880/1‑148 (MQ=255)
                                             aaaCAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCtt                                                                                                      <  2:540286/149‑1 (MQ=255)
                                              aaCAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTc                                                                                                     >  2:204178/1‑149 (MQ=255)
                                                    tAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTg                                                                                               >  1:402487/1‑149 (MQ=255)
                                                     aaaGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGc                                                                                              <  2:168450/149‑1 (MQ=255)
                                                      aaGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCt                                                                                             >  3:79102/1‑149 (MQ=255)
                                                      aaGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCt                                                                                             >  4:235152/1‑149 (MQ=255)
                                                         gATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTgg                                                                                          >  4:535111/1‑149 (MQ=255)
                                                              cGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGc                                                                                     <  1:177276/149‑1 (MQ=255)
                                                              cGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGc                                                                                     <  4:356576/149‑1 (MQ=255)
                                                               gAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCt                                                                                    >  3:74500/1‑149 (MQ=255)
                                                                aaaCGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGttcttc                                                                                                                          >  1:62658/1‑110 (MQ=255)
                                                                aaaCGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCCATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGttcttc                                                                                                                          <  2:62658/110‑1 (MQ=255)
                                                                 aaCGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAAACGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGc                                                                                  >  4:265222/1‑149 (MQ=255)
                                                                  aCGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCa                                                                                 >  3:217467/1‑149 (MQ=255)
                                                                    gtcgtcGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCaaa                                                                               >  3:262649/1‑149 (MQ=255)
                                                                     tcgtcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAg                                                                              >  3:201112/1‑149 (MQ=255)
                                                                       gtcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTATTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCa                                                                            >  3:358400/1‑149 (MQ=255)
                                                                          gATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATAc                                                                         >  2:435586/1‑149 (MQ=255)
                                                                          gATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATAc                                                                         <  3:23676/149‑1 (MQ=255)
                                                                               ggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCaaa                                                                    >  4:268668/1‑149 (MQ=255)
                                                                                gAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAAc                                                                   >  4:95319/1‑149 (MQ=255)
                                                                                gAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAAc                                                                   >  2:15689/1‑149 (MQ=255)
                                                                                  aaaTTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACaa                                                                 <  1:325869/149‑1 (MQ=255)
                                                                                    aTTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGc                                                               <  2:25260/149‑1 (MQ=255)
                                                                                    aTTATCTACCATACGATGTTCGACAACATGCAAAAAGCGAGTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGc                                                               >  4:359230/1‑149 (MQ=255)
                                                                                       aTCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGt                                                            <  4:79102/149‑1 (MQ=255)
                                                                                        tCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTg                                                           <  2:171959/149‑1 (MQ=255)
                                                                                          tACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATa                                                                          <  4:58321/132‑1 (MQ=255)
                                                                                          tACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATa                                                                          >  3:58321/1‑132 (MQ=255)
                                                                                          tACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGa                                                         >  2:315892/1‑149 (MQ=255)
                                                                                             cATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAaga                                                      <  1:386073/149‑1 (MQ=255)
                                                                                               tACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAgg                                                    >  3:224320/1‑149 (MQ=255)
                                                                                               tACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAgg                                                    <  3:605550/149‑1 (MQ=255)
                                                                                                  gATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCaa                                                 >  4:61115/1‑149 (MQ=255)
                                                                                                       tCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTc                                            <  1:15689/149‑1 (MQ=255)
                                                                                                         gACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGttcttc                                                                                                                          >  2:184042/1‑69 (MQ=255)
                                                                                                         gACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGttcttc                                                                                                                          <  1:184042/69‑1 (MQ=255)
                                                                                                           caacaTGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAg                                        <  2:402487/149‑1 (MQ=255)
                                                                                                            aacaTGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAgc                                       <  1:63109/149‑1 (MQ=255)
                                                                                                               aTGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGcc                                    >  1:282597/1‑149 (MQ=255)
                                                                                                               aTGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGcc                                    >  1:200244/1‑149 (MQ=255)
                                                                                                                 gCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCAc                                  <  4:358400/149‑1 (MQ=255)
                                                                                                                 gCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCAc                                  <  4:74500/149‑1 (MQ=255)
                                                                                                                              gACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGc                     <  3:359230/149‑1 (MQ=255)
                                                                                                                                  tGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTa                 <  3:504947/149‑1 (MQ=255)
                                                                                                                                    aaTCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGa               >  2:104140/1‑149 (MQ=255)
                                                                                                                                     aTCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGaa              >  3:178402/1‑149 (MQ=255)
                                                                                                                                           ctgcctgTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAAttttt        <  4:262649/149‑1 (MQ=255)
                                                                                                                                            tgcctgTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCg                    <  3:66382/136‑1 (MQ=255)
                                                                                                                                            tgcctgTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCg                    >  4:66382/1‑136 (MQ=255)
                                                                                                                                            tgcctgTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTCCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTg       >  2:528816/1‑149 (MQ=255)
                                                                                                                                                 gTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTGACGaa  <  3:54765/149‑1 (MQ=255)
                                                                                                                                                 gTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTGACGaa  <  3:421965/149‑1 (MQ=255)
                                                                                                                                                 gTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTGACGaa  >  4:62906/1‑149 (MQ=255)
                                                                                                                                                 gTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTGACGaa  <  3:268668/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTGACGAA  >  NC_000913/3206529‑3206822

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