Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,206,673 T→A 5.7% C71S (TGT→AGT) ‡ 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,6730TA 5.7% 245.2 / 5.2 70C71S (TGT→AGT) ‡ttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (39/27);  new base A (2/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

ACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTG  >  NC_000913/3206526‑3206817
                                                                                                                                                   |                                                                                                                                                
aCCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGAAAAAAGCGATTGACCTGAATCGTcctgcctg                                                                                                                                                 <  2:87191/149‑1 (MQ=255)
   gCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCa                                                                                                                                              >  4:54765/1‑149 (MQ=255)
        gATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacg                                                                                                                                         >  4:216902/1‑148 (MQ=255)
        gATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacg                                                                                                                                         >  4:574589/1‑148 (MQ=255)
         aTTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacga                                                                                                                                        <  4:616119/149‑3 (MQ=255)
         aTTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGACAGGCAGGacga                                                                                                                                        <  3:574589/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:135323/1‑149 (MQ=255)
                     aTGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTAt                                                                                                                            <  4:608230/149‑1 (MQ=255)
                        ccTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGtt                                                                                                                         <  2:268337/149‑1 (MQ=255)
                           gatgaCGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGttctt                                                                                                                      <  2:607087/149‑1 (MQ=255)
                                  tggtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTaaa                                                                                                               >  1:168450/1‑149 (MQ=255)
                                      ggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCg                                                                                                           <  1:484115/149‑1 (MQ=255)
                                      ggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCg                                                                                                           >  2:484115/1‑149 (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                                                                                        >  4:235152/1‑149 (MQ=255)
                                                         aaGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCt                                                                                        >  3:79102/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                                                                     >  3:58321/1‑132 (MQ=255)
                                                                                             tACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATa                                                                     <  4:58321/132‑1 (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:200244/1‑149 (MQ=255)
                                                                                                                  aTGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGcc                               >  1:282597/1‑149 (MQ=255)
                                                                                                                    gCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCAc                             <  4:74500/149‑1 (MQ=255)
                                                                                                                    gCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCAc                             <  4:358400/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               >  4:66382/1‑136 (MQ=255)
                                                                                                                                               tgcctgTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCg               <  3:66382/136‑1 (MQ=255)
                                                                                                                                               tgcctgTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTCCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTg  >  2:528816/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTGTCAGGACACCGGGGAGATTATGTTCTTCGTTAAAGTCGGTTCCCGCTTCCCACTGCTTGGCGAGCTGCAAAGCATACTCAAACAAGCCGTGGAAGAGGCAACCGTCAAAGCGCCACTACGTCACAATGCGGTAGAAATTTTTG  >  NC_000913/3206526‑3206817

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: ‑