Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,757 T→C 8.1% intergenic (‑140/+24) ylaB ← / ← ylaC putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7570TC 8.1% 349.5 / 17.2 125intergenic (‑140/+24)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (56/56);  new base C (6/4);  total (64/61)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 132 total aligned reads displayed.

GTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGATGG  >  NC_000913/478608‑478905
                                                                                                                                                     |                                                                                                                                                    
gtgtTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCtt                                                                                                                                                      <  2:1092866/150‑1 (MQ=255)
 tgtTCTCACCATCATATCCTTATTAACCCAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGGTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCttt                                                                                                                                                     <  2:709362/150‑1 (MQ=255)
 tgtTCTCACCAGCCTAACCTTATTACCCCAGTGTAGGGAAAGTATATCTGATTGTGGAAACAGGCTGGCTGGAGATAAGGAGGTAGAGCTTGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCttt                                                                                                                                                     <  2:1102121/150‑1 (MQ=255)
 tgtTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCttt                                                                                                                                                     <  2:390564/150‑1 (MQ=255)
  gtTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAATAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                                                                                    <  1:620294/150‑1 (MQ=255)
    tctcACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGCGCTTTACg                                                                                                                                                  >  2:1462770/1‑150 (MQ=255)
    tctcACCAGCATAACCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTTGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCTTTTGACCAGGCGCAGAGCTTTACg                                                                                                                                                  <  1:214375/150‑1 (MQ=255)
       cACCAGCATATCCTTCTTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCt                                                                                                                                               <  2:1684618/150‑1 (MQ=255)
       cACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                                                               <  2:557766/150‑1 (MQ=255)
        aCCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                                                                              <  2:85608/150‑1 (MQ=255)
         ccAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                                                                             <  2:436295/150‑1 (MQ=255)
         ccAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                                                                             <  2:582325/150‑1 (MQ=255)
         ccAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                                                                             >  2:1342415/1‑150 (MQ=255)
          caccATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                                                                            <  2:260285/147‑1 (MQ=255)
           aGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                                                                           <  2:794364/150‑1 (MQ=255)
            gCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTAt                                                                                                                                          >  2:787580/1‑150 (MQ=255)
             cATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATg                                                                                                                                         >  2:1640140/1‑150 (MQ=255)
                  ccTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATTGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                                                                    >  1:1701028/1‑150 (MQ=255)
                  ccTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                                                                    >  1:1025751/1‑150 (MQ=255)
                   cTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                                                                   >  1:1284492/1‑150 (MQ=255)
                   cTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                                                                   >  2:611510/1‑150 (MQ=255)
                   cTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCGGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                                                                   <  2:1056279/150‑1 (MQ=255)
                     tattaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCGGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                                                                 <  2:439445/150‑1 (MQ=255)
                       ttaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACTAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAAc                                                                                                                               <  1:696975/150‑1 (MQ=255)
                         aaCACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACtt                                                                                                                             >  2:952002/1‑150 (MQ=255)
                                 gtAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                                                                     <  2:470881/150‑1 (MQ=255)
                                   aGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGAATGTAACCCTAACTTAAGATGTTGa                                                                                                                   >  1:1103115/1‑150 (MQ=255)
                                    gggAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                                                                  >  1:503540/1‑150 (MQ=255)
                                     ggAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                                 >  1:1167643/1‑150 (MQ=255)
                                     ggAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                                 >  1:355338/1‑150 (MQ=255)
                                      gAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTc                                                                                                                <  1:1181789/150‑1 (MQ=255)
                                         aGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                                                                             >  1:312153/1‑150 (MQ=255)
                                         aGTATAGCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                                                                             <  1:353029/150‑1 (MQ=255)
                                            atatCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                                                                          >  1:1576989/1‑150 (MQ=255)
                                              atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                                                                        >  1:357698/1‑150 (MQ=255)
                                                  gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                                                                    <  1:1640140/150‑1 (MQ=255)
                                                    ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                                                                  <  2:66641/150‑1 (MQ=255)
                                                     tgtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTACCCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                                                                 >  1:333964/1‑150 (MQ=255)
                                                      gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                                                                >  2:1446386/1‑150 (MQ=255)
                                                      gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                                                                >  2:492128/1‑150 (MQ=255)
                                                         gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                             >  2:372704/1‑150 (MQ=255)
                                                           aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                                                                           >  1:287241/1‑150 (MQ=255)
                                                            aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                                                                          >  1:161640/1‑150 (MQ=255)
                                                             aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                                                                         <  2:357698/150‑1 (MQ=255)
                                                               ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                                                                       >  1:1110088/1‑150 (MQ=255)
                                                               ggctggctGGAGATAAGGAGGGAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                                                                       <  2:500415/150‑1 (MQ=255)
                                                               ggctggctGGAGAAAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                                                                       <  1:1569977/150‑1 (MQ=255)
                                                                 ctggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTa                                                                                     <  1:1404962/150‑1 (MQ=255)
                                                                  tggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                                                                    >  1:743096/1‑150 (MQ=255)
                                                                         agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTTGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                                                                             <  1:1342415/150‑1 (MQ=255)
                                                                            tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                                                                          >  2:1133172/1‑150 (MQ=255)
                                                                              aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTcc                                                                        >  2:1287237/1‑149 (MQ=255)
                                                                              aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                                                                        <  2:1592891/150‑1 (MQ=255)
                                                                               ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                                                                       >  2:218750/1‑150 (MQ=255)
                                                                                gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                                                                      <  1:901421/150‑1 (MQ=255)
                                                                                  ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                                                                    <  1:1349927/150‑1 (MQ=255)
                                                                                      gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                                                                <  2:1595543/150‑1 (MQ=255)
                                                                                        gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTAACCTTTAAGGa                                                              >  1:900708/1‑150 (MQ=255)
                                                                                         cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc                                                             >  2:1038111/1‑150 (MQ=255)
                                                                                          tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg                                                            <  1:1588376/150‑1 (MQ=255)
                                                                                           ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa                                                           >  2:978182/1‑150 (MQ=255)
                                                                                           ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa                                                           <  2:1178941/150‑1 (MQ=255)
                                                                                           ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa                                                           >  1:309078/1‑150 (MQ=255)
                                                                                              aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTAACCTTTACGGACGATCa                                                        >  2:1043215/1‑150 (MQ=255)
                                                                                                ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt                                                      >  1:38627/1‑150 (MQ=255)
                                                                                                      taGTGATGGGTGGAAACGATCATTTTAATTTGACCAGGCGACGAGCTTTACGCCTGGTATGTAAACCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGc                                                <  2:1701028/149‑1 (MQ=255)
                                                                                                         tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgt                                             <  1:372704/150‑1 (MQ=255)
                                                                                                               gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGTGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                                       <  2:1107314/150‑1 (MQ=255)
                                                                                                                gggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTTGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCt                                      <  1:222363/149‑1 (MQ=255)
                                                                                                                  gAAACGATCATTGTCAATTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTa                                    <  2:90947/150‑1 (MQ=255)
                                                                                                                   aaaCGATCATTGTCATTTGACCAGGTGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTAt                                   <  1:1006383/150‑1 (MQ=255)
                                                                                                                   aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTAt                                   <  2:216138/150‑1 (MQ=255)
                                                                                                                           aatGTCATTAGACCAGGCGCAGAGCTTTACGCCTGGTATGAAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATcc                           <  2:1097010/148‑1 (MQ=255)
                                                                                                                             tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCcgc                         <  1:787580/150‑1 (MQ=255)
                                                                                                                              gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCcgcg                        <  2:1403803/150‑1 (MQ=255)
                                                                                                                              gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCcgcg                        <  1:902229/150‑1 (MQ=255)
                                                                                                                               tCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGa                       <  1:1678101/150‑1 (MQ=255)
                                                                                                                                ttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac                      >  1:242923/7‑150 (MQ=255)
                                                                                                                                ttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac                      >  1:792558/7‑150 (MQ=255)
                                                                                                                                 aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAATATGTTGATTCGGCGCGAGCTAGGGTAAAAATATCGTAAAAAGACAGTACACCGTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGacc                     >  2:305549/1‑149 (MQ=255)
                                                                                                                                 aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGaca                     >  1:701171/1‑150 (MQ=255)
                                                                                                                                   cataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAGCCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcg                   >  1:373565/4‑150 (MQ=255)
                                                                                                                                     taCCATGAGCAGAGATTTACGCCTGGTATGTAACACTAACTTAAGATGTTGATTCGGCGCTAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg                 <  1:164376/149‑1 (MQ=255)
                                                                                                                                       ccAGGCGTAAAGCTCTGCGCCTGGTATGAAACCCTTACTTAAGATGTTGATTCGGCGAGAGCGAGGGAAAAAATATCTTAAAAAGCCACTACACCTTTACGGACGAGCATTATTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGga               >  2:213688/1‑150 (MQ=255)
                                                                                                                                       ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGga               >  2:743010/1‑150 (MQ=255)
                                                                                                                                       ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAATCGGCGCGAGCGAGGGTAAAAATATCTTATAAAGACAGTTCACCTTTACGGACGATCATTTTTTTCAGTTGTTCCTTATGAACATCCGCGACACGCGga               >  2:241062/1‑150 (MQ=255)
                                                                                                                                       ccAGGCGCAGAGCTGTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGCGCGAGGGGAAAAATATCGTAAAAAGACAGGTCACCGTTACGGACGCTCATTTTTTGCAGGTGTTCCTTATGAACATCCGCGACACGCGga               >  1:226332/1‑150 (MQ=255)
                                                                                                                                        cAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgag              >  2:1260273/1‑150 (MQ=255)
                                                                                                                                        cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgag              <  1:188601/150‑1 (MQ=255)
                                                                                                                                        cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgag              >  2:697068/1‑150 (MQ=255)
                                                                                                                                        cAGGCGCAGAGCTGTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgag              <  2:1160092/150‑1 (MQ=255)
                                                                                                                                         aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgaga             >  1:1597098/1‑150 (MQ=255)
                                                                                                                                          ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGagag            >  1:66391/1‑150 (MQ=255)
                                                                                                                                          ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGagag            <  1:1088722/150‑1 (MQ=255)
                                                                                                                                          ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGagag            >  2:840845/1‑150 (MQ=255)
                                                                                                                                           gcgcAGAGCTTTACGCCTGGTATGAAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGTGTCGCGgatgtt           <  1:611510/150‑5 (MQ=255)
                                                                                                                                              cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCa        <  2:1693131/150‑1 (MQ=255)
                                                                                                                                                  gCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGAt    >  2:743854/1‑150 (MQ=255)
                                                                                                                                                   cTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGATCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGATg   <  2:355338/150‑1 (MQ=255)
                                                                                                                                                    tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGAGTGGGTAAAAATATCGTATAAAGACTGTTCAACTTTACGTACGATCATTTTTTGAAGTTGTTCCTTATGAACATAAGCTACACGCGGAGAGTTGAAGATgg  >  2:236869/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTCAAGATGG  >  NC_000913/478608‑478905

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