Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,086,147 G→A A44A (GCG→GCA fdhD → formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,086,1470GA97.6% 289.2 / ‑3.1 85A44A (GCG→GCAfdhDformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (48/34);  total (51/34)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCGCTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGTAATGGT  >  NC_000913/4086003‑4086294
                                                                                                                                                |                                                                                                                                                   
tGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCGCTGGTc                                                                                                                                               >  1:663435/1‑151 (MQ=255)
   ttAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGAGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTa                                                                                                                                             >  2:932007/1‑150 (MQ=255)
      agagAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGTTGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAAc                                                                                                                                         >  1:319095/1‑151 (MQ=255)
          aaaGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGc                                                                                                                                      >  2:1367996/1‑150 (MQ=255)
          aaaGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGc                                                                                                                                      >  1:417797/1‑150 (MQ=255)
           aaGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCAt                                                                                                                                    <  1:936902/151‑1 (MQ=255)
            aGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACCACGGCAtt                                                                                                                                   >  2:866473/1‑151 (MQ=255)
             gTGAAAAAAACACAGCGAAAAGAAATCGAATATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCAttt                                                                                                                                  <  2:775021/151‑1 (MQ=255)
               gAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTc                                                                                                                                 <  2:99403/150‑1 (MQ=255)
                 aaaaaaCACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGc                                                                                                                               >  2:332186/1‑150 (MQ=255)
                 aaaaaaCACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGAAGATTTACAACACCCACGGTTAGTTGAGGTGGCGGAAGAAGCTCCCGGTGCACTGGTCTACAACGGAATTTTGc                                                                                                                               >  1:482030/1‑150 (MQ=255)
                   aaaaCACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCAtg                                                                                                                            <  1:217691/151‑1 (MQ=255)
                    aaaCACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCAtgt                                                                                                                           >  2:768686/1‑151 (MQ=255)
                    aaaCACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCAtgt                                                                                                                           <  2:785140/151‑1 (MQ=255)
                     aaCACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCAtgtg                                                                                                                          <  2:186992/151‑1 (MQ=255)
                         caGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATgtggtg                                                                                                                       <  2:582455/150‑1 (MQ=255)
                          aGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGtgat                                                                                                                     <  2:206848/151‑1 (MQ=255)
                            cGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGtgatga                                                                                                                   <  2:273178/151‑1 (MQ=255)
                            cGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGtgatg                                                                                                                    >  2:652022/1‑150 (MQ=255)
                            cGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGtgatg                                                                                                                    >  2:459741/1‑150 (MQ=255)
                             gAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTgatgat                                                                                                                  >  2:61135/1‑151 (MQ=255)
                               aaaGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATgg                                                                                                                <  2:399331/151‑1 (MQ=255)
                                aaGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATgga                                                                                                               <  1:315715/151‑2 (MQ=255)
                                aaGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGc                                                                                                               <  1:773304/151‑1 (MQ=255)
                                  gAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGt                                                                                                             <  2:933327/151‑1 (MQ=255)
                                     aTCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGc                                                                                                          >  2:769025/1‑151 (MQ=255)
                                       cGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGccc                                                                                                        >  1:659774/1‑151 (MQ=255)
                                        gAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGccc                                                                                                        >  2:347383/1‑150 (MQ=255)
                                        gAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCa                                                                                                       <  1:96470/151‑1 (MQ=255)
                                          aaaTGAGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGAAGATTTCCAAGACACACGGTTAGATGAGGCCACGGAAGGAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGccc                                                                                                        >  1:389632/1‑148 (MQ=255)
                                            aTGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGTATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGa                                                                                                   <  1:872196/151‑1 (MQ=255)
                                            aTGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCaa                                                                                                      >  2:685450/1‑148 (MQ=255)
                                              gtgACAAATATCACAGTTGTTCGTCAAATTGAGTTATGGCGGGGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTACCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCaaa                                                                                                     <  1:1162921/147‑1 (MQ=255)
                                              gtgACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGAcc                                                                                                 >  2:109066/1‑151 (MQ=255)
                                                gACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGTGTGGTGATTTACAAAACCCACAGATAGATGAGGTTGCGGAAAAAGAACCCGTTGCACTGGGCTACAACGGCATTTTGCATGTGGTgatgatcgggccgctca                                                                                                       >  1:1036894/1‑132 (MQ=255)
                                                gACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCtt                                                                                               >  2:723629/1‑151 (MQ=255)
                                                  cAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGa                                                                                             <  2:985225/151‑1 (MQ=255)
                                                   aaaTATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATATACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTg                                                                                              >  2:1061527/1‑149 (MQ=255)
                                                     atatCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTa                                                                                          <  1:775253/151‑1 (MQ=255)
                                                        tCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTAc                                                                                         >  1:1328418/1‑149 (MQ=255)
                                                          acaGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAAAGGCATTTCGCATGTGGGGATGATGGCGTCGCCCAAAGACCTTGAGGACttt                                                                                      >  1:779645/1‑150 (MQ=255)
                                                             ggTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGGGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTgcgc                                                                                  >  2:710414/1‑151 (MQ=255)
                                                                 ttCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCgg                                                                              >  1:949565/1‑151 (MQ=255)
                                                                 ttCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCgg                                                                              >  1:769839/1‑151 (MQ=255)
                                                                 ttCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGGTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCgg                                                                              >  2:1335643/1‑151 (MQ=255)
                                                                    gTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGttt                                                                           <  2:540628/151‑1 (MQ=255)
                                                                      cAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGttttt                                                                         >  2:760076/1‑151 (MQ=255)
                                                                       aaaTTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGttttt                                                                         <  1:822872/150‑1 (MQ=255)
                                                                        aaTTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCg                                                                       >  2:324906/1‑151 (MQ=255)
                                                                            gAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCttt                                                                   >  1:214414/1‑151 (MQ=255)
                                                                              gTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCCCTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTcc                                                                 >  2:1228161/1‑151 (MQ=255)
                                                                               ttATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCg                                                                >  1:73282/1‑151 (MQ=255)
                                                                                tATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCg                                                                >  2:629564/1‑150 (MQ=255)
                                                                                 aTGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGa                                                               <  1:929931/150‑1 (MQ=255)
                                                                                   ggcggcGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAgg                                                             <  1:1177487/150‑1 (MQ=255)
                                                                                      ggcgTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGAt                                                         <  1:816861/151‑1 (MQ=255)
                                                                                      ggcgTGATGATTTACAACACCCACGGTTAGATGAGGCCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGAt                                                         >  2:1368643/1‑151 (MQ=255)
                                                                                           gatgatTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCg                                                    >  1:220726/1‑151 (MQ=255)
                                                                                                 ttACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTc                                              >  1:221982/1‑151 (MQ=255)
                                                                                                  tACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTcc                                             >  1:1317708/1‑151 (MQ=255)
                                                                                                   acaacaCCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCcg                                            >  2:530054/1‑151 (MQ=255)
                                                                                                     aacaCCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCcgcg                                          <  2:597448/151‑1 (MQ=255)
                                                                                                       cacCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCGCTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCcgcgcg                                        >  2:164525/1‑151 (MQ=255)
                                                                                                       cacCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGCGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCcgcgcg                                        >  1:1177149/1‑151 (MQ=255)
                                                                                                           cACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGatat                                    <  1:110188/151‑1 (MQ=255)
                                                                                                           cACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGatat                                    <  2:1324228/151‑1 (MQ=255)
                                                                                                           cACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGatat                                    >  2:651231/1‑151 (MQ=255)
                                                                                                               gTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTc                                >  2:142375/1‑151 (MQ=255)
                                                                                                                  aGATGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCgg                              <  2:562174/150‑1 (MQ=255)
                                                                                                                    aTGAGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCAt                           <  1:476644/151‑1 (MQ=255)
                                                                                                                       aGGTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATgg                         >  1:710231/1‑150 (MQ=255)
                                                                                                                         gTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTGCGCGCGATATCTTCGGCATg                          >  1:857019/1‑147 (MQ=255)
                                                                                                                         gTCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCGTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATg                      <  1:178978/151‑1 (MQ=255)
                                                                                                                          tCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGt                     >  2:338253/1‑151 (MQ=255)
                                                                                                                          tCGCGGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGt                     >  1:811752/1‑151 (MQ=255)
                                                                                                                             cgGAAGAAGTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGt                  <  2:1114044/151‑1 (MQ=255)
                                                                                                                                  gaagTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGttcctt             <  1:618559/151‑1 (MQ=255)
                                                                                                                                    agTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGttccttcc           <  2:41439/151‑1 (MQ=255)
                                                                                                                                    agTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGttccttcc           <  1:655409/151‑1 (MQ=255)
                                                                                                                                     gTTCCCGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTtccttcct          >  1:1098796/1‑151 (MQ=255)
                                                                                                                                        cccGTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGt        >  2:557415/1‑150 (MQ=255)
                                                                                                                                           gTTGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGTAATg    <  2:124114/151‑1 (MQ=255)
                                                                                                                                            ttGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGTAATgg   >  2:1413834/1‑151 (MQ=255)
                                                                                                                                             tGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGTAATgg   >  1:1158399/1‑150 (MQ=255)
                                                                                                                                             tGCACTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGTAATGGt  >  1:72265/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGGTTAAGAGAAAGTGAAAAAAACACAGCGAAAAGAAATCGAAAATGTGACAAATATCACAGGTGTTCGTCAAATTGAGTTATGGCGGCGTGATGATTTACAACACCCACGGTTAGATGAGGTCGCGGAAGAAGTTCCCGTTGCGCTGGTCTACAACGGCATTTCGCATGTGGTGATGATGGCGTCGCCCAAAGACCTTGAGTACTTTGCGCTCGGTTTTTCGCTTTCCGAGGGGATTATCGAAAGTCCGCGCGATATCTTCGGCATGGATGTCGTTCCTTCCTGTAATGGT  >  NC_000913/4086003‑4086294

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