Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,027,258 G→A 6.1% L254L (CTG→TTG)  yedQ ← putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,027,2580GA 6.1% 175.7 / 2.9 69L254L (CTG→TTG) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (35/28);  new base A (2/2);  total (37/32)
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.26e-01

TATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAG  >  NC_000913/2027117‑2027400
                                                                                                                                             |                                                                                                                                              
tatTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCa                                                                                                                                         <  1:814288/149‑1 (MQ=255)
     ccTGTTGGACGGTCAGGCTTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCa                                                                                                                                    <  1:169895/149‑1 (MQ=255)
       tGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCAAGCACGCCATACCagt                                                                                                                                  >  1:484270/1‑149 (MQ=255)
             aTGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCCTCCCCAGCACGCCATACCAGTAGttat                                                                                                                            >  1:349337/1‑149 (MQ=255)
               ggctaGGATTCGAAGAGGTTAAAAATCTCAGCTTACTGTCATAGGGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGttatta                                                                                                                          <  1:502438/145‑1 (MQ=255)
                gTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTAc                                                                                                                         >  2:663574/1‑149 (MQ=255)
                   aGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAAGGTATTATCGATGGCGTTTCTTAAAAATTGCTGCATCTTACGCACGCCAATACTCATCCCCAGCACGCCCTACCCGACGTTATTACTGt                                                                                                                      >  2:275067/1‑149 (MQ=255)
                           aaGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGa                                                                                                              <  1:512140/149‑1 (MQ=255)
                                   aaaaaTCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTGTTATCGATGGCGTTGCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGAc                                                                                                      >  1:55536/1‑149 (MQ=255)
                                    aaaaTCGTAGCTTGCTGTCATAGAGCTGATACGCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGGACACTGTCg                                                                                                     <  2:224813/149‑1 (MQ=255)
                                     aaaTCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCCTAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACgg                                                                                                    >  1:515078/1‑149 (MQ=255)
                                                   tgacctAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCa                                                                                      <  1:2795/144‑1 (MQ=255)
                                                    gTCATAGAGCTGGTACTCACCACCGATGCTTTTAGGGATGTCGTGCGGTCTAAAGTGCTGCATAGCACGCACGGGACTACTTATCCCCAGCACGCCATACCAGTAGTCATTACAGTCTACCGCAACCCTGACGGTAAACTGCGGTTCAg                                                                                     <  2:593153/149‑1 (MQ=255)
                                                     tCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGt                                                                                    >  2:206788/1‑149 (MQ=255)
                                                     tCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTATTAAAAATTGCTGCATAGTACGCACGGGTATAGTCATCCCCAGCACGCCATACCCGTAGTTATTACTGTCTTCCGGAACACTGACGGTAACCTGCgggtactg                                                                                    >  2:537368/1‑143 (MQ=255)
                                                       aTAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTAt                                                                                  >  2:355953/1‑149 (MQ=255)
                                                       aTAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAACACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTAt                                                                                  >  1:1135679/1‑149 (MQ=255)
                                                       aTAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAACACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTAt                                                                                  >  1:492819/1‑149 (MQ=255)
                                                        tAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTCGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCatgatt                                                                                 >  2:797591/1‑144 (MQ=255)
                                                        tAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGGAtt                                                                                 >  2:1027772/1‑149 (MQ=255)
                                                                    tccccctCGGGGGTTTTATCGATGGCGTTTCTTAAAAACTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTc                                                                     <  2:613883/143‑1 (MQ=255)
                                                                         tgcggggTTTTTATGGAGGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGtt                                                                <  2:250850/144‑1 (MQ=255)
                                                                         aTCGCGGTTTTCATCGATGGCGTTTCTTACAAATTGCTCCATAGTACGCACGGGAATACTCAGCCCCAGCACGCACTACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTCTCCCgggt                                                                <  2:1200476/149‑3 (MQ=255)
                                                                          tCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTg                                                               >  1:533917/1‑149 (MQ=255)
                                                                          gcGAGGTCTGTAGCGATGGCGTTTCTTAAAAATTGCTCCATAGTACGCACGGGAATACTCAACCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTg                                                               <  1:157114/148‑1 (MQ=255)
                                                                            gAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCg                                                             >  2:362031/1‑149 (MQ=255)
                                                                                tttttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGt                                                         <  1:587241/149‑1 (MQ=255)
                                                                                 ttttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTg                                                        >  1:1110270/1‑149 (MQ=255)
                                                                                 ttttATCGATGGCGTTTCTCAAAAATGGCTGCATAGTACGCACGGGAATACTCATCCCCAGCCCGCCATACCAGAAGTTTTTACTTTCTTCCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTg                                                        <  2:334343/149‑1 (MQ=255)
                                                                                  ttcaTCGATGCCGTTTCTTAAAAATCGCTGCACAGCGCGCAGTGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGCCGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGa                                                       <  1:31058/146‑1 (MQ=255)
                                                                                   ttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCCGTCGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGGTCCGTATTGCTGGCGTGTTCCGGTTGCGAAGTGaa                                                      >  2:587315/1‑149 (MQ=255)
                                                                                   ttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGaa                                                      >  1:423949/1‑149 (MQ=255)
                                                                                   tcaTCGATGGCGTTTCTTAAAAATTGCTGCATACTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGaa                                                      <  1:114201/147‑1 (MQ=255)
                                                                                      tCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCa                                                   <  1:1136300/149‑1 (MQ=255)
                                                                                      tCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGGATTCCGGGCGTGTACCGGTCACGAAGTGCAcga                                                   >  2:237255/1‑147 (MQ=255)
                                                                                       ccaTGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCAGCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAg                                                  <  2:264934/147‑1 (MQ=255)
                                                                                       cGATGTCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAg                                                  <  1:1174291/149‑1 (MQ=255)
                                                                                       cGATGGCGTTTCGTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAg                                                  <  2:937053/149‑1 (MQ=255)
                                                                                        gATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGc                                                 >  1:494733/1‑149 (MQ=255)
                                                                                        gATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGAGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAGCCTGCGGAGCAGTATTGCTGGCGTGTTCCGGTTGCGGAGTGAagcagc                                                 >  1:388775/1‑149 (MQ=255)
                                                                                        gATGGCGGTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGATCCAGTTGCGAAATGAACCAGc                                                 >  1:19112/1‑149 (MQ=255)
                                                                                          tgccGTTTTTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGt                                               <  2:386884/146‑1 (MQ=255)
                                                                                            gcttttCTTAAAAATTGCTGCATAGTACGCACGGGAATACGCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTAc                                             <  1:275067/146‑1 (MQ=255)
                                                                                               tttCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGtccggg                                          >  1:428027/1‑144 (MQ=255)
                                                                                               tttCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCAGTTGCGAAATCACCCAGCGTACcgc                                          >  1:1133778/1‑149 (MQ=255)
                                                                                                ttCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGGGTGTTCCGGTTGCGAAGTGAACCAGCGTACcgcg                                         >  1:243991/1‑149 (MQ=255)
                                                                                                   ttAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCCGCGTCCcgcgcgc                                      >  2:564040/1‑148 (MQ=255)
                                                                                                   ttAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCCCCCgg                                      >  2:27803/1‑149 (MQ=255)
                                                                                                     aaaaaTTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTg                                    >  2:814714/1‑149 (MQ=255)
                                                                                                      aaaaTTGCTGCATATTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGGTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGCTTGCGAAGTGAACCAGCGTACCGCGCGGTgc                                   >  2:188015/1‑148 (MQ=255)
                                                                                                          ttGCCGCATAGTACGCATGGGATTACTCATCCCCAGCACTACATACCAGTAGTTATTACTGTCGACCGGAACTCTTACGGTAACCTTCGGTTCAGTATGGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACggt                               <  2:500238/149‑3 (MQ=255)
                                                                                                           tgctgcATAGTACGCACGAGAATACTCATAAACCGCCTGAGAGGCCAGTAGTTATTACTCTCTACCGGAACACTGACGGTACCCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGCGCACAAGCGTACCGCCCGGTGACgggt                              >  1:314216/1‑146 (MQ=255)
                                                                                                             ctgcATAGTACGCATGGGAATACTCATCCCCAGCACTCCATATCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGAtttt                            <  1:915628/149‑1 (MQ=255)
                                                                                                             ctgcATAGTACGCACGCGAATACACCTCCCCAGCACGCCATACCAGTAGTTATTACTGTCTGCCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGAtttt                            <  2:1152725/149‑1 (MQ=255)
                                                                                                               gcATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTtct                          <  1:355953/149‑1 (MQ=255)
                                                                                                               gcATAGTACGCACGGGAATACTCATCCCCAACACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTtct                          <  2:492819/149‑1 (MQ=255)
                                                                                                               gcATAGTACGCACGGGAATACTCATCCCCAACACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTtct                          <  2:1135679/149‑1 (MQ=255)
                                                                                                                       cGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGCTTTTCTCGTtgtta                  >  1:159617/1‑147 (MQ=255)
                                                                                                                        gCACGGGTATACTCAGCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGaa                 >  2:1087028/1‑149 (MQ=255)
                                                                                                                        gCACGGGAATACTCCTCACCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGaa                 <  2:377933/149‑1 (MQ=255)
                                                                                                                        gCACGGGAATACTCAGCCCCAGCACGCGATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGaa                 <  1:1206044/149‑1 (MQ=255)
                                                                                                                             ggAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGcc            >  1:359695/1‑149 (MQ=255)
                                                                                                                             ggAATACTCATCCCCACCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTAATGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGcc            <  1:40660/149‑1 (MQ=255)
                                                                                                                                  aCTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATaa       >  1:405091/1‑149 (MQ=255)
                                                                                                                                    tCATCCCCAGCACGCCATACCAGTATTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGTGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAAc     >  2:660225/1‑149 (MQ=255)
                                                                                                                                      aTCCCAAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCa   <  1:362031/149‑1 (MQ=255)
                                                                                                                                       tCCCCATCACGCCATTCCAGTTGTTATTACTGGCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTACGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGTCAATAAACCAg  <  1:576045/149‑1 (MQ=255)
                                                                                                                                       tCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAAccgg  >  1:23728/1‑147 (MQ=255)
                                                                                                                                       tCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGATAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCac  >  1:663901/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAG  >  NC_000913/2027117‑2027400

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