Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,515,731 C→A 5.8% E700D (GAG→GAT yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,515,7310CA 5.8% 189.1 / 3.5 70E700D (GAG→GATyfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (39/26);  new base A (2/2);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

TGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGC  >  NC_000913/2515585‑2515872
                                                                                                                                                  |                                                                                                                                             
ttgtccAACCCTGAAGTTAGCTCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGtctc                                                                                                                                             <  1:559471/144‑1 (MQ=255)
  gtccAACCCTTAAGTGAGCGCTTATGGGATAATTCCCCGGTTTTGTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGTTATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGa                                                                                                                                           <  1:82131/146‑1 (MQ=255)
   ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAACTACCCTTGGCGGTATTGAAACCCGGGCCTACGCGATACCGCCGACTGCTGTTGCGtcttgtc                                                                                                                                          >  1:283236/1‑145 (MQ=255)
   ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGAc                                                                                                                                          >  1:901208/1‑149 (MQ=255)
    tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTCGGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACa                                                                                                                                         <  1:775697/149‑1 (MQ=255)
    tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAAATATTGCACCCCGAGCTTATGCACTAGCGCCTGCTGCTGTTGCGTCTCGACa                                                                                                                                         >  2:220685/1‑149 (MQ=255)
          cTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCCATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTcc                                                                                                                                   >  2:216401/1‑149 (MQ=255)
          cTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTcc                                                                                                                                   >  2:185429/1‑149 (MQ=255)
            tAAGTTAGCGGTTATGGGATAATTCCCCGGTTATTTTACGCCTGGTAATCGCCGAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATGGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCcgc                                                                                                                                 <  2:1150332/149‑1 (MQ=255)
              aGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGa                                                                                                                               <  2:91657/149‑1 (MQ=255)
              aGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGa                                                                                                                               <  1:772828/149‑1 (MQ=255)
              aGTTAGCGCTTATGGGATAATGCCCCGGGTTTTTTACGCCTGTTAATCAGCCAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGa                                                                                                                               <  1:850755/149‑1 (MQ=255)
               gTTAGCGCTTATGGGATAATGCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGAc                                                                                                                              <  1:52055/149‑1 (MQ=255)
                ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTTCTTTTTTTAATCAGCTAATGGCTGCGGGCGAACAATCAAATACCCTTGCAGATATTGCACCCCGTGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACCCCGCGAcc                                                                                                                             >  2:864061/1‑149 (MQ=255)
                ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGCAAAATAACATACCATTGAAGATATTGCACCCCGAGCTGAAGCAATAACGCCTGCTACTGTTGCGTCTTTTCAATCTCCGCCAcc                                                                                                                             >  2:1274458/1‑149 (MQ=255)
                ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGAcc                                                                                                                             >  2:1004264/1‑149 (MQ=255)
                ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACACTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGAcc                                                                                                                             >  2:1156430/1‑149 (MQ=255)
                 tAGCGCTTATGGGATAATTCCCCGGATTTTTTAAGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCACATACCCTTGCAGATATTGCACCCCGACATCATGCAATACCGCCTGCTCCTGTTGCCTCCCGACAAACACcgccaaca                                                                                                                            >  2:551893/1‑144 (MQ=255)
                  aGCGCTTATGGGATAATTCCCCGGTTTGTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCGAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACcaa                                                                                                                           >  2:1023875/1‑148 (MQ=255)
                      cTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTc                                                                                                                       >  2:1313264/1‑149 (MQ=255)
                         aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGAGCCTGCTGCTGTTACGTCGCGACAAACTCCGCGACCACACTCAAt                                                                                                                    >  2:740599/1‑149 (MQ=255)
                              aTAATTCCCCGGTTTTTTTACGCCTGGTAATCAGCTAATGGCTGTGGGCGACCAATCAAATATCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCTCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACtt                                                                                                               <  2:334013/149‑1 (MQ=255)
                                  ttCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGGGCGTATGCAATAGCGCCTGCTGCTGTTGCGTCGCGACAAAATCCGCGACCACACTCAATGCCtgctgt                                                                                                           >  1:441207/1‑144 (MQ=255)
                                     cccGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCCAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTCCGCTTTc                                                                                                        >  1:1180552/1‑149 (MQ=255)
                                     cccGCTTGTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATCCCCTTGCAGATCTTGCACCCCGAGCTTATCCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTc                                                                                                        <  2:136129/149‑1 (MQ=255)
                                       cGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGc                                                                                                      <  2:312961/149‑1 (MQ=255)
                                         gTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCAACCCGAGCTTATGCAATAGCGCCTGCTGTTGTTGCTGCTCGACAAAGTCCGAGACCACACGCAATGACTTCGCTTTCGCCa                                                                                                    >  1:38624/1‑149 (MQ=255)
                                          tttttttACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAg                                                                                                   >  2:297607/1‑149 (MQ=255)
                                            tttttACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACCTCGCTTTCGCCAGAt                                                                                                 <  2:457408/149‑1 (MQ=255)
                                            tttttACGCCTGTTAATCAGCTAAGGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCCATGACTTCGCTTTCGCCAGAt                                                                                                 <  1:1207082/149‑1 (MQ=255)
                                                      tGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTTCTGACATTCTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGa                                                                                       >  1:647051/1‑149 (MQ=255)
                                                           aTCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCTCCTGCTGCTGTTGCGTATCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTca                                                                                  >  2:1249014/1‑149 (MQ=255)
                                                           aTCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTGCGCCAGATCGGTAATTGAtctcc                                                                                  >  1:9358/1‑148 (MQ=255)
                                                           aTCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTATCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTca                                                                                  >  2:133707/1‑149 (MQ=255)
                                                            tCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTcac                                                                                 >  1:280667/1‑149 (MQ=255)
                                                             cAGCTAATGGCTGCGGGCGACCCATCAAGGACCCTTGCAGATATTGCACCCCGAGCTTATACAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTccca                                                                                >  2:963146/1‑149 (MQ=255)
                                                              aGCTAATGGCTGCGGGCGACCAATCAACTACCCGTTCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGAtcttttaa                                                                               <  2:283236/149‑6 (MQ=255)
                                                                       cTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCCTCGCAt                                                                      >  2:303468/1‑149 (MQ=255)
                                                                       cTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCAt                                                                      >  1:1264985/1‑149 (MQ=255)
                                                                        tGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATc                                                                     <  1:189937/149‑1 (MQ=255)
                                                                          cGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGAACGGTAATTGATCTCACAATCATAGCATCCa                                                                   >  1:331699/1‑149 (MQ=255)
                                                                             gCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCg                                                                >  1:1256612/1‑149 (MQ=255)
                                                                             gCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCg                                                                <  1:1282019/149‑1 (MQ=255)
                                                                                   ttttaaaTACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCCGGTTCg                                                          >  2:838241/5‑149 (MQ=255)
                                                                                   aaTCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCg                                                          >  1:1094358/1‑149 (MQ=255)
                                                                                     tCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTCATTGATCTCACAATCATCGCATCCAGCGAGGTCGtc                                                        >  2:680904/1‑148 (MQ=255)
                                                                                     tCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGGAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTg                                                        >  2:851366/1‑149 (MQ=255)
                                                                                        aaTACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGCCa                                                     >  2:1118804/1‑149 (MQ=255)
                                                                                            cccTTGCAGATTTTGCACCCCGATCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAACCTCCGCGACCACACTCAATGTCCTCGCTTTCGCCTGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAatat                                                 <  1:1038910/149‑1 (MQ=255)
                                                                                               ttGCAGATATTGCACCCCGAGCTTATGCAATAGCGGCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCCCACTCATCCCATCCAGCGTGTTCGTGACAATATCtc                                              >  2:815083/1‑148 (MQ=255)
                                                                                                tGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCttt                                             <  1:713137/149‑1 (MQ=255)
                                                                                                tGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCttt                                             <  1:297607/149‑1 (MQ=255)
                                                                                                tGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCttt                                             >  1:1311630/1‑149 (MQ=255)
                                                                                                tGCAGATATTCCACCCCGATCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAACCTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGCAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCttt                                             <  1:157744/149‑1 (MQ=255)
                                                                                                    gATATTGCACCCCGAGCTTATGCAATAGTGCCTGCTGCTGTTGCGTCTCGGCAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCTAGATCGGTTTTTTATCTCACAATCATCGCATCCAGCGTGTTTGTGTCAATTTTTTTcacg                                         >  2:675795/1‑148 (MQ=255)
                                                                                                       atTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATGGGTAATTGATCTCACAATCATCGCTTCCAGCGTGTTCGTTACAATATCTTTCCCCAAg                                      >  1:822/1‑149 (MQ=255)
                                                                                                         tGCACCCCGAGCTTATGCAATAGCTCCTGCTGCTGTTGCGTATCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGAc                                    <  1:1249014/149‑1 (MQ=255)
                                                                                                         tGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTATCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGAc                                    <  1:133707/149‑1 (MQ=255)
                                                                                                          gCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACg                                   >  2:443633/1‑149 (MQ=255)
                                                                                                                    ttgagtcAATAGCGGCTGCTGCGGTTGCGTCTCGACACACGTCGCGACCACACTCACTGACTTCGCTTTCGCCAGATCGGTAATTCATCTCACCATCATCTCAGCCAGCGTGTTTGTGACAATATCTTTCACAAAGACGCCATCAAttt                         <  2:136454/143‑1 (MQ=255)
                                                                                                                     ttATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGCCTTCTCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGGGTTCGTGGCAATATCTTTCACAAAGCCGCCCTCAAtttt                        >  1:737698/1‑149 (MQ=255)
                                                                                                                      tATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTtg                       >  1:577820/1‑149 (MQ=255)
                                                                                                                      tATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTAAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCACCAATTTtg                       >  1:599683/1‑149 (MQ=255)
                                                                                                                       aTGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTtga                      >  2:1289324/1‑149 (MQ=255)
                                                                                                                       aTGCAACAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTtga                      <  1:423320/149‑1 (MQ=255)
                                                                                                                           aaTAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTgatgat                  >  1:208376/1‑149 (MQ=255)
                                                                                                                           aaTAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTgatgat                  <  1:220685/149‑1 (MQ=255)
                                                                                                                                 gcCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGc            <  1:351344/149‑1 (MQ=255)
                                                                                                                                 gcCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCCTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGc            <  1:342573/149‑1 (MQ=255)
                                                                                                                                 gcCGGCTGCTGGTGCGGCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGc            <  2:647051/149‑1 (MQ=255)
                                                                                                                                    tgctgctgTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGACCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTg         >  1:1050753/1‑149 (MQ=255)
                                                                                                                                           gTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGc  <  1:563970/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGC  >  NC_000913/2515585‑2515872

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