Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,027,173 T→C 5.4% Y282C (TAT→TGT)  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,1730TC 5.4% 230.6 / 2.9 74Y282C (TAT→TGT) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (38/32);  new base C (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAG  >  NC_000913/2027026‑2027317
                                                                                                                                                   |                                                                                                                                                
cAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATa                                                                                                                                                 <  2:347006/149‑1 (MQ=255)
cAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATa                                                                                                                                                 <  1:189835/149‑1 (MQ=255)
     aaCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCAt                                                                                                                                                  <  1:166771/143‑1 (MQ=255)
     aaCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCAt                                                                                                                                                  >  2:166771/1‑143 (MQ=255)
      aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGc                                                                                                                                             <  6:373413/147‑1 (MQ=255)
      aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCt                                                                                                                                            <  3:123824/148‑1 (MQ=255)
             gACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCa                                                                                                                                    >  1:314210/1‑149 (MQ=255)
                tGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCa                                                                                                                                 >  6:199337/1‑149 (MQ=255)
                  tCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATc                                                                                                                               >  7:333677/1‑149 (MQ=255)
                     aTACGAATGCCGCCCCGCGAGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAg                                                                                                                            >  3:219764/1‑149 (MQ=255)
                        cGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGt                                                                                                                          <  7:175586/148‑1 (MQ=255)
                                         gtCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTc                                                                                                                                     <  3:166865/120‑1 (MQ=255)
                                         gtCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTc                                                                                                                                     >  4:166865/1‑120 (MQ=255)
                                            aTGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCtt                                                                                                      >  8:299035/1‑148 (MQ=255)
                                             tGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTc                                                                                                                                     >  4:166161/1‑116 (MQ=255)
                                               ttCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTaaa                                                                                                   >  7:137684/1‑148 (MQ=255)
                                                    ttGCCTGCGCCAGCAAGGCCAGTTCACGATGATCAAAAATATTCCCTGTTTGATTGGCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgc                                                                                             <  4:184353/149‑1 (MQ=255)
                                                    ttGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgc                                                                                             >  8:205269/1‑149 (MQ=255)
                                                    ttGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgc                                                                                             <  8:137684/149‑1 (MQ=255)
                                                    ttGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgc                                                                                             >  7:157708/1‑149 (MQ=255)
                                                      gCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgct                                                                                            <  7:299035/148‑1 (MQ=255)
                                                        cTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCa                                                                                         >  1:265048/1‑149 (MQ=255)
                                                        cTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCa                                                                                         <  3:244441/149‑1 (MQ=255)
                                                        cTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCa                                                                                         >  6:142664/1‑149 (MQ=255)
                                                           cgcCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATGTCAGCTTACTGTCATAGAGCTGATACTCACCATTGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATa                                                                                       <  5:199337/148‑1 (MQ=255)
                                                             ccAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTa                                                                                    >  6:145292/1‑149 (MQ=255)
                                                             ccAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTa                                                                                    <  1:313044/149‑1 (MQ=255)
                                                                 cAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTCCGCa                                                                                >  3:83852/1‑149 (MQ=255)
                                                                    ggCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATc                                                                                                                   <  7:175576/111‑1 (MQ=255)
                                                                    ggCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATc                                                                                                                   >  8:175576/1‑111 (MQ=255)
                                                                    ggCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACgg                                                                             <  4:337591/149‑1 (MQ=255)
                                                                       cAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGaa                                                                          >  7:108171/1‑149 (MQ=255)
                                                                           tCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCACAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACt                                                                      >  3:278552/1‑149 (MQ=255)
                                                                           tCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCACAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACt                                                                      >  3:283497/1‑149 (MQ=255)
                                                                            cACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACt                                                                      >  3:139230/1‑148 (MQ=255)
                                                                             aCGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTc                                                                     <  8:333677/148‑1 (MQ=255)
                                                                             aCGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTc                                                                     <  8:165516/148‑1 (MQ=255)
                                                                             aCGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCa                                                                    >  8:217072/1‑149 (MQ=255)
                                                                               gAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATc                                                                  >  8:303345/1‑149 (MQ=255)
                                                                                   aTCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAATATCTCAGCTTACTGTGATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCa                                                              <  8:235152/149‑1 (MQ=255)
                                                                                    tCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATTCTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgct                                                                                            >  7:82724/1‑118 (MQ=255)
                                                                                    tCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATTCTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgct                                                                                            <  8:82724/118‑1 (MQ=255)
                                                                                      aaaaaTATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCa                                                           >  1:187239/1‑149 (MQ=255)
                                                                                      aaaaaTATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCa                                                           >  3:314317/1‑149 (MQ=255)
                                                                                            atTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCa                                                      >  2:57672/1‑148 (MQ=255)
                                                                                               cccTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATAc                                                   <  7:217072/148‑1 (MQ=255)
                                                                                                ccTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATAc                                                   >  5:338510/1‑147 (MQ=255)
                                                                                                ccTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATAc                                                   >  7:382123/1‑147 (MQ=255)
                                                                                                ccTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCACAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCa                                                 <  4:283497/149‑1 (MQ=255)
                                                                                                ccTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCACAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCa                                                 <  4:278552/149‑1 (MQ=255)
                                                                                                acTGTTGGATGGTCAGGATTGGAAGAGGGCAAAAATCTCAGCTTACTGTAATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCa                                                 <  2:265048/148‑1 (MQ=255)
                                                                                                 cTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCag                                                >  8:183954/1‑149 (MQ=255)
                                                                                                       gATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGtt                                           <  1:161832/148‑1 (MQ=255)
                                                                                                        aTGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGttat                                         <  8:328624/149‑1 (MQ=255)
                                                                                                           gTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGttatta                                       <  2:185049/148‑1 (MQ=255)
                                                                                                           gTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGttatta                                       <  7:183954/148‑1 (MQ=255)
                                                                                                           gTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTAc                                      >  4:367674/1‑149 (MQ=255)
                                                                                                           gTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGCAATACTCATCCCCAGCACGCCATACCAGTAGTTATTAc                                      >  1:393337/1‑149 (MQ=255)
                                                                                                            tCAGGATTTGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGGTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACt                                     <  7:103590/149‑1 (MQ=255)
                                                                                                            tCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCCTAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACt                                     <  5:142664/149‑1 (MQ=255)
                                                                                                             cAGGATTGGAAGAGGTCAAAAATCTCATCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGATTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTg                                    >  4:10505/1‑149 (MQ=255)
                                                                                                                     taaGAGTTCAAAATTCTCATCTTACTGTCATAGAGCTGATACTGACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCTCGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCg                             <  2:314210/147‑1 (MQ=255)
                                                                                                                     gAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCg                             >  4:177466/1‑148 (MQ=255)
                                                                                                                      aaGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCgg                            >  4:296577/1‑148 (MQ=255)
                                                                                                                         agGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAaca                        <  2:393337/149‑1 (MQ=255)
                                                                                                                              aaaaaTCTCAGCTTACTGTCATAGAGCTGATACTGACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCACAGCACGCCATACCAGTAGTTATTACTGTCTAc                               <  3:281206/137‑1 (MQ=255)
                                                                                                                              aaaaaTCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACTGTAATAATCATCCACAGCAAGCCATACCAGTAGTTATTACTGTCAAc                               >  4:281206/1‑137 (MQ=255)
                                                                                                                              aaaaaTCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGAc                   >  4:9819/1‑149 (MQ=255)
                                                                                                                               aaaaTCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACg                  <  8:65691/149‑1 (MQ=255)
                                                                                                                                       aGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTg          >  3:169780/1‑149 (MQ=255)
                                                                                                                                       aGCTTACTGTCATAGAGCTGATACTCACAATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGTAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTg          <  5:145292/149‑1 (MQ=255)
                                                                                                                                             cTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACTGTAACCTGCGGtt     <  6:338510/148‑1 (MQ=255)
                                                                                                                                              tGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGtt     >  5:100627/1‑147 (MQ=255)
                                                                                                                                              tGTCATAGAGCTGATACTCACCATAGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATACCCAGAACGCCATACAATTAGTTATTACAGTCTACCGGAAAACTGACGTTAACCTGCGGTTCa   >  4:298465/1‑149 (MQ=255)
                                                                                                                                                tCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAg  >  1:122502/1‑148 (MQ=255)
                                                                                                                                                tCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAg  >  1:103288/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAG  >  NC_000913/2027026‑2027317

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